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ABSTRACT 


The  cascade  algorithm  that  is  used  for  extended  surface  analysis  depends  on  a 
new  parameterization  called  the  thermal  transmission  matrix  to  represent  a  single  fin. 
This  thermal  transmission  matrix,  which  is  intended  to  replace  the  more  familiar  fin 
efficiency  as  a  design  and  analysis  parameterization,  is  a  linear  transformation  that 
maps  conditions  of  heat  flow  and  temperature  at  the  fin  tip  to  heat  flow  and 
temperature  conditions  at  the  fin  base.  The  cascade  algorithm  was  derived  by 
resorting  to  an  analogy  between  a  fin  and  the  electrical  transmission  line.  The 
cascade  algorithm  permits  a  fin  to  be  subdivided  into  many  subfins  each  having  a 
thermal  transmission  matrix  and  then  the  individual  transmission  matrices  for  each  of 
the  subfins  can  be  used,  via  a  simple  matrix  product  to  form  an  overall  equivalent 
thermal  transmission  matrix  for  the  entire  fin.  This  thesis  develops  a  thermal 
transmission  matrix  for  the  radiating  rectangular,  trapezoidal,  and  triangular  fins  both 
for  the  free  space  and  non-free  space  environments.  Test  cases  have  been  run  and 
their  solutions  exactly  match  those  contained  in  the  literature.  The  thesis  concludes 
with  optimization  studies  for  each  profile  considered  where  it  is  observed  that  simple 
algebraic  equations  can  be  employed  to  describe  the  optimum  geometry. 
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I.  INTRODUCTION 


Since  the  advent  of  the  space  age,  extended  surfaces  have  been  extensively 
studied  as  a  means  for  heat  rejection  from  all  sorts  of  vehicles.  Because  weight 
optimization  is  of  the  utmost  importance  in  spacecraft  design,  it  is  desirable  to  have 
an  accurate  technique  that  will  rapidly  optimize  the  geometry  of  the  extended  surface 
used  for  rejecting  the  heat  absorbed  on  or  generated  within  the  vehicle.  To  this  point 
in  time,  the  optimization  problem  involved  the  time  consuming  solutions  of  second 
order  non-linear  differential  equations,  by  numerical  methods  or  by  cumbersome  and 
possibly  inaccurate  graphical  methods.  A  more  advantageous  method  for  determining 
the  minimum  mass,  maximum  radiant  heat  rejection  extended  surface  geometry  has 
been  needed. 

The  purpose  of  this  study  is  to  employ  the  new  parameterization  called  the 
thermal  transmission  matrix  to  represent  a  single  fin  to  evaluate  the  performance  of 
several  longitudinal  fins.  Each  fin  studied  is  subdivided  into  100  subfins  and  the 
thermal  transmission  matrix  is  calculated  for  each  subfin.  Then,  by  a  simple  matrix 
multiplication  of  the  individual  subfin  thermal  transmission  matrices,  the  overall 
transmission  matrix  is  computed.  Not  only  will  this  method  produce  the  base  and  tip 
conditions  of  the  overall  fin  entity,  but  the  temperature  and  heat  flow  at  any  point 
between  the  fin  base  and  fin  tip  can  be  determined.  This  iterative  method  considers 
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the  spacecraft  radiant  heat  rejection  problem  in  both  the  free  and  non-free  space 
environments.  Optimization  must  be  considered  and  simple  algebraic  equations  that 
describe  the  fin  geometry  are  formulated  for  the  longitudinal  fins  of  rectangular, 
trapezoidal,  and  triangular  profiles. 

This  new  approach  to  fin  synthesis  and  analysis  permits  an  exact  solution  to  any 
degree  of  accuracy  that  the  user  wishes  to  specify  and  produces  the  result  in  a 
computational  efficient  manner. 
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n.  THEORETICAL  PRESENTATION 


A.  GENERALIZED  DIFFERENTIAL  EQUATION 

There  are  three  distinct  modes  by  which  heat  flows  from  one  point  to  another. 

A  rather  simplistic  definition  of  each  can  be  stated.  Conduction  is  the  mode  by  which 
heat  flows  through  a  material  by  molecular  motion.  If  the  heat  flows  from  a 
confining  surface  to  a  moving  fluid,  the  mode  is  known  as  convection.  If  no 
intervening  medium  is  present  and  the  thermal  energy  is  transported  from  one  point  to 
another  in  the  form  of  electromagnetic  waves,  the  mode  is  known  as  radiation  or 
radiant  heat  transfer. 

The  flow  of  heat  by  conduction  is  proportional  to  the  temperature  gradient  and 
cross  sectional  area  perpendicular  to  the  direction  of  heat  flow. 


-A(x) 


dm 

dx 


(2.1) 


Insertion  of  a  proportionality  constant  yields 


*  A(x)^  (2.2) 

ax 

where  the  minus  sign  assures  a  positive  heat  flow  in  the  presence  of  a  temperature 
gradient  which  must  be  negative.  This  serves  to  define  the  thermal  conductivity  of 
the  material  written  as  k(T)  because  it  is  a  material  property  almost  always  a  function 
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KT)  = 


(2.3) 


-A(x) 


dUx) 


of  temperature.  Thus,  the  heat  flow  at  x  normal  to  a  surface  of  area  A(x)  at  a 
temperature  T(x)  is  written  as 


^cd  =  ~  Aix) 


dJlx) 


In  convection  the  heat  flow  from  a  confining  surface  at  temperature  T(x)  to  a 
fluid  whose  bulk  temperature  is  will  be  proportional  to  the  temperature  difference 
and  the  amount  of  surface  area. 


«  S(x)  [  J{x)  -  r,  ] 


Insertion  of  a  proportionality  constant  permits  the  representation 


9„  -  K  Sir)  ( m  -  n  ] 


and  this  serves  to  define  the  convective  heat  transfer  coefficient  h,v  as 


^  _ 

^  "  S{x)  [  T(x)  -  r,  ] 
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and  thus  Newton’s  Law  of  Cooling  is  written  as 


=  *cv(^>7)  Six)  [  Tix)  -  r,  ]  (2.8) 

where  the  convective  heat  transfer  coefficient  h„  is  written  as  a  function  of  position 
and  temperature. 

Radiation  requires  no  intervening  medium  and,  for  a  perfect  emitter  and 
absorber  in  full  view  of  one  another,  the  familiar  Stefan-Boltzmann  law  shows  that  the 
heat  flow  is  proportional  to  the  surface  area  and  the  difference  of  the  fourth  powers  of 
the  absolute  temperature. 

=  a  Six)  [  r/  -  T*  ]  (2-9) 

The  proportionality  constant  used  here  is  the  Stefan-Boltzmann  constant, 
a  =  5.66961  x  10*  watt  /  m^  K^ 

If  the  surfaces  are  not  perfect  emitters  and  not  perfect  absorbers  and  they  are 
not  in  full  view  of  one  another,  the  departure  from  perfection  is  via  the  emissivity 
factor  F<  and  the  shape  or  arrangement  factor  F^. 

=  ^  Pa  P.  Six)  [  r/  -  ]  (2  .10) 

The  problem  of  a  finned  surface  radiating  to  outer  space  simplifies  the  foregoing 
equation  (2.10).  The  temperature  of  the  fin  becomes  T(x)  and,  because  outer  space  is 
at  about  2  K,  Tj  can  be  assumed  to  be  equal  to  0.  [Compare  Tfx)"*  =  200"  to  T/  = 
2"].  Moreover,  F^  is  equal  to  1,  the  value  of  a  small  body  in  a  large  enclosure  (small 
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Figure  2.1  Energy  Balance 


fin  in  a  large  enclosure,  outer  space).  Finally  the  emissivity  factor  for  this  case 
(small  body  in  a  large  enclosure)  is  F,  =  e(x,T).  The  emissivity  e  is  a  property  of  the 
radiating  fin  surface  and  temperature.  Thus  equation  (2.10)  is  rewritten  as 

=  a  e(x,7)  S(x)  T{xf  (2.11) 

The  difference  between  the  rates  of  heat  entering  and  leaving  the  element  in 
Figure  2, 1  by  conduction  (dq^^.b  -  dq^,,  J  must,  in  the  steady  state,  be  equal  to  the 
rates  of  heat  dissipated  from  the  face,  L  dx  by  convection  (dq„i  and  dq5v2)  and 
radiation  (dq,,  and  dqr2)  plus  any  rates  of  heat  that  arrive  from  external  sources  (dq^ 


and  dqa). 

The  presence  of  two  terms  for  dq„,  dq,,  and  dq,  permit  unequal  values  on  the 
two  faces.  Thus  with 
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where  E  consists  of  external  heat  inputs  from  solar,  Eg,  and/or  terrestrial,  Ep,  sources 
then  equation  (2.15),  the  steady  state  thermal  energy  balance  for  a  differential  element 
dx,  is  rewritten  as 


d_ 

dx 


m  5(jc)  L  ^ 
dx 


+  h^,ix,T)  ]  L  Tdx 


+  o  e(x.7)  LT*dx  +  2aLEdx  (2.17) 


B.  LIMITING  ASSUMPTIONS 

Further  assumptions  must  be  applied  to  the  generalized  differential  equation 
(2.17)  before  a  mathematical  solution  can  be  attempted.  The  assumptions  employed 
by  Murray'  and  Gardner^  are  used  for  a  starting  point  for  further  analysis.  The 
Murray-Gardner  Assumptions  listed  in  Table  1  are  essential  in  simplifying  the 
formulation  of  the  generalized  differential  equation  (2.17)  so  that  mathematical 
analysis  is  possible.  Consider  an  arbitrary  symmetrical  longitudinal  fin  shown  in 
Figure  2.2  of  length  (b-a),  depth  L,  and  width  6(x)  dissipating  radiant  heat  to  a 
constant  temperature  T,  surrounding  environment. 

The  profile  area  Ap  is  defined  as 

Ap  =  j  5(x)  dx  (2.18) 

x-h~a 
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The  initial  conditions  at  the  fin  tip  are 


T^  =  nx=a)  (2.19) 

and 

Applying  the  tenth  Murray-Gardner  assumption  (Table  1),  that  no  heat  leaves  from 
the  tip  of  the  fin,  equation  (2.20)  leads  to 

dTl  ^  Q  (2.21) 

dx 

Two  more  boundary  conditions  are  taken  at  the  fin  base 

nx=b)  =  T^,  (2.22) 

and 

q(  x=b  )  =  (2.23) 

Applying  the  fifth  and  ninth  Murray-Gardner  assumptions  (Table  1),  allows  one  to 
ignore  the  radiant  heat  flow  from  the  edges.  The  surface  area  S(x)  perpendicul  r  to 
the  radiant  heat  flow  is  given  by 

S(x)  =  2  L  Y  <2.24) 

The  face  area  S(x)  through  which  the  heat  is  dissipated  depends  on  the  arc  length  7,  a 
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function  of  the  fin  length  b. 


The  cross  sectional  area  A(x)  perpendicular  to  the  heat  flow  by  conduction  is 

A(x)  =  L  6ix)  (2-25) 


Because  the  second  Murray-Gardner  assumption  (Table  1)  states  that  the  material  is 
homogeneous,  both  the  emissivity  and  thermal  conductivity  are  assumed  to  be 
constant  and  uniform 
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TABLE  1 


THE  MURRAY-GARDNER  ASSUMPTIONS^ 


1.  The  heat  flow  into  the  fm  and  the  temperature  at  any  point  on  the  fm 
remains  constant  with  time. 

2.  The  fm  material  is  homogeneous,  its  thermal  conductivity  is  the  same  in 
all  directions  and  remains  constant. 

3.  The  heat  transfer  coefficient  between  the  fm  and  the  surrounding  medium 
is  uniform  and  constant  over  the  entire  surface  of  the  fm. 

4.  The  temperature  of  the  medium  surrounding  the  fin  is  uniform. 

5.  The  fm  width  is  so  small  compared  with  its  height  that  temperature 
gradients  across  the  fm  width  may  be  neglected. 

6.  The  temperature  at  the  base  of  the  fin  is  uniform. 

7.  There  are  no  heat  sources  within  the  fm  itself 

8.  Heat  transfer  to  or  from  the  fm  is  proportional  to  the  temperature  excess 
between  the  fm  and  the  surrounding  medium. 

9.  There  is  no  contact  resistance  between  fins  in  the  configuration  or  between 
the  fm  at  the  base  of  the  configuration  and  the  prime  surface. 

10.  The  heat  transferred  through  the  outermost  edge  of  the  fm  (the  fm  tip) 
is  negligible  compared  to  that  through  the  lateral  surfaces  (faces)  of  the  fm. 
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e(x,T)  =  e 


(2.26) 


k(T)  =  k  (2-27) 

The  contribution  from  convective  dissipation  from  the  faces  of  the  fin  is  ignored 
because  the  arbitrary  symmetrical  longitudinal  fin  is  operating  in  a  space  environment. 

C.  GENERALIZED  DIFFERENTIAL  EQUATION  FOR  A  FIN 

Combining  the  generalized  differential  equation  (2.17)  and  equations  (2.24) 
through  (2.27),  and  assuming  that  there  are  no  external  heat  inputs  (q,,  =  q,2  =  0), 
the  steady  state  thermal  energy  balance  for  the  differential  element  shown  in  Figure 
2.2  is 

AU  I  6(x)^^)  -  2  L  a  e  7(x)'‘  ^  =  0  (2.28) 

dx\  dx  )  dx 

The  temperature  excess  6{\)  is  defined  as  the  temperature  difference  between  the 
temperature  at  point  x  and  the  constant  temperature  of  the  surrounding 
environment  T,. 

e(x)  =  Tlx)  -  r,  (2.29) 

and  from  this,  it  is  seen  that 

dllx)  =  dd(x)  (2-30) 

Applying  the  temperature  excess  equations  (2.29)  and  (2.30),  the  equation  for  the 


24 


steady  state  thermal  energy  balance  for  the  differential  element  is  rewritten  as 


^  I  -  ^  ew  ^  =  0 

dx  \  k  ax 


(2.31) 


where  h,  is  defined  as  the  radiative  heat  transfer  coefficient.  Because  T,  »  0, 


=  a  e  6(jc)^ 


(2.32) 


A  fundamental  concept  of  calculus  is  that  the  length  of  arc  d7  in  a  cartesian 
coordinate  framework  is 


dy  =  [  dx^  +  dy^] 


(2.33) 


The  derivative  of  the  arc  length  is 


^  =  itoM  =  1  ♦ 

dx  Ax  \  (ox 


(2.34) 


The  +  or  -  sign  is  to  be  taken  according  as  the  arc  length  y  increases  or  decreases  as 


X  increases. 


Combining  the  equation  for  the  steady  state  thermal  energy  balance  differential 
element  (2.31)  and  the  arc  length  equation  (2.34),  the  generalized  differential  equation 
for  a  longitudinal  fin  of  arbitrary  symmetrical  profile  is 


^[6(;c)  Q(x)  [l  +  6(x)2]  =  0 

dx  dx  k 


(2.35) 


Equation  (2.35)  is  a  second  order,  linear,  homogeneous  differential  equation  with  the 
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boundary  conditions 


and 


dQj  x=a  )  ^  Q 
dx 


(2.36) 


e(  x-b )  =  e. 


(2.37) 


D.  LONGITUDINAL  FIN  OF  RECTANGULAR  PROHLE 

The  longitudinal  fin  of  rectangular  profile  is  illustrated  in  Figure  2.3.  It  should 
be  noted  that  the  origin  of  the  x  axis  is  located  at  the  tip  of  the  fin  and  not  at  the  base 
of  the  fin.  Because  the  width  is  constant  over  the  entire  profile,  the  width  function 
5(x)  is  simply 

^  (2.38) 

The  profile  area  Ap  is  written  as 


x-b 

=  [  *(*)  dx  =  b^b 
*-o 


(2.39) 


Inserting  the  width  function  (2.38)  into  the  generalized  differential  equation  for  a  fin 
equation  (2.35)  yields 


'  dx^ 


UIl 

k 


e(x)  =  0 


(2.40) 
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or 


where 


d(x)  =  0  (2.41) 

dx^ 


m 


(2.42) 


The  radiative  heat  transfer  coefficient  h,  is  assumed  to  be  a  constant  for  a  very  small 
subfin  where  Tb  *  T„. 

The  general  solution  of  the  generalized  differential  equation  (2.41)  is 


e(;c)  =  Cj  e’"*  +  Cj  e""*  (2.43) 

where  C,  and  C2  are  arbitrary  integration  constants  that  are  evaluated  by  applying  the 
initial  condition  of  equation  (2.22). 


6(6)  =  0^  =  C,  ^  +  q  c'"* (2.44) 

The  heat  flow,  qt,  into  the  base  of  the  fin  is  found  by  inserting  the  general 
solution  of  the  generalized  differential  equation  (2.43)  and  the  cross  sectional  area 
perpendicular  to  the  heat  flow  by  conduction  (2.25)  into  Fourier’s  Law  of 
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Figure  2.3  Longitudinal  Fin  of  Rectangular  Profile 

Conduction,  equation  (2.4). 

qix)  =  L  it  ^  ^[Cj  +  Cj  c"'"  ']j  (2-45) 

Simplifying  and  evaluating  at  x  =  b  with  the  initial  condition  of  equation  (2.23), 
equation  (2.45)  becomes 

=  L  6^  A:  m  (  C,  c"  *  -  Cj  c'"  ^  )  (2.46) 
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Solving  for  Cj  and  Cj  yields 


Cl 


-m  b 


^  - 


,-m  b 


2  y 


r  ^b 


(2.47) 


and 


where 


Ca 


.mb  mb 

—  ^b  -  —  Qb 
2  *  2  ^ 
b 


(2.48) 


=  L  [2  h^k  b^)  (2.49) 

is  defined  as  the  characteristic  thermal  admittance  Yi,  of  the  fm.  Likewise  the 
characteristic  thermal  impedance,  Zb,  of  the  fin  is  the  reciprocal  of  Yb 

4  =  ^  =  I  (  2  (2.50) 

An  examination  of  equations  (2.43)  through  (2.48)  reveals  that  only  two  of  the 
four  variables  are  independent  and  that  the  specification  of  any  two  variables  allows 
for  the  determination  of  the  remaining  two.  Combining  equations  (2.43)  through 
(2.48)  yields  the  two  port  inverse  thermal  transmission  matrix  for  a  longitudinal  fin  of 
rectangular  profile 


0. 

cosh  mb 

-Z^  sinh  mb 

-y^  sinh  mb 

cosh  mb 

(2.51) 
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where  the  inverse  thermal  transmission  matrix  is  defined  as 


T'  = 


cosh  mb 
-Yf,  sinh  mb 


-Zi^  sinh  mb 
cosh  mb 


(2.52) 


The  thermal  transmission  matrix  is  defined  as  the  inverse  of  the  inverse  thermal 


transmission  matrix. 


T  =  T‘> 


A  B 
C  D 


(2.53) 


where  A,  B,  C,  and  D  are  the  elements  of  the  thermal  transmission  matrix.  The 


thermal  transmission  matrix  representation  is 


0. 

cosh  mb 

sinh  mb 

0. 

sinh  mb 

cosh  mb 

(2.54) 


E.  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROnLE 

The  longitudinal  fm  of  trapezoidal  profile  is  shown  in  Figure  2.4.  Note  that  the 
origin  of  the  axis  is  not  at  the  tip  of  fm  but  located  to  the  right  of  the  tip 
where  x  =  0. 

The  width  function  5(x)  for  this  profile  is 

6(x)  =  2  X  tan  4>  (2.55) 
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where  4>  is  defined  as 


4>  =  arctan 


2  (b  -  a) 


(2.56) 


The  profile  area  Ap  is 


x-b 

=  /  aw  (6^  *  a,) 


(2.57) 


Inserting  the  width  function  (2.55)  into  the  generalized  differential  equation  for  a  fin 
[Equation  (2.35)]  yields 


or 


where 


2  X  tan(<{)) 


dx^ 


+  2  tan(<{)) 


dQjx) 

dx 


d(x)  =  0  (2.58) 

k  cos(<t>) 


0(x)  =  0 


(2.59) 


m 


k  sin(4>) 


i/z 


(2.60) 


Equation  (2.59)  represents  a  second  order,  linear,  homogenous  differential  equation. 
The  general  solution  is  expressed  in  terms  of  modified  Bessel  functions.  The  general 


31 


Figure  2.4  Longitudinal  Fin  of  Trapezoidal  Profile'^ 
solution  to  equation  (2.59)  is 


GW  =  Cj  /j  2  w  jc  2  m  X  )  (2.61) 

The  following  transformation  is  made 

u  =  2  m  (2.62) 

such  that  at  X  =  b 

=  (2.63) 
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and  at  X  =  a 


u=2ma^f^  (2.64) 

a 

After  applying  the  boundary  conditions  of  equations  (2.36)  and  (2.37)  and 
making  the  transformations  of  equations  (2.63)  and  (2.64),  the  general  solutions 
become 

9.  •  C,  W  »  C,  A:„(i,j)  (2.65) 

and 

1  k  b.  L 

9,  =  - - - [C,  7.(«^  -  C,  Ar,(a^]  (2.66) 

Simultaneous  solution  of  the  general  solution  equations  (2.65)  and  (2.66)  yield 

C,  =  *:,(»,)  0.  + 1  J  w  (2.67) 

2  k  b.Lmr 
\  »  / 

and 

q  =  7,(“*)  9,  -  A  W  9,  (2.68) 

2  k  b.L  rrr 
\  »  / 

Substituting  and  rearranging  the  general  solution  equations  (2.65)  and  (2.66)  and  the 
equations  for  C,  and  Cj  (2.67)  and  (2.68),  the  elements  of  the  thermal  transmission 
matrix  for  the  longitudinal  fin  of  trapezoidal  profile  are  written  as: 
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A  =u^  I  KJiu;)  + 


(2.69) 


B  = 


u. 


Ikh, 


— -  [W w  -  W  W] 

L  tn^ 


C  =  2  k  6L 


[  JC,(uJ  -  /,(£/„)  JC^Cupi 


D  =  «,  [  VuJ  +  I,(u^  K,(uJ] 


(2.70) 


(2.71) 

(2.72) 


F.  LONGITUDINAL  FIN  OF  TRIANGULAR  PROFILE 

The  longitudinal  fin  of  triangular  profile  shown  in  Figure  2.5  is  a  special  case  of 
the  longitudinal  fin  of  trapezoidal  profile.  There  are  many  similarities  between  the 
longitudinal  fin  of  triangular  and  trapezoidal  profiles.  The  differential 
equation  for  the  longitudinal  fin  of  triangular  profile  is  the  same  as  the  differential 
equation  for  the  longitudinal  fin  of  trapezoidal  profile  equation  (2.59). 


^  d^d(x)  ^  dd(x) 
dx^  dx 


-  e(x)  =  0 


(2.73) 
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where 


and 


(2.74) 


<t> 


arctan 


(2.75) 
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The  profile  area  Ap  is 


x=b 

.  /  SW  <fe  =  I  8, 
1*0  ^ 


(2.76) 


However  the  evaluation  of  the  arbitrary  constants  from  the  boundary  conditions 
differ  between  the  triangular  and  the  trapezoidal  profiles.  As  the  value  of  x 
approaches  zero  the  temperature  excess  at  the  tip  must  be  equal  to  some  finite 
value,  therefore 

/„(  w=0  )  =  1  (2.77) 

In  order  to  physically  obtain  a  finite  temperature  at  x,  =  u,  =  0,  Cj  in  equation 
(2.61)  must  be  equal  to  zero  because  K2(0)  is  unbounded.  The  unbounded  condition 
allows  cancellation  of  the  Ko(u)  terms  from  the  general  solutions  for  temperature 
excess  and  heat  flow. 

After  cancellation  and  rearranging  the  general  solutions  for  the  longitudinal  fin 
of  trapezoidal  profile  equations  (2.65)  and  (2.66),  the  general  solution  for  the 
longitudinal  fin  of  triangular  profile  becomes 


<ib  = 


2  k  L  I^iu^ 


(2.78) 


As  a  result  of  the  zero  tip  width,  the  edge  parameters  do  not  depend  on  the  base 
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parameters.  The  thermal  input  admittance  element  for  the  longitudinal  fin  of 
triangular  profile  is 


Y  =  —  - 


2  k  L  m} 


(2.79) 
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m.  CASCADE  ANALYSIS  ALGORITHM 


A.  GENERALIZED  TWO  PORT  TRANSMISSION  LINE  MODEL 

The  radiating  fin  is  a  distributed  heat  flow  configuration  which  is  analogous  to 
the  transmission  line  which  is  a  distributed  electrical  network.  The  general  two  port 
transmission  line  model  is  shown  in  Figure  3.1.  The  model  is  described  by  two  input 
terminals.  At  point  "a"  there  is  an  applied  voltage  V,  and  a  current  I,  flows  into  the 


Transmission 

Lina 

V. 

Vb 

! 

Figure  3.1  Two  Port  Transmission  Line  Model 
model.  There  is  also  a  pair  of  terminals  at  point  "b"  where  a  voltage  Vb  occurs  and  a 
current  Ib  flows  out  of  the  terminals.  The  two  port  model  of  the  transmission  line  can 
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consist  of  any  arrangement  of  passive  elements  within  the  two  port  but  with  no  active 
current  or  voltage  sources  within  the  two  port. 

It  can  be  shown  from  network  theory  that  the  two  port  transmission  line  model 
in  Figure  3.1  can  be  represented  by  the  following  equations: 

V^=AV^^BI^  (3.1) 


=  C  *  D 

where  the  determinant  of  the  ABCD  matrix  defined  by 


A  B 

Vt> 

C  D 

is  equal  to  unity. 


(3.2) 


(3.3) 


A  D  -  B  C  =  1  (3.4) 

The  values  A,  B,  C,  and  D  are  constants  and  are  defined  as  the  elements  of  the 
transmission  matrix. 

The  electrical  transmission  line  analogy  can  be  directly  applied  to  an  equivalent 
thermal  transmission  line  model  where  the  heat  flow  is  analogous  to  current  and  the 
temperature  is  analogous  to  voltage.  Equations  (3.1)  and  (3.2)  are  rewritten  in  the 
equivalent  thermal  transmission  line  model  as: 


=  A  B 


(3.5) 
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C  T^*  D 


(3.6) 


B.  CASCADE  ALGORITHM 

The  thermal  transmission  parameters  are  very  useful  in  describing  two  port 
networks  connected  in  a  cascaded  arrangement.  It  is  an  exercise  in  elementary  two 
port  theory.  Examine  the  two  equations  for  the  two  networks  shown  in  Figure  3.2. 


and 


9.; 

C,  D, 

'^2 

.^2 

^2. 

(3.7) 


(3.8) 


Continuity  for  the  configuration  in  Figure  3.2  dictates  that 


06,  =  K 


and 


(3.9) 


^6,  = 


(3.10) 


which  can  be  represented  as 


0^. 


e. 


(3.11) 


The  overall  thermal  transmission  parameter  matrix  for  the  two  networks 


40 


cascaded  in  Figure  3.2  is  by  definition 


A  B 

A,  s; 

Aj  B2 

C  D 

^2  ^2, 

Figure  3.2  Cascaded  Network 

and  the  thermal  transmission  matrix  for  a  cascaded  network  of  n  networks  is 

A  B 
C  D 

Equation  (3.13)  is  an  expression  of  what  is  referred  to  as  the  cascade  algorithm. 

C.  COMPUTER  ANALYSIS 

The  cascade  algorithm  expression  (3.13)  can  be  applied  to  a  fin  divided  into  n 
subfins  each  having  the  same  length.  Converting  the  temperature  excess  notation 


A,  B, 

Cl 
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back  to  the  temperature  notation  (T,  *  0),  the  cascade  algorithm  expression  (3.13)  is 
written  as 


A,  B, 

B2 

Cl  Di 

C2  D2 

The  base  temperature  and  the  base  heat  flow  qt  are  the  conditions  at  the  base  and 
the  tip  temperature  T,  and  tip  heat  flow  q,  are  the  conditions  at  the  edge. 

A  problem  arises  because  the  tip  temperature  T,  and  the  base  heat  flow  qt  are 
not  known  or  not  specified.  The  cascade  algorithm  expression  (3.14)  can  not  be 
solved  directly  because  it  is  essentially  a  linear  transformation  that  maps  tip  conditions 
to  base  conditions.  However,  a  solution  can  be  calculated  indirectly.  If  one  makes 
an  initial  estimation  for  the  tip  temperature  T,,  it  can  be  used  to  calculate  the  radiant 
heat  flow  q,  out  of  the  tip  of  the  fin. 

g,  =  o  e  63  L  r/  (3.15) 

Thus  T,  and  q,  are  the  starting  values  for  the  computations  using  the  cascade 
algorithm  expression  (3.14).  At  the  conclusion  of  the  process,  a  comparison 
involving  all  n  subfins  is  then  made  between  the  known  or  specified  base  temperature 
Tb  and  the  calculated  base  temperature.  If  the  comparison  indicates  an  error,  the  tip 
temperature  T,  is  changed  and  the  iterative  process  is  repeated  and  continues  until  the 
error  between  the  presumed  known  base  temperature  Tb  and  the  calculated  base 
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temperature  are  within  an  acceptable  range.  Figure  3.4  provides  a  complete  flow 
chart  of  the  cascade  algorithm. 

Further  amplification  concerning  the  radiative  heat  transfer  coefficient  needs  to 
be  made.  Recalling  equation  (2.32)  where 

hj.  =  o  e.  T(x)  ^  (3.16) 

it  is  not  clear  as  to  what  temperature,  T(x),  should  be  applied  to  the  computation  of  h, 
and  hence  to  the  cascade  algorithm  expression  (3.14).  The  left  and  right  faces  of  the 
differential  element  in  Figure  3.3  have  an  average  temperature. 

T  + 

T  =  ^  (3.17) 

-‘av  2 


Figure  3.3  Differential  Element 
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and  as  long  as  and  Tt  are  not  vastly  different,  equation  (3.16)  can  be  written  using 


T... 


=  o  e  Tj  (3.18) 

The  environmental  parameter  V  is  introduced  for  the  case  in  which  there  are 
portions  of  the  surrounding  environment  at  some  temperature  other  than  T,  *  0.  An 
example  of  this  is  in  the  interception  of  thermal  radiation  from  a  nearby  solar  panel  at 
some  temperature  not  near  absolute  zero.  The  environmental  parameter  is  defined  as 


where  the  parameter  Kj  allows  for  the  application  of  external  heat  inputs  from  solar, 
Ej;,  or  terrestrial,  Ep,  sources  or  both. 

=  E  a  (3.20) 

The  environmental  parameter  is  equal  to  zero  for  the  free  space  case  and  to  some 
number  greater  than  zero  but  less  than  one  for  the  non  free  space  case.  Substituting 
equation  (3.19)  into  (3.18),  the  radiant  heat  transfer  coefficient  becomes 


hr  = 


O  € 


■‘av  2  T 


(3.21) 


External  heat  inputs  (qj,  and  qj  from  solar  input  Eg  and  terrestrial  input  Ej  are  now 
applied  to  the  cascade  algorithm  expression  (3. 14)  through  the  application  of  the 
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environmental  parameter  V. 


D.  PROGRAM  DESCRIPTION 

The  computer  code  consists  of  three  analysis  programs  for  each  of  the  profiles 
considered,  and  a  MAIN  Program  that  is  used  for  general  housekeeping.  The 
PERFORMANCE  Analysis  Program  evaluates  an  entered  set  of  parameters  returning 
the  characteristics  (volume,  profile  area  Ap,  profile  number  PN,  base  heat  flow  q^, 
edge  temperature  T„  and  the  fin  efficiency  tj)  of  the  fin  along  with  the  thermal 
transmission  parameter  matrix.  The  PROFILE  AREA  Analysis  Program  also 
evaluates  an  entered  set  of  parameters  but  it  returns  a  data  set  of  100  fins  having  a 
constant  profile  area  Ap  with  the  base  fin  width  S  and  length  b  which  are  used  as  the 
input  variables.  The  OPTIMIZE  Analysis  Program  is  similar  to  the 
PERFORMANCE  Analysis  Program  with  the  exception  that  it  returns  the  optimum 
fin  geometry  which  results  in  the  minimum  mass  and  maximum  base  heat  flow  fin. 

The  MAIN  Program  provides  menu  driven  access  to  the  various  analysis 
programs.  The  user  has  the  ability  to  set  the  computational  main  program  error 
coefficient.  The  user  also  has  the  ability  to  set  the  save  data  switch  which  controls 
data  flow  into  an  ASCII  formatted  disk  file.  Commas  are  used  for  the  data  delimiter 
which  allows  file  importation  into  MATLAB,  Word  Perfect,  Draw  Perfect,  and  Math 
CAD  for  additional  analysis  or  for  plotting.  An  additional  feature  provided  is  the  unit 
conversion  module.  Table  2  list  the  available  units  that  are  provided  to  the  user  by 
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the  MAIN  Program.  Figure  3.4  provides  a  flow  chart  of  the  cascade  algorithm 
expression  (3.14)  which  is  common  to  all  of  the  analysis  programs. 

The  PERFORMANCE  Analysis  Program  requires  the  following  input  conditions 
to  be  entered: 

•  base  temperature 

•  density 

•  emissivity 

•  environmental  parameter 

•  fin  base  width 

•  fin  edge  width' 

•  fin  depth 

•  fin  length 

•  thermal  conductivity 

The  program  calculates  and  returns  the  following  values; 

•  base  heat  flow 

•  edge  temperature 

•  fin  efficiency 

•  profile  area 

•  profile  number 

•  thermal  transmission  parameter  matrix 

•  total  mass 

The  fin  efficiency  rj  and  the  profile  number  PN  are  extensively  refeired  to  in 
literature  and  are  provided  for  use  in  making  comparisons.  The  fin  efficiency  and  the 
profile  number  are  defined  as 


^1 


2  a  e  b  L  Tt,* 


(3.22) 


'Applies  for  only  the  longitudinal  fin  trapezoidal  profile. 
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Figure  3.4  Cascade  Algorithm  Flow  Chart 
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TABLE  2 


UNITS  AVAILABLE  IN  THE  MAIN  PROGRAM 


Admittance 

•  watts  /  °C,  BTU  /  hr  °F 

Area 

•  m^,  cm^,  fF,  in^,  mm^ 

Heat 

•  watts,  kw,  BTU  /  sec,  BTU  /  hr 

Impedance 

•  "C  /  watts,  °F  hr  /  BTU 

Length  and  Depth 

•  m,  cm,  ft,  in,  mm 

Mass 

•  kg,  grams,  Ibm,  oz 

Temperature 

•  K,  °C,  °F,  °R 

Thermal  Conductivity 

•  watts  /  m  K,  BTU  /  ft  hr  °F 

Width 

•  m,  cm,  ft,  in,  mm 

O  fi  J-)2  'T-  3 

PN  =  - —  (trapezoidal)  (3.23) 

k  6. 


(rectangular)  (3.24) 
the  save  switch  is  switched  on,  the  following  differential  element  conditions  are 


48 


written  to  a  disk  file: 


•  base  heat  flow 

•  base  temperature 

•  distance  from  the  tip  of  the  fin 

The  OPTIMIZE  Analysis  Program  requires  the  following  input  conditions  to  be 
entered: 


•  base  temperature 

•  density 

•  emissivity 

•  environmental  parameter 

•  fin  edge  width^ 

•  fin  width 

•  thermal  conductivity 

The  program  calculates  and  returns  the  following  optimum  values: 

•  base  heat  flow 

•  edge  temperature 

•  fin  base  width 

•  fin  efficiency 

•  fin  length 

•  profile  area 

•  profile  number 

•  thermal  transmission  parameter  matrix 

•  total  mass 

The  program  saves  the  same  quantities  that  the  PERFORMANCE  Analysis  Program 
saves  depending  on  the  condition  of  the  save  switch. 

The  PROFILE  AREA  Analysis  Program  requires  the  following  input  conditions 


^Applies  for  only  the  longitudinal  fin  trapezoidal  profile. 
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to  be  entered; 


•  base  temperature 

•  density 

•  emissivity 

•  environmental  par?  meter 

•  profile  area 

•  starting  fin  base  width  value 

•  ending  fin  base  width  value 

•  taper  ratio^ 

•  fin  depth 

•  thermal  conductivity 

The  program  calculates  and  writes  to  a  disk  file  the  following  parameters  for  100  fins 
having  the  same  profile  area  Ap  within  the  specified  width  range: 

•  base  heat  flow 

•  base  width 

•  base  temperature 

•  edge  temperature 

•  fin  efficiency 

•  thermal  transmission  parameter  matrix 

•  total  mass 

The  analysis  programs  all  use  the  thermal  transmission  matrices  computed  in 
Chapter  2  in  the  cascade  algorithm  with  the  exception  of  the  longitudinal  fin  of 
triangular  profile.  The  triangular  and  trapezoidal  profile  analysis  are  very  similar. 

The  analysis  program  uses  the  thermal  admittance  equation  (2.79)  to  compute  the  base 
conditions  for  the  first  subfin  having  a  thickness  of  5,.  The  subfin  is  located  at  the  tip 
of  the  fin.  The  base  conditions  are  then  applied  and  used  as  the  tip  conditions  in  the 
thermal  transmission  matrix  for  the  trapezoidal  profile  (2.69)  through  (2.72).  The 


^Applies  for  only  longitudinal  fin  of  trapezoidal  profile. 
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cascade  algorithm  is  then  used  to  calculate  the  thermal  transmission  matrices  for  the 
remaining  99  subfins  in  order  to  finally  determine  the  overall  thermal  transmission 
matrix  and  the  base  conditions.  The  trapezoid  that  is  being  used  in  the  calculation  has 
an  almost  zero  tip  width  compared  to  the  base  width  such  that 

=  100  6^  (3.25) 


Borland  Turbo  C++  was  used  to  write  the  computer  code.  The  computer  code 
listing  appears  in  Appendix  A.  Full  compiler  optimization  was  selected  in  order  to 
decrease  the  run  times.  Average  run  times  on  an  IBM  486/487  machine  appear  in 
Table  3. 


TABLE  3 

AVERAGE  PROGRAM  RUN  TIMES 


PERFORMANCE 

(sec) 

PROFILE  AREA 
(sec) 

OPTIMIZE 

(sec) 

Rectangular 

4 

97 

4 

Trapezoidal 

4 

600 

4 

Triangular 

5 

625 

5 
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IV.  PROGRAM  ANALYSIS  RESULTS  AND  DISCUSSION 


A.  PROGRAM  ERROR  TOLERANCE 

In  order  to  allow  for  meaningful  interpretation,  an  understanding  of  the  accuracy 
of  the  results  produced  from  the  cascade  algorithm  is  vital.  A  test  to  show 
convergence  is  used  to  illustrate  the  approximate  relative  error. 

Table  4  provides  an  illustration  of  how  the  error  coefficient  in  the  MAIN 
Program  affects  the  accuracy  of  the  program  results.  Table  4  is  produced  by 

TABLE  4 


PROGRAM  ERROR  TOLERANCE 


Program 

Error 

Ta 

«  % 

Qb 

*  % 

Coeff  % 

(K) 

error 

(watts) 

error 

1.000 

282.18 

0.0106 

926.28 

2.0918 

0.500 

282.24 

0.0106 

931.83 

1.5052 

0.100 

282.19 

0.0071 

946.36 

0.0307 

0.050 

282.22 

0.0035 

946.94 

0.0920 

0.010 

282.22 

0.0035 

946.10 

0.0032 

0.005 

282.21 

0.0000 

946.05 

0.0021 

0.001 

282.21 

0.0000 

946.07 

0.0000 

running  the  same  input  parameters  many  times  while  decreasing  the  error  coefficient 
in  the  MAIN  Program.  The  PERFORMANCE  Analysis  Program  for  the  longitudinal 
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fin  of  rectangular  profile  was  used.  The  Table  4  results  are  applicable  for  any  of  the 
analysis  programs  because  of  the  commonality  of  the  cascade  algorithm.  As  the  error 
coefficient  in  the  MAIN  Program  decreases,  the  output  value  converges.  The 
approximate  relative  error  e,  is  defined  as 

True  Value  -  Approximate  Value 
True  Value 

The  exact  true  value  is  unknown,  however  the  value  produced  when  the  error 
coefficient  in  the  MAIN  Program  is  equal  to  0.001  can  be  assumed  to  be 
approximately  equal  to  the  exact  true  value. 

Examining  the  results  illustrated  in  Table  4,  one  can  conclude  that  when  the 
error  coefficient  in  the  MAIN  Program  is  less  than  or  equal  to  0. 1  the  approximate 
relative  error  observed  in  the  output  is  generally  much  less  than  the  error  coefficient 
in  the  main  program. 

B.  GRAPHICAL  EXAMPLES 

The  purpose  of  this  section  is  to  describe  and  provide  examples  on  how  each  of 
the  three  analysis  programs  can  be  used.  The  computer  code  listings  are  provided  in 
Appendix  A. 

1.  Example  One:  Longitudinal  Fin  of  Rectangular  Profile  Analysis 

A  longitudinal  fin  of  rectangular  profile  is  0.5  meter  long,  1.0  centimeter 
thick,  and  1.0  meter  deep.  The  fin  is  made  of  aluminum  (k  =  209.4  watt  /  m  K, 


100  (4.1) 
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p  =  2.6  kg  /  dm^,  and  e  =  0.85).  If  the  base  temperature  is  400  K,  what  are  the  tip, 
midpoint,  and  base  conditions  (error  <  0.1%)  both  in  the  free  space  and  non  free 
space  environments? 

Because  the  geometry  of  the  fin  is  already  known,  the  PERFORMANCE 
Analysis  Program  is  used  for  the  evaluation.  The  save  switch  in  the  MAIN  Program 
must  be  placed  in  the  save  position  in  order  to  later  retrieve  the  midpoint  conditions 
from  the  disk  file.  The  computer  output  display  for  the  free  space  environment 
(V=0)  appears  in  Table  5. 

Figures  4.1  and  4.2  provide  a  detailed  graphical  representation  of  the 
required  conditions  both  in  the  free  space  and  non  free  space  environments.  The  base 
conditions  are  read  from  where  the  length  is  equal  to  50  centimeters  and  the  tip 
conditions  are  read  from  where  the  length  is  equal  to  0  centimeter.  In  the  free  space 
environment  the  tip,  midpoint,  and  base  conditions  for  this  particular  fin  geometry  are 
listed  in  Table  6. 
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TABLE  5 


PERFORMANCE  ANALYSIS  OUTPUT  DISPLAY 


THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE 


Density 
Fin  Depth 
Total  Mass 
Fin  Length 
Fin  Width 
Emissivity 
Profile  Area 
Profile  Number 
Base  Heat  Flow 
Fin  Efficiency 
Base  Temperature 
Edge  Temperature 
Thermal  Conductivity 
Environmental  Parameter 


2.60  kg  /  dm^3 
1.00  m 
13.00  kg 
0.50  ro 
1.00  cm 
0.85 

50.00  cm'‘2 
0.607 

732.85  watts 
0.594 
400.07  K 
324.45  K 

209.40  watts  /  m  K 
0.00 


Transmission  Parameter  Matrix 


A 

B 


1.2287 

2.2374  watt  /  deg  C 


B  =  0.2583  deg  C  /  watt 
D  =  1.2842 


TABLE  6 


FIN  CONDITIONS  IN  FREE  SPACE 


Location 

Temp 

q 

(K) 

(watts) 

Base 

400.07 

732.85 

Midpoint 

341.66 

292.33 

Tip 

324.45 

10.68 

Figure  4.1  Heat  Dissipation  From  A  Rectangular  Profile  Fin 
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Length  (cm)  Distance  From 


Figure  4.2  Temperature  Distribution  For  A  Rectangular  Profile  Fin 
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Equation  (3.4)  is  used  below  to  verify  that  the  determinant  of  thermal 
transmission  matrix  listed  in  Table  5  is  indeed  equal  to  unity. 

1.2287  0.2583 

det 

2.2374  1.2842 

2.  Example  Two:  Longitudinal  Fin  of  Rectangular  Profile  Area  Analysis 

A  longitudinal  fin  of  rectangular  profile  1.0  meter  deep  is  made  of 
aluminum  (k  =  209.4  watt  /  m  K,  p  =  2.6  kg  /  dm\  and  e  =  0.85).  It  is  required 
to  remove  a  minimum  of  500  watts  of  heat  from  the  base  of  a  fin  at  a  temperature  of 
400  K.  Determine  the  geometry  of  the  fin  using  several  values  for  the  profile  area 
A,. 

The  PROFILE  AREA  Analysis  Program  is  used  because  the  geometry  of 
the  fin  is  unknown.  Figure  4.3  provides  a  graphical  representation  which  shows  the 
constant  profile  area  Ap  relationship  as  a  function  of  width  versus  heat  flow.  Any  of 
the  dimensions  listed  in  Table  7  would  meet  the  minimum  requirement  to  remove  the 
500  watts  of  heat  from  the  base  of  the  fin. 

Examining  Figure  4.3,  one  can  see  that  the  maximum  heat  flow  occurs  at 
a  particular  width  value.  More  heat  is  dissipated  as  the  profile  area  Ap  is  increased. 
The  value  of  the  width  at  the  point  of  maximum  heat  flow  also  increases  as 
the  profile  area  Ap  increases.  The  point  of  maximum  heat  flow  is  considered  as  the 
optimum  geometry  (maximum  heat  flow,  minimum  mass)  for  the  particular  profile 
area  Ap.  A  summary  of  the  optimum  geometry  points  extrapolated  from  the  data  files 


1.0000 


(4.2) 
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(s»wi)  MOid  1WH 


FIGURE  4.3  Profile  Area  Analysis 


TABLE  7 


PROHLE  AREA  ANALYSIS  POINTS 


Ap 

b 

qb 

(cm^) 

(cm) 

(m) 

(watts) 

25 

0.334 

0.748 

507.4 

2.041 

633.7 

2.762 

771.2 

0.334 

8.982 

523.1 

produced  from  the  PROFILE  AREA  Analysis  Program  appear  in  Table  8. 

TABLE  8 

OPTIMUM  GEOMETRY  POINTS  FOR  THE  RECTANGULAR  PROHLE 


Ap 

bo 

MasSo 

Qbo 

Qbo/m 

(cm^) 

(cm) 

(m) 

(kg) 

(watts) 

(w/kg) 

25 

0.646 

0.387 

6.50 

579.76 

89.19 

100 

1.543 

0.648 

26.00 

924.53 

35.56 

200 

2.440 

0.819 

52.00 

1165.83 

22.42 

300 

3.103 

0.967 

80.00 

1334.40 

16.68 

3.  Example  Three;  Comparison  of  Various  Profiles 

A  longitudinal  fin  1.0  meter  deep  having  a  profile  area  Ap  equal  to  50  cm^ 
is  made  of  aluminum  (k  =  209.4  watt  /  m  K,  p  =  2.6  kg  /  dm’,  and 
e  =  0.85).  The  base  temperature  is  held  at  400  K.  Compare  the  relative 
performance  between  the  longitudinal  fins  of  rectangular,  trapezoidal,  and  triangular  profiles. 


The  same  approach  is  used  as  in  the  previous  example.  A  longitudinal  fin 
of  trapezoidal  profile  having  a  taper  ratio  TR  =  0.5  is  used.  The  taper  ratio  is 
defined  as 


Figure  4.4  provides  a  graphical  representation  which  shows  the  relationship  between 
the  three  profiles  being  considered. 

The  curve  in  Figure  4.4  produced  by  the  longitudinal  fin  of  trapezoidal 
profile  is  midway  between  the  curves  produced  by  the  longitudinal  fins  of  rectangular 
and  triangular  profiles.  Of  more  importance  is  the  fact  that  at  the  location  of  the 
optimum  geometry  points,  the  longitudinal  fin  of  triangular  profile  produces  the 
largest  heat  flow  while  the  longitudinal  fin  of  rectangular  profile  produces  the  smallest 
heat  flow.  The  optimum  geometry  points  extrapolated  from  the  disk  file  generated 
from  the  PROFILE  AREA  Analysis  Program  app)ear  in  Table  9.  Not  only  is  there  an 
optimum  fin  geometry  for  a  particular  profile  area  Ap  but  there  is  also  an  optimum 
profile  shape  which  produces  maximum  heat  flow. 

4.  Example  Four:  Optimum  Geometry  Comparison 

A  longitudinal  fin  2.0  meters  deep  (€  =  0.95)  is  required  to  remove  at 
least  1000  watts  from  a  spacecraft  in  a  free  space  environment  (V=0).  The  base 
temperature  is  held  at  400  K.  Using  the  list  of  materials  apj)earing  in  Table  10 
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Figure  4.4  Comparisons  of  Various  Profiles 
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TABLE  9 


OPTIMUM  GEOMETRY  POINTS  FOR  THE  VARIOUS  PROFILES 


Taper 

5b 

b 

Mass 

Qb 

Qb/m 

Profile 

Ratio 

(cm) 

(m) 

(kg) 

(watts) 

(W/kg) 

RECT 

1.0 

0.958 

0.522 

13.00 

732.47 

56.343 

TRAP 

0.5 

1.270 

0.524 

13.00 

772.93 

59.456 

TRI 

0.0 

1.309 

0.763 

13.00 

826.78 

63.598 

conduct  a  heat  flow  per  unit  mass  comparison  using  the  longitudinal  fins  of 
rectangular  and  triangular  profiles. 


TABLE  10 

MATERIALS  USED  IN  THE  OPTIMUM  COMPARISON* 


Chemical 

P 

k 

Material 

Symbol 

kg/dm^ 

w  /  m  K 

Aluminum 

A1 

2.6 

209.4 

Copper 

Cu 

8.8 

384.0 

Iron 

Fe 

7.2 

58.0 

Magnesium 

Mg 

1.7 

157.0 

Molybdenum 

Mo 

10.2 

145.0 

Silicon 

Si 

2.3 

83.0 

Silver 

Ag 

10.5 

407.0 

Zinc 

Zn 

6.8 

140.0 
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The  OPTIMIZE  Analysis  Program  is  used  for  the  analysis.  The 
parameters  were  entered  individually  into  the  OPTIMIZE  Analysis  Program  and  the 
results  are  listed  in  Table  11. 


TABLE  1 1 

OPTIMUM  COMPARISON  RESULTS 


5bo 

Mas  So 

qbo 

qbo/m 

Material 

cm 

cm'' 

kg 

watts 

w  /  kg 

Ag 

RECT 

0.40 

21.73 

45.63 

997.86 

21.86 

Ag 

TRI 

0.54 

15.09 

31.69 

995.45 

31.41 

A1 

RECT 

0.79 

42.24 

21.96 

998.96 

45.49 

A1 

TRI 

1.05 

29.33 

15.25 

995.47 

65.27 

Cu 

RECT 

0.43 

23 . 03 

40.54 

997 . 94 

24 .61 

Cu 

TRI 

0.57 

15.99 

28.15 

995.45 

35.36 

Fe 

RECT 

2.84 

152.49 

221.11 

1003.99 

4 . 54 

Fe 

TRI 

3.79 

105.89 

153.54 

995.75 

6.48 

Mg 

RECT 

1.05 

56.33 

19.60 

999.74 

51.01 

Mg 

TRI 

1.40 

39.12 

13.61 

995.49 

73 . 14 

Mo 

RECT 

1.14 

61 . 00 

124.43 

1000.00 

8.04 

Mo 

TRI 

1.51 

42.36 

86.41 

995.54 

11.52 

Si 

RECT 

1.99 

106.56 

49.66 

998.41 

20.11 

Si 

TRI 

2 . 65 

73.99 

34.48 

995.59 

28.87 

Zn 

RECT 

1.18 

63 . 17 

85.92 

1000.12 

11.64 

Zn 

TRI 

1.57 

43.87 

59.66 

995.50 

16.68 
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Figure  4.5  shows  a  bar  graph  representation  of  the  heat  flow  per  unit  mass 
relationship  using  the  various  materials  listed  in  Table  10.  It  is  evident  from  Figure 
4.5  that  the  optimized  longitudinal  fm  of  triangular  profile  made  of  magnesium  is  the 
most  efficient  material  considered  of  those  listed  in  Table  10,  with  aluminum  being  a 
close  second.  The  melting  temperature  of  both  materials  (*  900  K)  is  far  above  the 
base  temperature  (400  K)  stated  in  the  problem  statement.  Comparing  Tables  10  and 
11  materials  having  a  small  value  for  the  density  (p  <  3.0  kg  /  dm^)  and  a  large 
value  for  the  thermal  conductivity  (  k  >  150  watt  /  m  K)  produce  large  values  for 
heat  flow  per  unit  mass. 

C.  DERIVATION  OF  W  ORKING  OPTIMUM  GEOMETRY  EQUATIONS 

As  shown  in  Tables  8  and  9  optimum  geometry  points  can  be  determined  using 
the  PROFILE  AREA  Analysis  Program.  The  overall  success  of  finding  the 
dimensions  for  the  optimum  geometry  using  the  PROFILE  AREA  Analysis  Program 
has  several  limitations.  The  first  limitation  is  that  the  starting  and  ending  fin  base 
width  values  must  be  known.  And  it  can  be  argued  that  a  second  limitation  is  that 
excessive  analysis  program  run  times  which  might  exceed  ten  minutes. 

Relationships  for  the  optimum  geometry  for  the  longitudinal  fins  of  rectangular, 
trapezoidal,  and  triangular  profiles  would  remove  both  limitations.  Kern  and  Kraus’ 
and  Chung  and  Nguyen'”  stated  that  the  dimensions  for  optimum  fin  geometry  can 
be  expressed  by 
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Figure  4.5  Optimum  Profile  Comparison  Using  Various  Materials 
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(4.4) 


K  =  — — — 

L  o  e  r/ 


koeT,^ 


isi' 


(4.5) 


^Po  =  P^ 


1 

k 


(4.6) 


where  Pi(TR),  PaCTR),  and  P3(TR)  are  defined  as  polynomials  in  which  the  taper  ratio 
TR  is  the  independent  variable. 

In  general  the  polynomial  P  is  expressed  in  the  form  of 


P(JR)  =  Cj  TR  +  q  TR  "-2  +  +  C„  7/?  °  (4-7) 

where  n  is  the  degree  of  tne  polynomial. 

The  polynomials  P,,  P^,  and  Pj  are  determined  from  running  the  PROFILE 
AREA  Analysis  Program  using  a  slightly  modified  code.  The  computer  code  listing 
appears  in  Appendix  B.  A  random  number  generator  is  used  to  generate  the  input 
parameters.  The  random  number  generator  produces  input  parameters  in  the  range 
listed  in  Table  12.  Murray-Gardner  (Table  1)  assumption  number  five  (the  fin  width 
is  small  compared  with  its  depth  and  length)  is  checked  in  order  to  insure  that  the 
original  assumptions  are  kept  valid. 
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TABLE  12 


RANGE  OF  THE  RANDOM  INPUT  PARAMETERS 

0.01  <  e  <  0.99 
20  <  L  (  cm  )  <  220 
10  <  k  (  watt  /  m  K)  <  210 

100  <  Tb  (  K  )  <  600 

25  <  qb  (  watts  )  <  1025 
0.01  <  Taper  Ratio  TR  <  0.99 


The  program  enters  into  a  loop  calculating  the  heat  flow  Qb  using  the  random 
input  parameters.  The  program  remains  in  this  loop  until  the  maximum  heat  flow  qb 
is  found.  Once  found  the  program  writes  the  random  input  parameters  and  the  value 
determined  to  be  the  maximum  heat  flow  qb  to  a  disk  file.  Another  set  of  random 
input  parameters  is  generated  and  the  maximum  heat  flow  qb  once  again  is  calculated. 
This  sequence  continues  until  100  optimum  data  sets  are  written  to  a  disk  file. 

Three  data  sets  using  the  random  input  parameters  are  produced,  one  each  for  the 
longitudinal  fins  of  rectangular,  trapezoidal,  and  triangular  profiles. 

The  last  step  in  determining  the  polynomials  P,,  P2,  and  P3  is  to  fit  the  collected 
data  stored  in  the  disk  file  to  a  "best"  polynomial  fit  having  the  smallest  degree  d 
using  the  least  square  principle.  Equations  (4.4)  through  (4.6)  are  rewritten  as 
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P(TR\  = 


(4.8) 


Lb^oe  7/ 
<lb 


P(TR\  =  k  a  eT^ 


(4.9) 


P{TR), 


=  Ap^  k 


-T 


(4.10) 


Equations  (4.8)  through  (4.10)  are  placed  into  a  MATLAB  "m"  file.  Each  of  the 
collected  data  files  are  read  into  MATLAB  and  were  evaluated  using  the  MATLAB 
polynomial  fitting  function  polyfit  with  the  taper  ratio  TR  as  the  independent  variable 
using  varying  degrees  d  in  order  to  find  the  polynomials  P,,  Pj,  and  Pj  which  "best" 
approximate  the  previously  collected  data  from  the  disk  file.  The  results  appear  in 
Table  13.  The  listing  of  the  MATLAB  evaluation  m  files  appear  in  Appendix  C. 

The  standard  error  £,  is  defined  as 

1/2 

(4.11) 

where  yj  is  the  computed  value  of  y^  and  there  are  N  values. 


E.  = 


E  Cv,  -  V,) ' 

N  -  1 
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TABLE  13 


COEFTICIENTS  OF  THE  OPTIMUM  POLYNOMIALS 


Degree 

Standard 

Profile 

d 

c, 

Q 

C3 

Error 

PI  RECT 

2 

0.8675 

0.0000 

0.0000 

3.54e-3 

PI  TRAP 

2 

-0.0290 

0.8138 

0.0000 

2.12e-3 

PI  TRI 

2 

0.0000 

0.9040 

0.0000 

2.55e-3 

P2  RECT 

2 

1.8648 

0.0000 

0.0000 

1.63e-2 

P2  TRAP 

3 

0.3881 

-1.0324 

2.7522 

1.22e-2 

P2  TRI 

2 

0.0000 

2.4852 

0.0000 

1.40e-2 

P3  RECT 

2 

1.6178 

0.0000 

0.0000 

2.06e-2 

P3  TRAP 

2 

0.5147 

1.1382 

0.0000 

1.24e-2 

P3  TRI 

2 

0.0000 

1.1234 

0.0000 

9.46e-3 

The  equation  for  the  optimum  profile  area  for  either  ti.  'ectangular, 
triangular,  or  trapezoidal  profiles  is  written  as 


=  (  1.1382  +  0.5147  TR  ) 


1 


k 


L  T 


(4.12) 
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V.  CONCLUSIONS 


This  thesis  has  described  a  new  technique  that  uses  a  parameterization  called  the 
thermal  transmission  matrix  to  represent  a  single  fm  to  evaluate  the  performance  of 
several  fins.  The  algorithm  is  based  on  the  simple  ideal  of  subdividing  the  fin  into 
100  subfins  and  then  the  thermal  transmission  matrix  is  calculated  for  each  subfin. 
Then,  by  a  simple  matrix  multiplication  of  the  individual  subfin  thermal  transmission 
matrices,  the  overall  thermal  transmission  matrix  is  computed.  Not  only  will  this 
produce  the  base  and  tip  conditions  of  the  overall  fin  entity,  but  the  temperature 
distribution  and  heat  flow  at  any  point  between  the  fin  base  and  the  fin  tip  can  be 
easily  determined. 

Several  conclusions  can  be  made  from  the  Chapter  4  Program  Analysis  Results 
and  Discussion.  It  was  observed  that  the  profile  area  Ap,  increases  as  the  cube  of  the 
heat  flow  qt.  In  order  to  double  the  heat  flow  q^,  the  profile  area  Ap  must  be 
increased  eight  times. 

Material  selection  is  very  important.  The  weight  of  the  fin  is  proportional  to 
the  density  and  the  thermal  conductivity  is  inversely  proportional  to  the  profile  area. 
Materials  having  a  thermal  conductivity  greater  than  150  watt  /  m  K  and  a  density 
less  than  3  kg  /  dm^  are  generally  the  best  suited  materials.  For  the  same  material, 
environmental  conditions,  heat  flow,  and  base  temperature,  the  optimum  longitudinal 
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fin  of  triangular  profile  requires  only  69  percent  as  much  material  as  the  optimum 
longitudinal  fin  of  rectangular  profile,  with  the  optimum  longitudinal  fin  of 
trapezoidal  profile  falling  in  betw'een. 

Simple  working  relationships  to  compute  the  optimum  fin  dimensions  are  found. 
Table  16  compares  the  findings  from  this  thesis  with  the  results  from  previous  work 
by  others.  The  differences  between  the  optimum  dimension  polynomial  coefficients 
listed  in  Table  16  are  very  small,  which  verifies  the  accuracy  of  the  thermal 
transmission  matrix  parameterization  technique. 

TABLE  16 

OPTIMUM  DIMENSION  POLYNOMIAL  COEFFICIENTS 


PI 

P2 

P3 

RECTANGULAR 

.8675 

1.8648 

1.6178 

This  Thesis 

.8846 

1.8473 

1.6341 

Chung  and  Nguyen" 

.8842 

1.8579 

1.6056 

Liu" 

.8846 

1.8473 

1.6341 

Bartas  and  Sellers" 

TRIANGULAR 

.9040 

2.4852 

1.1234 

This  Thesis 

.9485 

2.3104 

1.0957 

Chung  and  Nguyen 

.9546 

2.3027 

1.0991 

Nilson  and  Curry" 

.9091 

2.4255 

1.1025 

Kern  and  Kraus" 

This  new  approach  to  fin  synthesis  and  analysis  permits  an  exact  solution  to  any 
degree  of  accuracy  that  the  user  wishes  to  specify  and  produces  the  result  in  a 
computational  efficient  manner. 
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APPENDIX  A  COMPUTER  CODE  LISTING 


/* 


Main  Program 


*/ 


^include  <process.h> 
#include  <conio.h> 
#include  <stdlib.h> 
^include  <stdio.h> 

#define  Yes  1 
#define  No  0 


char 

•SI  [] 

= 

< 

II 

watt  /  deg  C",  " 

BTU  /  hr  deg  F"  >; 

char 

*S2[] 

s 

{ 

II 

m  ■  2",  "  cm  '  2" 

"  ft  ‘  2",  "  in  ‘ 

2",  "  im 

char 

*S3[] 

= 

< 

II 

watts",  "  kw",  " 

BTU  /  sec",  "  BTU  / 

hr"  >; 

char 

•S4n 

= 

< 

II 

deg  C  /  watt".  " 

deg  F  hr  /  BTU"  ); 

char 

*S5[] 

= 

{ 

II 

m",  "  cm",  "  ft". 

"  in",  "  mm"  }; 

char 

•S6[) 

= 

i 

II 

kg",  "  grams",  " 

tbm",  "  02"  >; 

char 

•s7n 

= 

i 

II 

deg  K",  "  deg  C", 

"  deg  F",  "  deg  R" 

>; 

char 

•ssn 

= 

{ 

II 

watt  /  m  deg  k", 

"  Blu  !  ft  hr  deg  F' 

■  ); 

char 

•S9C] 

( 

M 

No",  "  Yes"); 

char 

•on 

= 

{ 

“<1>  watt  /  deg 

C  <2>  BTU  /  hr  deg 

f", 

"<1>  m  "  2  <2>  cm  '  2  <3>  ft  *  2  <4>  in  '  2  <5>  rm  ‘  2", 

"<1>  watts  <2>  kw  <3>  BTu  /  sec  <4>  BTU  /  hr", 

"<1>  deg  C  /  watt  <2>  deg  f  hr  /  BTU", 

"<1>  m  <2>  cm  <3>  ft  <4>  in  <5>  mm", 

"<1>  kg  <2>  grams  <3>  tbm  <4>  oi", 

"<1>  deg  K  <2>  deg  C  <3>  deg  f  <4>  deg  R", 

"<1>  watt  /  m  deg  K  <2>  BTu  /  ft  hr  deg  F", 

"<1>  No  <2>  Yes"  >; 

void  print_main_menu(),  print_menu_one( ),  print_menu_two( ); 

void  setupO,  overlaysO,  update_uni  ts(  ); 

float  error=0.1; 

int  T-0,  L=0,  M=0,  S=0,  H=0,  d=4,  K=0,  A=1,  2=0,  Y=0; 

char  main_key=0; 

int  mainOl 

FILE  *data_in; 

if  (  Idata  in  =  f  openC'uni  ts.  ext",  "r") )  !=  NULL){ 

fscanf(data_in,"%d  %d  %d  %d  %d  Xd  Xd  Xd  Xd  %d  %f",&T,&L,SM,SS,SH,Sd,&K,&A,&Z,&Y,Serror),- 
fclose(data_in);  > 

pr  i nt_mai n_menu( ) ; 

do( 

V  if  (main_key  ==  '!')  setupO; 

if  (main_key  >=  '2'  &&  main_key  <=  '4')  overlaysO;  ) 

while  ((  main_key  =  getchl)  )  i=  '5'); 

windowC 1 , 1 ,80,25); 

clrscrO; 

_setcursortype(_NORHALCURSOR ) ; 
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return  0;> 


void  overlays( void){ 
char  key=0; 

int  Quit  =  No,  result; 
pr  i  nt_inenu_one( )  ; 
do 


switch  (key){ 

case  'A': 

clrscrO; 

Qui  t  = 

Yes; 

if 

(niain_key 

=  =  'ZO 

result 

= 

spawnl(P_WAlT, 

••RECT.EXE",  NULL); 

if 

(main_key 

==  •30 

result 

= 

spawnl(P_UAIT, 

“TRAP.EXE",  NULL); 

if 

(main  key 

==  'A') 

result 

= 

spawnl(P_WAIT, 

••TRI.EXE",  NULL); 

break; 

case  'B': 

clrscrO; 

Quit  = 

Yes; 

if 

(main_key 

==  '2') 

result 

= 

spawnl(P_WAlT, 

••ARECT.EXE",  NULL); 

if 

(main_key 

==  •3') 

result 

= 

spawnl(P_WAIT, 

••ATRAP.EXE",  NULL); 

if 

(main  key 

==  'A') 

result 

= 

spawnl(P_WAIT, 

"ATRI .EXE",  NULL); 

break; 

case  'C  : 

clrscr(  ); 

Qui  t  = 

Yes; 

if 

(main_key 

==  >2') 

resul t 

= 

spawnl(P_WAIT, 

"0RECT.EXE",  NULL); 

if 

(niain_key 

==  •3') 

result 

spawnl(P_WAIT, 

"0TRAP.EXE",  NULL); 

if 

(main_key 

==  'A') 

resul t 

= 

spawn 1 (P_UAI T , 

"01RI.EXE",  NULL);) 

while  ((  key  =  toupper(getch( ) )  )  •=  'D'  &&  Quit  ==  No); 

print_main_menu< ); 
return;) 

void  setupCvoidX 
char  key=0; 


print_menu_two( ); 
do 

switch  (key){ 


case  'A': 


case  'B' : 


case  'C: 


case  'O': 


case  'E': 


case  'F': 


case  'O': 


gotoxy(2, 17);  clreoU);  cprintf ("/'.s",Ot0]  ); 
while  (key  1=  *1'  &&  key  !=  '2')  key  =  getchO; 

Y  =  key  -  ' 1 ' ;print_menu_two( );  break; 

gotoxyC  1 , 17);  clreoU);  cprintf ("%s'', Oil)  ); 

while  (key  !=  '1'  &&  key  !=  '2'£&  key  •=  '3'£&  key  !=  key  !=  '5') 

key  =  getch( ); 

A  =  key  -  '  1  '  ;print_nienu_two( );  break; 
gotoxy(  1 , 17);  clreolO;  cprintf ("%s", 0(2]  ); 

while  (key  '=  '1'  &&  key  '=  '2'&S  key  1=  '3'&&  key  !=  'A')  key  =  getchO; 
M  =  key  -  '  1 '  ;print_(iienu_two( );  break; 

gotoxy(2, 17);  clreoK);  cprintf  ("%s", 0(3)  ); 
while  (key  !=  '1'  &&  key  !=  '2')  key  =  getchO; 

Z  =  key  -  '  1 ' ;print_nienu_two( );  break; 

gotoxy(2, 17);  clreolO;  cprintf ("Xs", 0(4]  ); 

while  (key  !=  '1'  &&  key  !=  '2'£&  key  !=  '3'&&  key  !=  'A'&&  key  !=  '5') 
key  =  getchO; 

L  =  key  -  '  1  ‘ ;print_nienu_two();  break; 
gotoxy(3, 17);  clreolO;  cprintf ("%s", 0(5)  ); 

while  (key  !=  '1'  &&  key  !=  •2'S&  key  !=  '3'&&  key  !=  'A')  key  =  getchO; 
M  =  key  -  ' 1 ' ;print_menu_two();  break; 

gotoxy(  1 , 17);  clreolO;  cprintf ("%s", 0(6]  ); 

while  (key  !=  '1'  &&  key  !=  '2'&&  key  !=  '3'&&  key  !=  'A')  key  =  getchO; 
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T  =  key  -  ' 1 ‘ ;pri nt_menu_tHO( );  break; 

case  'H';  90toxy(1,17);  clreoU);  cprinrf("Xs‘',0[7]  ); 

while  (key  !=  '1'  &&  key  !=  '2')  key  =  getchO; 

K  =  key  -  ‘ 1 ' ;print_menu_two( );  break; 

case  'I':  gotoxy(2, 17);  clreoK);  cprintf  (“%s'‘,0[4)  ); 

while  (key  !=  '1'  &£  key  •=  •2‘&S  key  !=  ‘3‘i&  key  !=  key  !=  '5') 

key  =  getch( ); 

d  =  key  -  '  t ‘ ;print_menu_two();  break; 

case  'J*:  gotoxy(12, 17);  clreolO;  cprintf (•'%&“, 0(8]  ); 

while  (key  !=  '1‘  &&  key  !=  *2*)  key  =  getchO; 

S  =  key  -  ' 1 ' ;print_menu_two( );  break; 

case  '<■:  gotoxy(22, 12);  clreolO;  scanf ("%(", ierror); 

pr i nt_menu_two( ) ; ) 

while  ((  key  =  toupper(getch( ) )  )  •=  'L'); 

update_uni ts( ); 
pr  i nt_(na i  n_menu( )  ; 
return; > 

void  update_uni ts(void)( 

FILE  *unit; 

unit  =:  f  openC'uni  ts  .ext",  "w" )  ; 

fprintf(unit,"%d  %d  %d  %d  lid  7=d  %d  %d  >io  ‘4d  %f ",  T ,  L  ,M,  S,  H,  d,k.  A,  2,  Y,  error )  ; 
f closeluni t ) ; 
return; ) 

void  print_main_menu(void)( 

_setcursortype(_NOCURSOR ) ; 
w indow( 1,1, 80,25 ); 
clrscr(  ); 

gotoxy(27,3);  cprintf ("EXTENDED  SURFACE  ANALYSIS"); 
window(20,6,80,25); 

cprintf("<1>  Setup/  Initialization\r\n"); 
cprintf("<2>  Fin  of  Rectangular  Prof i le\r\n"); 
cprintf("<3>  Fin  of  Trapezoidal  Prof i le\r\n"); 
cprintf("<A>  Fin  of  Triangular  Prof i le\r\n"); 
cprintf(“<5>  Return  to  DOS"); 
return; ) 

void  print_nienu_one(void)( 

c I rscr( ); 
gotoxy( 1,1); 

cprintf  (" . \r\n"); 

cprintfC  <A>  Predict  the  Performancc\r\n"); 

cprintfC  <B>  Analyze  the  Profile  Area\r\n"); 

cprintfC  <C>  Optimize  the  Oimensions\r\n"); 

cprintfC  <D>  Return  to  Main  Menu\r\n"); 

cprintfC . "); 

return;) 

void  print_nienu_two(void)( 

clrscrO; 

gotoxy(1,1); 

cprintfC . \r\n"); 

cprintfC<A>  Admittance  :  %s\r\n",  S1[Y]); 

cprintf("<B>  Area  :  /ls\r\n",  S2(A]); 

cprintfC<C>  Heat  :  %s\r\n",  S3[H]); 

cprintf ("<D>  Impedance  :  %s\r\n",  SA[2]); 

cprintf C<E>  Depth  :  %s\r\n",  S5[L]); 
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cprintf(“<F>  Mass  :  Xs\r\n“,  S6[M]); 

cprintf (‘‘<G>  Tanperature  :  Xs\r\n'‘,  S7tT]); 

cprintf("<H>  Thermal  K  :  Xs\r\n‘*,  S8[K]); 

cprintf{"<I>  Uidth  :  Xs\r\n",  S5[d]); 

cprintf (“<J>  Save  Data  :  Xs\r\n",  S9[S]); 

cprintf {"<IC>  Percent  Error:  Xf\r\n",  error); 
cprintf{"<L>  Return  to  Main  Menu\r\n“); 

cprintf  (" . 

return;) 
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/ 


PERFORMANCE  Analysis  Program 
Longitudinal  Fin  of  Rectangular  Profile 


*/ 


#include 

<stdio.h> 

#include 

<stdl ib.h> 

tinclude 

<math.h> 

lYinclude 

<conio.h> 

#include 

"uni ts.h" 

#def ine 

Yes  1 

#def ine 

No  0 

#def ine 

N  100 

#def i ne 

sigma  5.66961e-8 

char 

file  name [20) ; 

double 

*_<).  B_(),  C_(),  D_(),  m(),  1 

double 

length,  width,  depth; 

double 

inc,  overall  diff  =  100,  Tmax 

double 

qb=0,A1,  B1,“C1,  01,  Ta,  Tb, 

float 

error,  V,  epsilon,  k,  Tbase,  ( 

int 

T,  L,  M,  S,  H,  d,  K,  i,  X,  2, 

int 

overflow,  possible_overf low  = 

void 

inputo,  computeO,  resultsO 

Yb(),  A,  B.  C,  D; 

Tbv_min  =  100000,  data [N+1] [6] ; 
Tave,  temp,  qa,  segment_di f f ; 
ty,  mass,  K2,  eta,  PN; 

nutnber_of_t  r  i  es=0; 


/* 


MAIN  PROGRAM 


************»«**«******»*****«****»»**«*«*«•*****««**«•**•****************  * ! 

int  main(){ 

FILE  •data_in,  *out; 

data  in  =  fopen("unitb.ext","r"); 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  'A  %d  %f“,iT,SL,&M,&S,SH,Sd,&<,SX,S2,&Y,&error); 

fclose(data”in); 

error  =  error  /  100; 

input( ); 

K2  =  2  *  epsilon  *  sigma  *  v  *  pow(Tbase,'i); 

Tave  =  Tbase  •  pow(V,  0.25)  ♦  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 

else 

Ttip  =  .0001  +  Tave; 
inc  =  .5  *  Ttip;  i  =  1; 

while  (fabs(overal l_di f f )  >  (Tbase  *  error)  jj  i  <=  N  )< 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  width  •  epsilon  ♦  pow(Ta,  4)  •  depth; 

while  (i  <■-  N  &&  overflow  ==  No)( 

segment_diff  =  10; 

while  (fabs(segment_diff )  >  ((Tb  •  error)  /  N)  &&  overflow  ==  No){ 

A  =  A_();  B  =  B_();  C  =  C_();  0  =  D_(); 
temp  =A*Ta+B*  qa; 
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qb  =  C  •  Ta  ♦  D  *  qa; 
seginent_dif f  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_diff; 

if  (Tb  >  Tbase  +50  ||  Tb  <  0)< 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_niin  =  Trip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;>> 

if  (overflow  ==  No)< 

datati]  tO]  =  qb;  data[i](1]  =  Tb;  data[i)[2]  =  A; 
data(i][3]  =  B;  datati]  (4]  -  C;  datati]  [5]  =  D; 

Ta  =  Tb;  qa  =  qb;  ♦+];>> 

if  (nuniber_of_tries  >  29)( 

clrscrO;  cprintf ("Fatal  Error;  Unable  to  Confute"); 
getchO;  exit(O);) 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (  Tb  <  Tbase  )  Tmax  =  Ttip; 

overall_diff  =  (Ib  -  Tbase)  /  2; 

if  (  possible_overf low  ==  Yes){ 
inc  =  inc  /  2; 

if  (overal l_di f f  >0)  Ttip=  Ttip  -inc; 
else  Ttip  =  Ttip  ♦  inc;) 

else 

Ttip  =  Ttip  -  overal l_di ff; 

numbor_of_tr i es++; 

if  (Ttip  <  Tave)  Ttip  =  Tave; 

gotoxy(27,24); 

cprintfC'N  =  Xd  Ta  =  %4.1f  Tb  =  %4.1f  ",number_of_tries,w_T(T,Ttip),w_T(T,Tb));>> 
for  (i  =  N  •  1;  i  >=  1;  i--){ 

A1  s  A  •  datati]  (2]  ♦  B  •  datati]  t4];  B1  =  A  •  datati]  13]  ♦  B  •  datati]  tS]; 

Cl  =  C  •  datati]  12]  ♦  D  •  datati]  14);  D1  =  C  *  datati]  13]  +  D  •  datati]  tS); 

A  =  A1;  B  =  B1;  C=C1;  D=D1;} 

if  (S  ==  1)( 

out  =  f open( f i le_name,  "w"); 
for  (i  =  1;  i  <=  N;  i**) 

fprintf(out,  "Xe,  Xe,  Xe\n",  w_L(d,  i  •  length  /  N),  w_T (T .data ti ]  11] ) ,  w  H(H, 
datati]  tO] )  ); 

f closet  out);) 


Ta  =  Ttip; 

mass  =  width  •  length  *  depth  *  density  *  1000; 

PN  =  length  •  sqrt(  sigma  *  epsilon  *  pow(Tbase,3)  /  (  k  •  width)); 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

results( ); 

exit(O); 

return;} 

/• . . . 


OUTPUT  PROCEDURE 


*/ 


void  results(void)( 


clrscrO; 

cprintf("\a  THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE"); 

cprintf  ("\r\n  =  ===  =  irrz=rr  =  =  ==  =  z  =======:  =  r===:ii==  =  =========z:i==m===s  ===  =  ===r  =  ==\r\n"); 


78 


cprintfC  Density  =  %4.2f  kg  /  dm''3\r\n",densi ty); 

cprintfC"  Fin  Depth  =  p_L(L,  depth); 

cprintfC  Total  Mass  =  p_M(M,  mass); 

cprintfC  Fin  Length  =  ");  p~L(L,  length); 

cprintfC  Fin  Width  =  '■);  p_L(d,  width); 

cprintfC  Emissivity  =  %4.2f\r\n",  epsilon); 

cprintfC  Profile  Area  =  ");  p_A(X,  width*length); 

cprintfC  Profile  Number  =  5i4.3f\r\n“,PN); 

cprintfC  Base  Heat  Flow  =  •');  p_H(H,  qb>; 

cprintfC  Fin  Efficiency  =  %A.3f\r\n",eta); 

cprintfC  Base  Temperature  =  '*);  p_T(T,Tb); 

cprintfC  Edge  Temperature  =  ");  p_T(T,Ta); 

cprintfC  Thermal  Conductivity  =  ");  p_k(K,  k); 

cprintfC  Environmental  Parameter  =  5iA.2f\r\n",  V); 

cprintfC  ===============rr======r====r=======s============r=============s==\r\n"); 

cprintfC  Transmission  Parameter  Matrix  \r\n\r\n"); 

cprintfC  A  =  %4.Af  B  =  ",  A);  p_Z(Z,B); 

cprintfC  C  =  ");  p_T(Y,C); 

cprintfC  D  =  %4.Af",  D); 

getchl ); 

return;) 


y*  *********************************************************************** 

INPUT  PROCEDURE 

A*************************************************************************  It  j 


void  input(void){ 

c I rscr< ) ;_setcursortype(_NORMALCURSOR ) ; 


cprintfC\r\n  THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE\n\r"); 

cprintf  C\r\n  =  =  =  =  =  :r  =  zz  =  ii  =  zzz  =  ii  =  z  =  =  =  =  ===zz===r  =  =  =  =====  ===z  =  z  =  ir  =  zi:=ii=r=mi  =  \r\n") 

cprintfC  Fin  Depth  ");  depth  =  i_L(L); 

cprintfC  Fin  Length  ");  length  =  i_L(L); 

cprintfC  Fin  Width  ");  width  =  i_L(d); 

cprintfC  Emissivity  =  ");  scanfC%f",  Sepsilon); 

cprintfC  Density  (  kg  /  dm*3  )  =  ");  scanf l"*4f'',  ^density); 

cprintfC  Base  Temperature  ");  Tbase  =  i_T(T); 

cprintfC  Thermal  Conductivity  ");  k  =  i_K(IO; 

cprintfC  Environmental  Parameter  (0  to  .8)  =  ");  scanfC%f",  SV); 

if  <S  ==  1)  { 

cprintfC  Output  Filename  is  =  "); 

scanfC%s",  file_name);> 

cprintf  C\r\n================= -=======^===============:=========="); 


if  (  epsilon  >  1.0  ]]  epsilon  <=  OC 

cprintf C\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchl);  exit(O);) 

if  (  V  <  0  ; 1  V  >  0.8){ 

cprintf C\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchl);  exitlO);) 

_setcursortypel_NOCURSOR ) ; 
return; > 


global  FUNCTIONS 


•/ 


double  hrlvoidX 

float  Tav  =  iTb  +  Ta)  /  2; 

returnlsigma  •  epsilon  ♦  powlTav,  3)  -  IIC2  /  12  *  Tav)));> 
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double  m(void)( 

return(  sqrt(  (  2.0  *  hrO)  /  (  k  *  width)  )  );) 


double  YbCvoidX 

return(  sqrt{  2.0  *  hr()  *  k 

double  A_{void){ 

return(  cosh<  tnO  *  length  / 

double  B_{void)< 

return  {  sinhC  mO  •  length  / 

double  C_(void)< 

return  (  sinh(  m()  •  length  / 

double  D_(void){ 

returnC  cosh(  mO  *  length  / 


*  width)  *  depth  );> 
N  )  );> 

N  )  /  VbO  );> 

N  )  •  YbO  );) 

N  )  );> 
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y*  ***************!»*************************************************»********* 


PERFORMANCE  Analysis  Program 
Longitudinal  Fin  of  Trapezoidal  Profile 

****************************************************************************  * ! 


#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
^include  <conio.h> 
#include  "bessel.h" 
#include  "units.h" 


#def ine 
#def ine 
#def i ne 
#def ine 


Yes  1 

No  0 

N  100 

sigma  5. 66961 e- 8 


char  f  i  le_name[20] ; 

double  A_(),  B_(),  C_0,  D_(),  m(),  hr(),  uCdouble  x).  A,  B,  C,  D; 

double  length,  depth,  edge_Hidth,  base_width,  fin_taper; 

double  b1 ,  db,  a,  b,  ua,  ub,  Ta,  Tb,  a_width,  b_Hidth; 

double  db=0,  inc,  Tave,  overaU_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000,  data[N+1H6]  ; 

double  A1,  81,  Cl,  01,  taper_ratio,  Ttip,  temp,  qa,  segment_di f f ; 

float  V,  epsilon,  k,  error,  Tbase,  density,  mass,  K2,  Ap,  eta,  PN; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  2,  Y; 

int  overfloy,  possible_overf low  =  No,  number_of_tries=0; 

void  input! ),  results!); 


*********************************************************************** 


MAIN  PROGRAM 


************************************************************************** 


int  mainO< 

FILE  •data_in,  •out; 

data_in  =  fopen!"uni ts.ext","r"); 

fscanf!data_in,"%d  %d  %d  %d  %d  %d  /id  %d  %d  %d  %f",&T,&L,SM,&S,&H,&d,&K,&X,SZ,&Y,&error); 

fclose!data”in); 

error  =  error  /  100; 

input!); 

taper_ratio  =  edge_width  /  base_width; 
b1  =  length  /  !  1  -  taperratio  ); 
fin_taper  =  atanZ!  base_Hidth,  2*b1  ); 
db  =  length  /  N; 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow!Tbase,A); 

Tave  =  Tbase  •  pow!V,  0.25)  ♦  0.0001; 

if  !V==0) 

Ttip  =  .8625  *  Tbase  •  sqrt!  2  *  V  )  +  !  !  Tbase  -  LTbase  •  sqrt!  2  •  V  )))/2); 

else 

Ttip  =  Tave  ♦  0.0001; 
inc  =  .5  •  Ttip;  i=1; 

while  !fabs!overall_diff)  >  !Tbase  *  error)  j j  i  <=  N  ){ 


overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  •  edge_width  *  epsilon  •  pow!Ta,  A)  *  length; 

a  =  b1  -  length; 
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while  (i  <=  N  &&  overflow  ==  No)C 


segment_diff  =  10; 

a_width  =  edge_width  •  a  /  (  b1  -  length); 
b_width  =  edge_width  •(a+db)/(b1-  length); 

while  <fabs(segment_dif f )  >  (  Tb  *  error)  /  N  iS  overflow  ==  No){ 

ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  =  A_{);  B  =  B_();  C  =  C_();  D  =  D_(); 
temp  =  A  •  Ta  ♦  B  •  qa; 
qb  =  C  •  Ta  ♦  0  •  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_dif f; 

if  (Tb  >  Tbase  ♦  50  i;  Tb  <  0)( 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  -Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)> 


if  (overflow  ==  ko)< 

datalillO)  =  qb;  data[i][l]  =  Tb;  dataCi][21  =  A; 
data[i][3]  =  B;  data[i)[4]  =  C;  data[i][51  =  D; 
Ta  =  Tb;  qa  =  qb;  a  =  a  »  db;  ♦«i;)) 


if  (number_of_tries  >  29)( 

clrscr();  cprintf ("Fatal  Error;  unable  to  Compute"); 

getchl);  e^itfO);) 

« 

if  (overflow  ==  ko)( 

i f  (TmaA  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  •  Tbase)  /  2; 

if  (  possible_overf low  ==  Yes)( 
inc  =  inc  /  2; 

if  (overat l_di f f  >  0  )  Ttip  =  Ttip  •  inc; 
else  Ttip  =  Ttip*  inc;) 

else 

Ttip  =  Ttip  -  overal l_diff ; 
number_of _t  r i es+* ; 
gotoxyT27, 24 ) ; 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",number_of_tr i es, w_T (T , T t ip) , w_T(T ,Tb) ); 

if  (Ttip  <  Tave)  Ttip  =  Tave;)} 

for  ( i  =  N  -  1 ;  i  >=  1 ;  i -  - )( 

A1  =  A  *  dat3li][21  +  B  *  data[i)[4];  B1  -  A  *  data[i)[3)  ♦  B  •  dataliltS); 

Cl  =  C  ♦  data[i][2]  +  0  *  data[n[4);  01  =  C  *  data[i)[3]  +  0  *  data[i)[53; 

A  =  A1;  B  =  61;  C=C1;  0=01;) 


If  (S  ==  1)< 

out  =  f open( f i le_name,  "w"); 
for  (i  =  1;  i  <=  N;  i++) 

fprintf(out,  "%e,  %e,  %e\n",  w_L(d,  i  *  length  /  N),  w_T(T,data[i)  [1] ),  w_H(H, 

datali]  [0] )  ); 

fclose(out);) 

T  a  =  T 1 1  p; 

Ap  =  .5  *  length  •  (  edge_width  ♦  base_width  ); 
mass  =  Ap  •  length  *  density  *  1000; 

PN  =  2  *  sigma  •  epsilon  *  pow; length,  2)  *  pow(Tbase,3)  /  (k*  basewidth); 
eta  =  qb  /  (  2  •  sigma  *  epsilon  *  length  *  length  *  pow(Tbase,4)); 
resui ts( ); 
e« ' t(0); 
return; ) 
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/ 


OUTPUT  PROCEDURE 

. . . . . . .  ./ 

void  resul ts(void){ 


clrscrO; 

cprintfC  THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFILE"); 

cpr intf ("\r\n  =======r=================is=====i======i=r ======= ====z===z=z=====r\r\n") 

cprintfC  Density  =  X4.2f  kg  /  dm"3\r\n",densi ty); 

cprintfC  Fin  Depth  =  ");  p_L(L,  depth); 

cprintfC  Total  Mass  =  ");  p_M(M,  mass); 

cprintfC  Fin  Length  =  ");  p_L(L,  length); 

cprintfC  Emissivity  =  %4.2f\r\n",  epsilon); 

cprintfC  Profile  Area  =  ");  p_A(X,Ap); 

cprintfC  Profile  Number  =  5i4.3f\r\n",PN); 

cprintfC  Base  Heat  Flow  =  ");  p_H(H,  qb); 

cprintfC  Fin  Efficiency  =  %A.3f\r\n",eta); 

cprintfC  Fin  Edge  Width  =  ");  p_L(d,  edge_Midth); 

cprintfC  fin  Base  Width  =  ");  p_L(d,  base_width); 

cprintfC  Base  Temperature  =  ");  p_T(T,Tb); 

cprintfC  Edge  Temperature  =  ");  p_T(T,Ta); 

cprintfC  Thermal  Conductivity  =  *');  p_K(K,  k); 

cprintfC  Environmental  Parameter  =  %A.2f\r\n",  V); 

cpri ntf ( "  =============================:====================================\r\n") ; 

cprintfC  Transmission  Parameter  Matrix  \r\nSr\n"); 

cprintfC  A  =  %4.4,f  G  =  ",  A);  p_Z(Z,B); 

cprintfC  C  =  ");  p  Y(Y,C); 

cprintfC  D  =  %4.4f",  D); 

getchC  ); 

return; > 


/*  . . * . . . . . . . . 

INPUT  PROCEDURE 

void  input(void){ 

clrscrl );  _se t cursor type(_NORMAL CURSOR ); 


cprintfC\r\n  THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFlLE\r\n"); 

cpr intf {"\r\n  =========================================================-======== \r\n") 

cprintfC  Fin  Depth  ");  depth  =  i_L(L); 

cprintfC  Fin  Length  ");  length  =  i_L(L); 

cprintfC  Fin  Edge  Width  ");  edge_width  =  i_L(d); 

cprintfC  Fin  Base  Width  ");  base_width  =  i_L(d); 

cprintfC  Emissivity  =  ");  scanfCXf",  Sepsilon); 

cprintfC  Density  (  kg  /  dm'3  )  =  ");  scanfC/if",  idensity); 

cprintfC  Base  Temperature  ");  Tbase  =  i_T(T); 

cprintfC  Thermal  Conductivity  ");  k  =  i_K(k); 

cprintfC  Environmental  Parameter  (0  to  .8)  =  ");  scanfcXf",  &V); 

if  (S  ==  1)  { 

cprintfC  Output  Filename  is  =  "); 

scanfC%s",  file_name);) 

cprintf{"\r\n  =================================================================="); 

if  (  edge_width  >=  base_width)( 


cprintf C\r\n\r\nFATAL  ERROR;  Fin  Edge  Width  >=  Fin  Base  Width"); 
getchO;  exit(O);) 

if  (  epsilon  >  1.0  \\  epsilon  <=  0){ 

cprintf C\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchO;  exit(0);> 
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if  (  V  <  0  II  V  >  0.8)< 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(O);) 

_setcursortype(_NOCURSOR); 
return; > 

y*  A********************************************************************** 


GLOBAL  FUNCTIONS 


double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

returnlsigma  *  epsilon  •  powCTav,  3)  -  CK2  /  (2  *  Tav)));) 
double  mCvoidX 

returnC  hr()  /  (k  *  sin( f i n_taper ) ) ) ; ) 

double  uCdouble  x)( 

return(  2  *  sqrt(x  *  m( )  )  );) 

double  A_(void){ 

return(  ua  *  ( 1 1 (ua)*K0(ub)  ♦  1 0(ub)*Kl (ua) )  );) 

double  B_(void)<. 

return  (  ub*ub*( I0(ub)*K0(ua)- I0{ua)*X0(ub))  /  (  2*k*b_width*depth*m( )));) 
double  C_(void)< 

return  (  2*k*a_width*depth*m( )’ub*( IT (ub)*Kl (ua)  -  1 1 (ua)*K1 (ub) )/(ua)  );) 
double  0  (voidX 

return(  ub*<  1 0<ua)»k1 (ub)+ 1 1 (ub)»k0{ua) )  );) 
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PERFORMANCE  Analysis  Program 
Longitudinal  Fin  of  Triangular  Profile 


(Kinclude  <stdio.h> 
#include  <stdl ib.h> 
^include  <math.h> 
ffinclude  <conio.h> 
^include  "bessel.h" 
#include  "uni ts.h" 


iddefine  Yes  1 

lYdefine  No  0 

#define  N  99 

#define  sigma  5.66961e-8 

char  f i le_name [20] ; 

double  A_(),  B_(),  C_(),  D_{),  m( ) ,  hr(),  ufdouble  x).  A,  B,  C,  D; 

double  length,  depth,  edge_Hidth,  base_Hidth,  fin_taper; 

double  b1 ,  db,  a,  ua,  ub,  Ta,  Tb,  a_width,  b_width; 

double  qb=0,  inc,  Tavc,  overall  dtff  =  100,  Tmax  =  0,  Tbv_min  =  100000,  data[N+1][6] 

double  A1,  B1,  Cl,  01,  Ttip,  temp,  qa,  segment_dif f ; 

float  V,  epsilon,  k,  error,  Tbase,  density,  mass,  K2,  Ap,  PN,  eta; 

int  T,  L,  M,  S,  H,  d,  k,  i,  k,  2,  Y; 

int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

void  inputC),  resultsl); 


MAIN  PROGRAM 


int  rnainlX 

FILE  •data_in,  *out; 

data  in  =  f openCuni ts .ext" , "r" ) ; 

fscanf(data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f",ST,fiL,iM,SS,&H,Sd,&K,&X,SZ,iY,Serror>; 

fclose(data_in); 

error  =  error  /  100; 

inputl ); 

db  =  length  /  (  N  ♦  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  basewidth,  2*bl  ); 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbasc,L); 

Tave  =  Tbase  *  pow(V,  0.25]  ■»  0.0001; 

if  (V==0) 

Ttip  =  .8625  •  Tbase  *  sqrt(  2  *  V  >  ■*  L  (  Tbase  -  Ubase  *  sqrtT  2  *  V  )))/2); 

else 

Ttip  =  Tave  ♦  0.0001; 
inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overall_diff )  >  (  Tbase  *  error)  j)  i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
segment_diff  =  10; 

while  C f abs(segment_di f f )  >  (  Tb  •  error)  /  N  S&  overflow  ==  NoX 
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ub  =  u(db}; 

temp  =  Ta  /  lO(ub); 

C  =  2  *  k  *  edge_Hidth  •  depth  *  pow(m(),2)  *  IKub)  /  (ub  *  lO(ub)); 
qb  =  C  *  temp; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  ♦  segment_di f f ; 

if  (Tb  >  Tbase  +50  )|  Tb  <  0){ 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  ■  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;>) 

if  (overflow  ==  No){ 

datatiJCO]  =  qb;  dataliHl]  =  Tb;  data[i][2]  =  1; 
data[i][3]  =  0;  data[i3[41  =  C;  datali]  (53  =  1; 

Ta  =  Tb;  qa  =  qb;  ++i;) 

a  =  b1  -  length; 

while  (i  <=  k  &&  overflow  ==  ko){ 
segment_diff  =  10; 

a_width  =  edgewidth  *  a  /  (  bl  -  length); 
b_width  =  edge_width  *(a+db)/(b1-  length); 

while  ( fabs(segment_di t f )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No)< 

ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  =  A_();  B  =  6_0;  C  =  C_();  0  =  D_(); 
t  emp  =  A*Ta  +  8*qa; 
qb  =C*Ta»D*  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_dif f ; 

if  (Tb  >  Tbase  +  50  1 1  Tb  <  OX 

overflow  =  Yes;  poss>ble_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  •  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  •  inc); 
if  (Ttip  <  TaveT  Ttip  =  Tave;)} 

if  (overflow  ==  No){ 

data[il[0]  =  qb;  data(i3(1)  =  Tb;  datati][2)  =  A; 
data(il(3]  =  B;  data[i)t43  =  C;  datali]  (5)  =  D; 

Ta  =  Tb;  qa  =  qb;  a  =  a  +  db;  ++i;>> 

if  (number_of_tries  >  29)( 

clrscrO;  eprintf ("Fatal  Error:  Unable  to  Compute”); 
getchO;  exit(O);) 

if  (overflow  ==  NoX 

i f  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overalldiff  =  (Tb  -  Tbase)  /  2; 

if  (  possi bleoverf low  ==  Yes)( 
me  =  me  /  2; 

if  (overal l_di f f  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  +  inc ;> 

else 

Ttip  =  Ttip  -  overal l_di ff ; 
number_of _t  r i es++  ; 
gotoxy(27,2A); 

cprintfC'k  =  %d  Ta  =  %4.1f  Tb  =  %A.1f  ",nuniber_of_tries,w_T(T,Ttip),w_T(T,Tb)) 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 
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for  ( i  =  N  -  1;  i  >=  1;  i -  - ){ 

A1  =  A  •  data[i]t2]  ♦  B  •  data[i][A];  B1  =  A  *  datati]  [3]  ♦  B  •  datati]I5]; 
Cl  =  C  *  datati]  [2]  +  D  *  datati]  t4];  D1  =  C  *  datati]  t3]  ♦  D  *  datati]  t5]; 
A  =  A1;  B  =  B1;  C=C1;  D=D1;) 

length  =  length  +  db; 


if  (S  ==  1){ 

out  =  fopenC f i le_name,  "w"); 
for  <i  =  1;  i  <=  N;  i»+) 

fprintfCout,  "%e,  Xe,  %e\n*',  w_L(d,  i  *  length  /  N),  w  T(T, datati]  t1] ),  h  H(H, 

datati]  tO]  )); 

fclose(out);) 

Ta  =  Ttip; 

Ap  =  .5  *  length  •  base_width; 
mass  =  Ap  *  depth  *  density  *  1000; 

PN  =  2  *  sigma  *  epsilon  *  powl length, 2)  *  pon(Tbase,3)  /  (k*base_width>; 
eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  •  depth  *  pow(Tbase,4)); 
resul ts( ) ; 
ex i t ( 0 )  ; 
return; > 


OUTPUT  PROCEDURE 


*/ 


void  resul tsCvoi d)( 


clrscrC ); 

cprintf("\a  THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROF  I LE \r\n" ) ; 

Cprintf  ("  =  =  =  =  =  =  :  =  =  =  =  =  s  =  z  =  =  =  r  =  =  =  =  =  =  =  =  m  =  = 

cprintfC  Density  =  %4.2f  kg  /  dnr3\r\n",dens i ty) ; 

cprintfC  Fin  Depth  =  ");  p_L(L,  depth); 

cprintfC  Total  M..SS  =  ");  p~M(M,  mass); 

cprintfC  Fin  Length  =  ");  p”l(L,  length); 

cprintfC  Emissivity  =  %4.2fV\r>",  epsilon); 

cprintfC  Profile  Area  =  ");  p_A(X,Ap); 

cprintfC  Profile  Number  =  %4.3f \r\n",PN); 

cprintfC  Base  Heat  Flow  -  ");  p_H(H,  qb); 

cprintfC  Fin  Efficiency  =  %A.3f\r\n",eta); 

cprintfC  Fin  Base  Width  =  ");  p_L(d,  base_width); 

cprintfC  Base  Temperature  =  ");  p_T(T,Tb); 

cprintfC  Edge  Temperature  =  ");  p_T(T,Ta); 

cprintfC  Thermal  Conductivity  =  ");  p_K(K,  k); 

cprintfC  Environmental  Parameter  =  /!«.2f\r\n",  V); 

cprintfC  zzzzzzzzzzzzzzz=zz:zzzzzzzzzzzzzzzzzzrzzzzzzzzzzzzzzzzzzzzzzzzzzzz\r\n"); 

cprintfC  Transmission  Parameter  Matrix  \r\n\r\n"); 

cprintfC  A  =  6  =  ",  A);  p_2(2,B); 

cprintfC  C  =  ");  p_Y(Y,C); 

cprintfC  D  =  %4.4f",  D); 

getchC ); 

return; ) 


/* 


INPUT  PROCEDURE 


*/ 


void  input(void)< 

clrscrC );  _setcursor type{_NORMAL CURSOR); 


cprintfC\r\n  THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROFILE\r\n"); 

cpr intf ("\r\n  zzzrzzzrzzzzzzzzzzzzzzzzzzzzzzzzzzrzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzz\r\n"); 
cprintfC  Fin  Depth  ");  depth  =  i_L(L); 

cprintfC  Fin  Length  ");  length  =  i_L(L); 
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cprintfC  Fin  Base  Width  •');  base_width  =  i_L(d); 

cprintfC  Emissivity  =  scanf("%7",  Sepsilon); 

cpnntfC  Density  (  kg  /  dm*3  )  =  ");  scanf("%f",  idensity); 

cprintfC  Base  Temperature  Tbase  =  i_T(T); 

cp''intfC  Thermal  Conductivity  k  =  i_k(K); 

cprintfC  Environmental  Parameter  (0  to  .8>  =  '■>;  scanfCitf",  &V); 

if  (S  ==  1)  { 

cprintfC  Output  filename  is  = 

scanfCXs",  f  i  le_name);> 

cpr  intf  C\r\n  ===============  ========;=;======z=:=======m========i===============ii); 


if  (  epsilon  >  1.0  ||  epsilon  <=  0)( 

cprintf C\r\n\r\nFATAL  ERROR;  Emissivity  is  not  between  0  and  1"); 
getchO;  exit(O);) 

if  (  V  <  0  1 1  V  >  0.8){ 

cprintf C\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8") 
getchO;  exit(O);) 

edge_width  =  base_width  /  1e2; 

_setcursortype{_NOCURSOR ) ; 
return; > 

/.  . . . . . . 

GLOBAL  FUNCTIONS 

. . . . . . . . . 

double  hr(void){ 

float  Tav  =  (Tb  ♦  Ta)  /  2; 

returnCsigma  *  epsilon  *  powCTav,  3)  •  (k2  /  (2  *  Tav)));} 
double  m(void){ 

returnl  hr()  /  (k  *  si n< f in_tapor ) ) ) ; ) 

double  uCdouble  x}< 

returnC  2  *  sqrt(x  ’  m( )  )  );) 

double  A  (voidX 

returnl  “ua  *  ( 11  (ua)*K0(ub)  +  1 0(ub)*pCl (ua) )  );) 

double  B_(void){ 

return  (  ub*ub*( I0(ub)*K0(ua)- I0(ua)*k0{ub))  /  (  2*k*b_width*depth*m( ) ) ); ) 
double  C_lvoid){ 

return  (  2*k*a_width*depth*m(  )*ub*(  I  Uub)*Kl  (ua)  -  n(ua)*K1(ub))/(ua)  );> 
double  0_(void){ 

returnl  ub*(  10(ua)*k1 (ub)»l 1 {ub)*K0(ua))  );) 


/* . . . . . . 

profile  area  Analysis  Program 
Longitudinal  Fin  of  Rectangular  Profile 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  <conio.h> 

#include  "units. h" 

Hfdefine  Yes  1 

•define  No  0 

•define  N  100 

•define  sigma  5.66961e-8 

char  f i le_name [20] ; 

double  length,  width,  depth,  start ing_width,  ending_width,  width_inc; 

double  A,  B,  C,  0,  Ta,  Tb,  qb=0; 

float  V,  K2,  epsilon,  k,  Tbase,  density,  mass,  Ap,  error,  eta; 

int  T,  L,  M,  S,  H,  d,  tC,  i,  X,  2,  Y; 

void  input{),  computeC); 

/.  . . . . . . . . . 

MAIN  PROGRAM 

. . . . . . . . .  ./ 

int  main(){ 

FILE  *data_in,  ‘out; 

data  in  =  fopen<"units.ext","r"]; 

fscanf(data_in,"%d  %d  %d  Xd  Xd  Xd  Xd  Xd  Xd  XM  5;f",iT,SL,SM,SS,&H,id,&iL,£:x,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

inputC ); 

width_inc  =  (ending_width  -  start i ng_width  )  /  100; 
width  =  starting_width; 

while  (width  <  ending_width){ 

compute( ); 

while  (qb  ==  0){ 

width  =  1 .001  *  width; 
compute( ) ; ) 

gotoxy(19,23); 

cprintf ("width  =  Xit.M  length  =  %A.1f  heat  =  %A.1f",  w  L(d, width),  w  L(L, length), 
w_H(H,  qb)); 

if  (S  ==  1){ 

out  =  fopen(f i le_name,  "a"); 

fprintf(out,  "%e,  %e,  %e,  %e,  %e,  %e,  %e,  Xe,  Xe,  Xe\n",  w_L(d, width), 

w_T(T,Tb),  w_T(T,Ta),  eta,  w_H(H,  qb),  w_M(M,mass),  A,  w_Z(Z,B),  w_Y(Y,C),  0); 
f c lose(out ) ; ) 

width  =  width  ♦  width_inc;) 

gotoxy(21 ,24);clreol ( );cpr intf ("\aCompleted  \r\n"); 

clreolO; 

getch( ); 
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ex)t(O); 
return  0;> 

/•  . . . . . 

COMPUTE  ROUTINE 

. . . . . . . . .  ./ 

void  computeCvoidX 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  YbO; 

double  inc,  overall_diff  =  100,  Tmax  =  0,  Tbv_nijn  =  100000,  data [N+1]  [6] ; 

double  A1,  Bl,  Cl,  01,  Ttip,  temp,  qa,  Tave,  segment_di f f ; 

int  overflow,  poss ible_overf low  =  No,  number_of_tries=0; 

length  =  Ap  /  width; 

IC2  =  ?  *  epsilon  *  sigma  *  V  *  pow(Tbase,4i); 

Tave  =  Tbase  *  powlV,  0.25)  *  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  ♦  {  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2) 

else 

Ttip  =  .0001  ■»  Tave; 
inc  =  .5  *  Ttip;  i  =  1; 

while  (fabsCoveral l_di f f )  >  (  Tbase  *  error)  [ |  i  <=  N  )( 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  •  widtn  ’  epsilon  *  powCTa,  A)  *  depth; 

while  (i  <=  N  &&  overflow  ==  No)( 

segment_diff  =  10; 

while  (fabs(segment_dif f )  >  (  Tb  *  error)  /  N  overflow  ==  No)( 

A  =  A_();  B  =  B_();  C  =  C_();  0  =  D_(); 
temp  =A*Ta«B*  qa; 
qb  :C*Ta»D*  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  •►50  |  |  Tb  <  0)( 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_mir'  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  ire); 
if  (Ttip  Tave)  Ttip  ^  Tave;)) 

if  (overflow  ==  No)< 

datati][0]  =  qb;  data[i)[1)  =  Tb;  data[i][2]  =  A; 
data[i)t3]  =  B;  data[ij(6]  =  C;  data[i)(5]  =  D; 

Ta  =  Tb;  qa  =  qb;  ■►■►i;)) 

if  (nuniber_of_tries  >  29)< 
qb=d^;  return; ) 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (Tb  s  Tbase  )  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf tow  ==  Yes){ 
me  =  inc  /  2; 

if  (overal l_di f f  >  0  )  Ttip  =  Ttip  -  inc; 
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else  Ttip  =  Trip  +  inc;) 

else 

Ttip  =  Ttip  ■  overall_diff; 
number_of _t r i es++ ; 

90toxyT27,24); 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  544. If 
if  (Ttip  <  Tave)  Ttip  =  Tave;>) 


" , number_of _t  r i es , w_T ( T , T  t i p) , w_T ( T , Tb) ) 


for  (i  =  N  -  1;  i  >=  1;  i--){ 

AT  =  A  •  data[i][2]  +  B  *  data[i][4];  B1  =  A  *  data[i)[3]  +  B  •  data[i][5]; 
Cl  =  C  *  data[i][2]  ♦  D  •  data[i)[4];  D1  =  C  *  data[i][3)  +  D  *  datatiJCS); 
A  =  A1;  B  =  B1;  C=C1;  D=D1;) 


Ta  =  T  t i p; 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

mass  =  width  *  length  •  depth  •  density  *  1000; 

return;) 

y *  *********************************************************************** 


INPUT  PROCEDURE 


•/ 


void  input(void){ 

c I rscrc ) ;  _set cur sortypc(_NORMAL CURSOR ) ; 


cprintf ("XrVn  THE  LONGITUDINAL  TIN  Of  RECTANGULAR  PROFILE"); 

cpr  intf  <"\r\n  =;  =  =  =  =  =  =  =  =  =  i  ===  =  =  =  ==i==  =====  =  =  =  =  =  =:  ===r======  =  =  =  :::  =  =  =  =  =  =  i========:=  =  \r\n"); 

cpr\ntf("  fin  Depth  ");  depth  =  i_L(L); 

cprintfC"  Starting  fin  Width  ");  start ing_width  =  i_L(d>; 

cprintf("  Ending  fin  Width  ");  ending_width  =  i_Lrd); 

cprintfC  Profile  Area  ");  Ap  =  i_ACX); 

cprintfC  Emissivity  =  '*);  scanfC%f",  Sepsilon); 

cprintfC  Density  (  kg  /  dm*3  )  =  ");  scanfC5,f",  idensity); 

cprintfC  Base  Temperature  ");  Tbase  =  i_T(T); 

cprintfC  Thermal  Conductivity  k  =  i_K(K); 

cprintfC  Environmental  Parameter  (0  to  .9)  =  ");  scanfCXf",  &V); 

if  (S  ==  1)  { 

cprintfC  Output  Filename  is  =  ");  scanfC/Cs",  file_name);} 

cpr  intf  ("\n  =  =  =  i  =  =  =  =  =  =  =  =  =  Z  =  =  =  I  =  =  sir  =  =  =  zz  =  =  =  =  =  ::=r  =  ri  =  =  =  =  =  =  ===  =  =  =  =  ism  =  =  =  =  =  i  =  =  \r\n"); 


if  (  epsilon  >  1.0  ]]  epsilon  <=  0)< 

cprintf C\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchO;  exitCO);) 

if  (  V  <  0  1 1  V  >  0.8){ 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchC);  exit(O);) 

_setcursortype(_NOCURSOR )  ; 
return;) 


GLOBAL  FUNCTIONS 


•/ 


double  hr(void){ 

float  Tav  =  (Tb  ♦  Ta)  /  2; 

return(sigma  •  epsilon  •  powCTav,  3)  ■  (tC2  /  (2  *  Tav)));) 
double  m(void)< 

returnC  sqrt(  (  2.0  *  hr())  /  (  k  *  width)  )  );) 
double  Yb(void){ 
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return(  sqrt(  2.0  •  hr()  *  k  •  width)  *  depth  );) 
double  A_(void)< 

returnC  cosh(  mO  *  length  /  N  )  };} 

double  B_(void)< 

return  (  sinh<  m()  •  length  /  N  )  /  Yb()  );) 
double  C_(void)< 

return  (~sinh(  m()  *  length  /  N  )  *  YbO  );> 
double  D_(void){ 

return(  cosh(  mO  *  length  /  N  )  );) 
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PROFILE  AREA  Analysis  Program 
Longitudinal  Fin  of  Trapezoidal  Profile 

****************************************************************************  -k j 

#include  <stdio.h> 

#include  <stdlib.h> 

Uinclude  <math.h> 

#include  <conio.h> 

#include  "units. h" 

Oinclude  "bessel.h" 

#define  Yes  1 

(Ydefine  Ro  0 

Adefine  N  100 

#define  sigma  5.66961e-8 

char  f i le_name [20] ; 

double  length,  edge_width,  basewidth,  a_width,  b_Hidth,  depth,  starting_width,  ending_Hidth 

double  A,  B,  C,  0,  Ta,  Tb,  qb=0,  ua,  ub,  fin_taper,  a,  width_inc; 

float  V,  K2,  epsilon,  k,  Tbase,  taper_ratio,  density,  mass,  Ap,  error,  eta; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  2,  Y; 

void  inputC),  computef); 


int  mainCX 

file  *data_in,  ’out; 

data_in  =  f open( "uni ts . e» t" , "r" ) ; 

f scanf (data_i n, "5id  %d  %d  %a  %d  %d  *40  %d  i;a  %d  xf ",iT ,SL ,£M,&S,&H,Sd, SK,&X,&2 , SY,Serror); 

fclose(data_in); 

error  =  error  /  100; 

i nput ( ); 

Hidth_inc  =  (ending_width  •  start ing_width  )  /  100; 
base_width  =  start i ng_width; 

while  (base_width  <  ending_width)( 

computeC ) ; 

while  (qb  ==  0){ 

base_width  =  1.001  *  base_width; 
computeC ); ) 

gotoxy(19,23); 

cprintf ("width  =  %4.1f  length  =  %4.1f  heat  =  %A.lf",  w_L(d,base  width),  w_L(L, length), 
w_H(H,  qb)); 

if  (S  ==  1){ 

out  =  f openC f i le_name,  "a"); 

fprintfCout,  "Xe,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e,  %e\n",  w  L(d,base_width), 

w_T(T,Tb),  w_T(T,Ta),  eta,  w_H(H,  qb),  w_M(M,mass),  A,  w_Z(Z,“B),  w_Y(Y,C),  D); 
fclose(out);) 

base_width  =  base_width  +  width_inc;) 

gotoxyCZI ,24 ) ;c I reol C ) ;cpr i nt f ( "\aCompleted  \r\n"); 
clreoU); 
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getchC ); 
exit(O); 
return  0;) 

y*  *********************************** ********************************* 

COMPUTE  ROUTINE 


void  computeCvoidX 

double  A_(),  B_{),  C_(),  D_(),  m(),  hr{),  uCdouble  x); 

double  len9th_Tnc,  b1 ,  db; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000, 

double  Al ,  Bl,  Cl,  01,  Ttip,  temp,  qa,  se9nient_diTf ; 

int  overflow,  possibte_overf low  =  No,  number_of_tries=0; 

edge_width  =  base_width  •  taper_rat io; 
length  =  2  •  Ap  /  (  base_width  *  edge_width  ); 
bl  =  length  /  (  1  -  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2*b1  ); 
db  =  length  /  N; 

K2  =  2  *  epsilon  *  sigma  *  V  *  powCTbasc.v); 

Tave  =  Tbase  *  powCV,  0.25)  *  0.0001; 

If  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  ♦  (  (  Tbase  -  (Tbase  *  sqrt( 

else 

Ttip  =  T ave  »  0 . 0001  ; 
inc  =  .5  *  Ttip;  i=l; 

while  ( f abs(overal l_di f f )  >  (  Tbase  *  error)  ||  i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  ed9e_width  •  epsilon  *  pow(Ta,  A)  ■*  depth; 

a  =  bl  •  length; 

while  (i  <=  N  ii  overflow  ==  No){ 
segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  bl  •  length); 
b_width  =  edge_width  *(a»db)/(b1-  length); 

while  ( f  abs(segnient_di  f  f )  >  <  Tb  *  error)  /  N  &&  overflow  ^ 

ua  =  u(a)  ;  ub  =  u(a  ♦  db); 

A  =  A_();  B  =  B_();  C  =  C_();  0  =  D_(); 
temp  =A*Ta+B*qa; 
qb  =C*Ta'»D*  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  »  Tbase  *50  j]  Tb  <  0){ 

overflow  =  Tes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbvmin  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)} 

if  (overflow  ==  No){ 

dataliJIO]  =  qb;  datati](1)  =  Tb;  datali)  [2]  =  A; 
datainiS]  =  B;  dataCi]  [A)  =  C;  dataCi]  [5)  =  0; 

Ta  =  Tb;  qa  =  qb;  a  =  a  +  db;  ♦♦i;)) 

if  (number_of_tries  >  29)( 
qb=0;  return;) 


•/ 


data[N+1]  [6]; 


2  *  V  )))/2); 


=  No)( 
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if  (overflow  ==  No)( 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overaU_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  Yes){ 
inc  =  inc  /  2; 

if  (overal l_di f f  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  +  inc;} 

else 

Ttip  =  Ttip  -  overal l_di ff ; 
nu^^be^_of_t^ies♦♦; 

90toxyT27,24); 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",nijmber_of_tries,w_T(T,Ttip),w_T(T,Tb)) 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

for  ( i  =  N  -  1 ;  i  >=  1 ;  i -  - ){ 

A1  =  A  *  datali]  [2]  +  B  *  datali]  (41;  B1  =  A  •  datali]  [3]  +  B  •  data[i)  [53; 

C1  =  C  •  datali]  [23  ♦  0  *  datali]  [4];  D1  =  C  •  datali]  [3]  +  D  *  datali]  153; 

A  =  Al;  B  =  Bl;  C=C1;  D=D1;} 

Ta  =  Ttip; 

mass  =  .5  *  (  edge_width  ♦  basc_widtri  )  *  length  *  depth  *  density  *  1000; 
eta  =  qb  /  (  2  *  sigma  *  epsilon  ’  length  *  depth  ’  pow(Tbase,  4)); 
return; > 

/*  . . . . . . . 


INPUT  PROCEOURE 


void  input(votd)( 


cl rscr( );  _se t cur sort ype(_NORMAL CURSOR); 


cprintf (" 
cprintf ("\r\n 
cprintf (" 
cprintf  (" 
cprintf  (" 
cprintf  (" 
cprintfC 
cprintf  (" 
cprintfC 
cprintfC 
cprintfC 
cprintfC 


THE  longitudinal  fIN  OF  TRAPEZOIDAL  PROFILE"); 


Fin  Depth  ");  depth  =  i_L(L); 

Starting  Base  Fin  Width  ");  start i ng_wi dth  =  i_L(d); 

Ending  Base  Fin  width  ");  endin9_width  =  i_LTd); 

Taper  Ratio  =  ");  scanf("%f",  &taper_ratio); 
Profile  Area  ");  Ap  =  i_A(X); 

);  scanfC%f",  Sepsilon); 

");  scanfC/if",  Sdensity); 


=\r\n“); 


Emissivity  = 

Density  (  kg  /  dm'3  )  = 

Base  Temperature  ");  Tbase  =  i_T(T); 

Thermal  Conductivity  ");  k  =  i_K(k); 

Environmental  Parameter  (0  to  .8)  =  ");  scanfCXf",  &V); 


cprintfC  Output  Filename  is  =  ");  scanf("7.s",  file_name);} 

cprintf  ("\n  =  =  =  =  =  =  =  =  =  r=:=;  =  =  :  =  i  =  =  =  =  =  =  z  =  z  =  :r  =  =  =-r  =  =  =  =  =  =  =  :  =  =  =  =  =  =  =  =  =  z  =  =  =  =  =  =  =  s  =  rr  =  =  =  \r\n"); 


if  (  taper_ratio  >  .99  ||  taperratio  <=  0.0)  ( 

cprintf C\r\n\r\nFATAL  ER^ROR:  Taper  Ratio  is  not  between  0  and  1"); 
getchO;  exit(0);) 

it  (  start ing_width  >=  ending_width)( 

cprintf ("\r\n\r\nPATAL  ERROR;  Starting  Base  Fin  Width  >=  Ending  Base  Fin  Width"); 
getchO;  exit(O);) 

if  (  epsilon  >  1.0  ||  epsilon  <=  0)( 

cprintf ("\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchO;  exit(O);) 

if  (  V  <  0  I  I  V  >  0.8)( 

cprintf C\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(O);) 
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_setcursortype(_NOCURSOR); 

return;) 

GLOBAL  FUNCTIONS 


*/ 


double  nr(void}{ 

float  Tav  =  (Tb  ♦  Ta)  /  2; 

returnCsigma  •  epsilon  *  pou(Tav,  3)  -  (K2  /  (2  *  Tav)));} 
double  mCvoidX 

returnC  hrO  /  (k  *  sin{f in_taper)));) 

double  uCdouble  x){ 

returnC  2  •  sqrt(x  •  mO  )  );> 

double  A  (voidX 

return(  ~ua  •  (11(ua)*K0(ub)  ♦  I0(ub)»KXua))  );) 
double  B_(voidX 

return  ub*ub*< I0(ub)*K0(ua)- I0(ua)*k0(ub))  /  {  2*k*b_width*depth*m()));> 
double  C_(voidX 

return  (  2*k*a_width*depth*m( )*ub*( 1 1 (ub)*Kl (ua)  -  l1(ua)*K1(ub))/(ua)  );) 
double  D_(voidX 

returnC  ub*<  I0(ua)*Kl (ub)*l 1 (ub)*kO<ua) )  );> 
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PROFILE  AREA  Analysis  Program 
Longitudinal  Fin  of  Triangular  Profile 

. . . . ! 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <inath.h> 

#include  <conio.h> 
ffinclude  "units. h“ 

#include  “bessel.h" 

#define  Yes  1 

#define  No  0 

#def i ne  N  99 

#define  sigma  5.66961e-8 

char  f i le_name [20] ; 

double  length,  width_edge,  base_Htdth,  a_uidth,  b_width,  depth,  starting_Hidth,  ending_Hidth; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  ua,  ub,  fin_taper,  a,  inc_width; 

float  V,  K2,  epsilon,  k,  Tbase,  density,  mass,  Ap,  error,  eta; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  2,  Y; 

void  input! ),  computel); 

/*  . . 

MAIN  PROGRAM 

. . ./ 

int  main(){ 

FILE  •data_in,  ‘out; 

data_in  =  fopenr'units.exi'',''r’'); 

fscanf(data_in,"%d  %d  'M  %d  %d  %d  !4d  %d  %d  %d  %f",ST,SL,SM,SS,SH,&d,&K,SX,S2,SY,Serror); 

fclose(data_in); 

error  =  error  /  100; 

inputO; 

inc_width  =  (ending_width  -  start ing_width  )  /  100; 
base_width  =  starting_width; 

while  (base_width  <  ending_widthj{ 

computeO; 

while  (qb  ==  0)( 

base_width  =  1.001  *  base_width; 
computeO;) 

gotoxy(19,23);  cprintf ("width  =  %4.1f  length  =  XA.If  heat  =  X4.1f",  w_L(d,base_width), 
w_L(L, length),  w_H(H,  qb)); 

if  (S  ==  1){ 

out  =  fopen(f i le_name,  "a"); 

fprintf(out,  "%e,  %e,  Xe,  %e,  %e,  Xe,  Xe,  Xe,  Xe,  Xe\n",  w_L(d,base_uidth), 

w_T{T,Tb),  w_T(T,Ta),  eta,  w_H(H,  qb),  w_M(M,mass),  A,  w  ZCzTB),  w  Y(Y,C),  D); 
fclose(out);} 

base_width  =  base_width  ♦  inc_width;) 

gotoxy(21 ,2A);clreol( );cprintf ("\aCompleted  \r\n"); 

clreolO; 

getchO; 
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exjt(0}; 
return  0;> 

/•  •••••** 


COMPUTE  ROUTINE 


•/ 


void  conipute(void){ 

double  A_(),  B_(),  C_(),  0_(),  tn(),  hr(),  uCdouble  x); 

double  length_Tnc,  bT,  db; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000,  data[N+1] [6] ; 

double  Al,  B1,  Cl,  01,  Ttlp,  tenp,  qb,  segment_di f fj 

int  overflow,  possible_overf low  =  No,  njiiber_of_tries='0; 

length  =  2  •  Ap  /  base_width; 

width_edge  =  base_width  /  100; 

db  =  length  /  (  N~+  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*b1  ); 

K2  =  2  •  epsilon  •  sigma  •  V  *  pow(Tbase,4); 

Tave  =  Tbase  *  powlV,  0.25)  ♦  0.0001; 

if  (V==0) 

.  Ttip  =  .8625  *  Tbase  *  sqrtC  2  *  V  )  ♦  (  (  Tbase  -  (Tbase  •  sqrtC  2  *  V  )))/2); 

else 

Ttip  =  Tave  ♦  0.0001 ; 
inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overaU_diff )  >  (  Tbase  •  error)  |  i  i  <=  N  X 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
se9ment_diff  =  10; 

while  (fabs(segment_diff)  >  <  Tb  •  error)  /  N  SS  overflow  ==  No){ 

ub  z  u(db); 

temp  =  Ta  /  lO(ub); 

C  =  2  •  k  •  width  edge  •  depth  *  pow(n(),2)  •  iKub)  /  (ub  •  lO(ub)); 
qb  =  C  •  Tb; 

segnient_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_diff; 

if  (Tb  >  Tbase  +50  jj  Tb  <  0)( 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  TaveT  Ttip  =  Tave;)) 

if  (overflow  ==  No)( 

data[ij[0]  =  qb;  data[i)C1)  =  Tb;  data(i)[2]  =  1; 
data[i][3]  =  0;  data[i][4]  =  C;  datati)t5)  =  1; 
a  =  b1  -  length; 

Ta  =  Tb;  qa  =  qb;  ++i;) 


while  (i  <=  N  &&  overflow  =:=  No)( 
segiiient_di  f  f  =  10; 

a_width  =  width_edge  •  a  /  (  b1  -  length); 
b_width  =  width“edge  •  (  a  ♦  db  )  /  (  b1  -  length); 

while  (fabs(segfnent_diff )  >  (  Tb  *  error)  /  N  &&  overflow  ==  No)( 
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ua  =  u(a)  ;  ub  =  u(a  *  db); 

A  =  A_();  8  =  B_();  C  =  C_();  D  =  D_(); 
tentp  =  A  *  Ta  +  B  •  qa; 
qb  =  C  •  Ta  ♦  D  *  qa; 
seg(nent_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  terfip  ♦  segment_diff ; 

if  (Tb  >  Tbase  +  50  Tb  <  0){ 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_mTn  =  Ttip; 
inc  =  (Tbv_min  •  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  TaveT  Ttip  =  Tave;>> 


if  (overflow  =-  No)( 

data[i)[0)  =  qb;  data[i][1]  =  Tb;  data[i)[2]  =  A; 
dataCi]l3J  =  B;  data[i)(41  =  C;  datati)  [5)  =  D; 
Ta  =  Tb;  qa  =  qb;  a  =  a  ♦  db;  ♦♦i;» 

if  (nuiiber_of_tries  >  29)< 
qb=0;  Teturn;) 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 


overaU_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possi bte_overf low  ==  Yes){ 

1nc  =  inc  /  2; 

if  (overal l_diff  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  ♦  inc;) 

else 

Tt\p  =  Ttip  •  overaU_diff; 
number_of_tries+'f; 
gotoxyT27,24); 

cprintfC'N  *  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",nutnber_of_tries,w_T(T,Ttip),w_T(T,Tb)) 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 


for  (i  =  N  •  1;  i  >=  1;  i--)( 

A1  =  A  •  datati)  12)  ♦  B  *  datati]  14);  B1  =  A  *  datati)  t3)  ♦  B  *  datati)  t5); 
Cl  =  C  •  datati)  12)  ♦  D  *  datati)  t4);  D1  =  C  *  datati)  t3)  +  0  *  datati)  15); 
A  =  A1;  B  =  B1;  C=C1;  D=01;) 


length  =  length  +  db; 

Ta  =  Ttip; 

mass  =  .5  *  base_width  •  length  *  depth  *  density  *  1000; 

eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 

return; ) 

***♦♦*♦***#******♦*********•••••*#*******♦*******♦****♦♦♦********»***** 


INPUT  PROCEDURE 


*/ 


void  input(void)( 


clrscrt ); 
cprintf (' 
cprintf (' 
cprintf  (* 
cprintf (' 
cprintf  (’ 
cprintf  (' 
cprintf (' 
cprintf (' 
cprintf (' 
cprintf (' 


setcursortype(  NORMALCURSOR); 

the  longitudinal  flN  OF  TRIANGULAR  PROF ILE\r\n") ; 


=\r\n"); 


Fin  Depth  '•);  depth  *  i_L(L); 

Starting  Base  Fin  Width  ••);  start ing_width  =  i_L(d); 
Ending  Base  Fin  Width  ");  ending_wTdth  =  i_LTd); 

Profile  Area  ");  Ap  =  i_A(X); 

Emissivity  =  "j.-'scanf ("Xf",  iepsilon); 
Density  (  kg  /  dm*3  )  =  ");  scanfC'Xf",  &density); 
Base  Temperature  ");  Tbase  =  i_T(T); 

Thermal  Conductivity  ");  k  =  i_K(K); 
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cprintfC"  Environmental  Parameter  (0  to  .8)  =  •');  scanfC'Xt",  &V); 

if  (S  ==  1)  < 

cprintfC  Output  Filename  is  =  '•);  scanfCXs",  file_name);> 

cprintf ("\n  ===s======================================a===========s===========\r\n"); 


if  (  startin9_Hidth  >=  ending_Midth){ 

cprintr("\r\n\r\nFATAL  ERROR:  Starting  Base  Fin  Width  >=  Ending  Base  Fin  Width") 
getchO;  exit(0);> 

if  (  epsilon  >  1.0  |[  epsilon  <=  0)( 

cprintf("\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchO;  exit(0};> 

if  (  V  <  0  1 1  V  >  0.8){ 

cprintf ("\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8") 
getchO;  exittO);) 

_setcursortype(_NOCURSOR); 

return;) 


GLOBAL  FUNCTIONS 


•/ 


double  hr(void}< 

float  Tav  =  (Tb  +  Ta)  /  2; 

returnCsigma  •  epsilon  *  pow(Tav,  3)  -  (K2  /  (2  *  Tav)));> 
double  mCvoidX 

returnC  hr()  /  (k  *  sin(f in_taper)));) 

double  uCdouble  x){ 

returnC  2  •  sqrt<x  *  mO  )  );) 

double  A_(void){ 

return{  ”ua  *  ( 1 1  (ua)*K0(ub)  ♦  t0(ub)*IC1  (ua) )  );) 

double  B_(void)t 

return  (  ub*ub*(10(ub)*K0(ua)-10(ua)*K0(ub))  /  (  2*k*b_width*depth*mO));) 

double  C_(void){ 

return  (~2*lf*a_width*depth*mO*ub*(I1(ub)*K1(ua)  •  n(ua>*K1(ub))/(ua)  );) 
double  D_(void)( 

returnC  ub*(  l0(ua)*K1(ub)+I1(ub)*K0(ua))  );) 

y*  A************************************************************************** 


OPTIMIZE  Analysis  Program 
Longitudinal  Fin  of  Rectangular  Profile 

****************************************************************************  * ! 


^include 

<stdio.h> 

#include 

<stdl ib.h> 

#include 

<math.h> 

#include 

<conio.h> 

iKinclude 

"units.h" 

#def ine 

Yes  1 

#def ine 

NO  0 

#def ine 

N  100 

((define 

sigma  5.66961e-8 

double 

length,  width,  depth.  A,  8,  C 

100 


Tbase,  heat,  PN,  eta,  v,  K2,  epsilon,  k,  density,  mass,  Ap,  error; 
int  T,  L,  M,  S,  H,  d,  K,  i,  X,  2,  Y; 

void  inputO,  computeO,  resultsO; 

. . . 


MAIN  PROGRAM 


*/ 


int  inain(){ 

FILE  *data_in; 

data_in  =  fopen("units.ext'‘,‘'r"); 

fscanf(data_in,«%d  %d  %d  %d  Xd  5ld  %d  %d  %d  Xd  Xf,ST,&L.&M,&S,SH,&d,&K,&X,&Z,SY,aerror); 

fclose(data_in); 

error  =  error  /  100; 

inputO; 

heat  =  heat  /  (  1  -  V); 

Ap  =  (1.6178  *  pow(heat/depth  ,  3))  /  (  k  *  powCepsi lon,2)  •  pow(Tbase,9)  *  pow(sigma,2)) 
uidth  =  (1.8648*  pow(heat/depth  ,  2))  /  (  k  •  si9(na*epsi lon*pow(Tbase,5)); 

computeO; 
resultsO; 
exi t(0); 
return  0;> 

y*  ************»**»***««»***********»ik***»##*«iit**#**»***»**»*«**»**#*****dr 

COMPUTE  ROUTINE 


**•«**»***»***********«  *|»**#*  ****«(  ^  ^ 


void  compute(void){ 

double  A_0,  B_0,  C_0,  D_0,  mO,  hrO,  YbO; 

double  inc,  overall_diff  *  100,  Tmax  =  0,  Tbv_min  =  100000,  datalN+l] t6] ; 

double  A1,  B1,  Cl,  01,  Ttip,  Tave,  temp,  qa,  segment  diff; 
int  overflow,  possible_overf low  =  No,  number_of_tries=0; 

length  *  Ap  /  width; 

K2  =  2  *  epsilon  *  sigma  »  v  •  pow(Tbase,4); 

Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbaso  *  sqrt(  2  *  V  )  ♦  (  (  Tbase  -  (Tbase  •  sqrt(  2  *  V  >>)/2); 

else 

Ttip  =  .0001  +  Tave; 
inc  =  .5  *  Ttip;  i  =  1; 

while  (fabs(overall_diff)  >  (  Tbase  *  error  )  | |  i  <=  N  )( 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  •  width  *  epsilon  •  pow(Ta.  4;  •  depth; 

while  (i  <=  N  &&  overflow  ==  NoX 

segment_diff  =  10; 

while  (fabs(segment_diff)  >  (  Tb  •  error)  /  N  &S  overflow  ==  No)< 

A  =  A_();  B  =  B_();  C  =  C_();  0  -  0_(); 
temp  =A*Ta*B*  qa; 
qb  =C*Ta+D*  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 
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Tb  =  temp  ♦  segment_diff ; 

if  (Tb  >  Tbase  +50  1 j  Tb  <  0){ 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2  ; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  TaveT  Ttip  =  Tave;}> 

if  (overflow  ==  No){ 

data[i][0]  =  qb;  dataCildl  =  Tb;  datalil  [2]  =  A; 
dataCi]C31  =  B;  datalil  (4]  =  C;  datalillSl  =  D; 

Ta  =  Tb;  qa  =  qb;  ++i;)> 

if  (nutnber_of_tries  >  29){ 

clrscrO;  cprintf ("Fatal  Error;  Unable  to  Confute"); 
getchO;  exit(0);> 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  Yes)( 
inc  =  inc  /  2; 

if  (overall_diff  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  +  inc;) 

else 

Ttip  =  Ttip  -  overall_diff; 
number_of_tries++; 
gotoxy('27724); 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",number_of_tries,w_T(T,TTip),w_T(T,Tb) 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

for  (i  =  N  -  1;  i  >=  1;  i--)< 

A1  =  A  •  datalil  12)  +  B  •  datali]l41;  B1  =  A  •  datali)  13)  +  B  •  datalil  (5); 

0=0'  datali]  (21  +  0  •  datali)  14);  01  =  C  •  datali)  (3)  +  D  •  datali)  [51; 

A  =  A1;  B  =  B1;  C=C1;  0=01;) 

Ta  =  Ttip; 

eta  =  qb  /  (2  *  sigma  *  epsilon  *  length  *  depth  •  pow(Tbase,4)); 

PN  =  length  •  sqrt(sigma  *  epsilon  *  pow(Tbase,3)  /  (k  •  width)); 
mass  =  width  *  length  *  depth  •  density  *  1000; 
return; > 

*******  ***tk'*«*  ************  **************************************  ******* 

OUTPUT  PROCEDURE 

**************************************************************************  * ^ 


void  results(void)( 


clrscrO; 

cprintf("\a  OPTIMUM  PARAMETERS"); 

cprintf("\r\n  THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE"); 

cprintf  ("\r\n  r===============m=======i=========r=r========s====s=============\r\n"); 

cprintfl"  Density  =  %4.2f  kg  /  dm"3\r\n", density); 

cprintfl"  Fin  Depth  =  ");  p_L(L,  depth); 

cprintfl"  Total  Mass  =  ");  p_M(M,  mass); 

cprintfl"  Fin  Height  =  ");  p^LlL,  length); 

cprintfl"  Fin  Width  =  ");  p_L(d,  width); 

cprintfl"  Emissivity  =  X4.2f\r\n",  epsilon); 

cprintfl"  Profile  Area  =  ");  p_A()(,  Ap); 

cprintfl"  Profile  Number  =  %4.3f\r\n",PN); 

cprintfl"  Base  Heat  Flow  =  ");  p_H(H,  qb); 

cprintfl"  Fin  Efficiency  =  X4.3f\r\n",eta); 

cprintfl"  Base  Temperature  =  ");  p_T(T,Tb); 
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cprintfC  Edge  Temperature  =  ");  p_T(T,Ta); 

cprintfC  Thermal  Conductivity  =  p”lC(K,  k); 

cprintfC"  Environmental  Parameter  =  5t4.2f\r\n",  V>; 

cprintfC"  =========================s=======i===============================r\r\n"); 

cprintfC"  Transmission  Parameter  Matrix  \r\n\r\n">; 

cprintfC"  A  =  %4.4f  B  =  ",  A);  p  2CZ,B); 

cprintfC"  C  =  p  TCY.C); 

cprintfC"  0  =  %A.Af«,  D); 

getchC); 

return;) 


jtt  *********************************************************************** 

INPUT  PROCEDURE 

«*********************1k************»lk*** **********************<***«********  * j 


void  inputCvoidX 

clrscrC );  _setcursortypeC_NORMALCURSOR); 


cprintfC"\r\n  THE  LONGITUDINAL  FIN  OF  RECTANGULAR  PROFILE\r\n"); 

cprintf C"\r\n  ============:======:=============i==r==========r==i================\r\n") 

cprintfC"  Fin  Depth  depth  =  i_LCL); 

cprintfC"  Heat  Dissipation  “);  heat  =  i_HCH); 

cprintfC"  Emissivity  =  ");  scanfC"%f",  Sepsilon); 

cprintfC"  Density  C  kg  /  dm“3  )  =  ");  scanfC"%f",  Sdensity); 

cprintfC"  Base  Temperature  ");  Tbase  =  i_TCT); 

cprintfC"  Thermal  Conductivity  ");  k  =  i_KCIC); 

cprintfC"  Environmental  Parameter  CO  to  .8)  =  ");  scanfC"%f",  SV); 


if  C  epsilon  >  1.0  j|  epsilon  <=  0){ 

cprintfC"\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchC);  exitCO);) 

if  C  V  <  0  I  1  V  >  0.8){ 

cprintf C"\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchC);  exitCO);) 

_set  cursor typeC_N0CURS0R ) ; 
return;) 

j*  *«*««**»*w*********ir*»*ik«»«**»#*«*»«*«»*ik]k«ik**«*1k***«*«**«************« 

GLOBAL  FUNCTIONS 

**************************************************************************  * ! 


double  hrCvoidX 

float  Tav  =  CTb  ♦  Ta)  /  2; 

returnCsigma  *  epsilon  *  powCTav,  3)  -  CK2  /  C2  •  Tav)));) 
double  mCvoidX 

returnC  sqrtC  C  2.0  *  hrO)  /  C  k  *  width)  )  );) 
double  YbCvoid)( 

returnC  sqrtC  2.0  *  hrC)  *  k  *  width)  •  depth  );) 
double  A_CvoidX 

returnC  coshC  mC)  *  length  /  N  )  );) 

double  B_Cvoid){ 

return  C  sinhC  mC)  •  length  /  N  )  /  YbC)  );) 
double  C_CvoidX 

return  C  sinhC  mC )  *  length  /  N  )  *  YbC)  );) 
double  0_CvoidX 
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rvtunK  coth(  •()  *  Itngth  /  N  )  >;} 
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OPTIMIZE  Analysis  Program 
Longitudinal  Fin  of  Trapezoidal  Profile 

*********************«««**«***«««  ****************  1^ ! 

#include  <stdio.h> 

#include  <stdl ib.h> 

#include  <math.h> 

^include  <conio.h> 

#include  "units. h" 

#include  "bessel.h" 

#define  Yes  1 

#define  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

char  f i I e_name [201 ; 

double  length,  edge_width,  a_width,  b_Midth,  base_width,  depth,  inc_uidth; 

double  A,  8,  C,  0,  Ta,  a,  Tb,  qb=0,  ua,  ub,  fin_taper,  Constant_3; 

Tloat  112,  taper_ratio,  V,  eta,  PN,  epsilon,  heat,  k,  Tbase,  density,  mass,  Ap,  error; 

int  T,  L,  M,  S,  H,  d,  k,  i,  X,  Z,  Y; 

void  inputO,  computeC),  resultsO; 

FILE  *data_in,  *out; 

A******************************************************-**************** 

MAIN  PROGRAM 

**************************************************************************  k j 

int  main(K 

data_in  =  fopenCuni  ts.ext",  "r"); 

fscanf(data_in,‘'%d  %d  %d  %d  %d  %d  %d  %d  %d  %d  %f",&T,SL,£M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error  =  error  /  100; 

inputO; 

heat  =  heat  /  (1  •  V); 

Constant_3  =  0.3881*pow( taper_rat i o,2)  -  1 .0324*t3per_rat io  +  2.7522; 

Ap  =  ((0.5147  •  taper_ratio  +  1.1382)  *  pow(heat/depth  ,  3))  /  (  k  •  pow(epsi lon,2)  * 

pow(Tbase,9)  *  pow(sigma,2)); 

base_width  =  (Constant_3  *  pow(heat/depth  ,  2))  /  (  k  *  sigma*epsi lon»pow(Tbase,5)); 
edge_width  =  taper_ratio  *  base_width; 

conipute( }; 
resultsO; 
exit(O); 
return  0;> 

fk  kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkiikkkk-kkkkkkkkkkitkkkkkiikk 

COMPUTE  ROUTINE 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkiiiikk  ♦  f 

void  compute(void)( 

double  A_(),  B_(),  C_(),  0_(),  m(),  hrO,  u(double  x); 

double  length_inc,  b1,  db; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000,  data [N+1]  [6] ; 

double  A1,  B1,  Cl,  D1,  Ttip,  temp,  qa,  segtnent_dif f; 

int  overflow,  possible_overf low  =  No,  number_of_tr:es=0; 
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length  =  2  *  Ap  /  (  base_Midth  +  edge_Hidth  ); 
b1  =  length  /  (  1  ■  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2*b1  ); 
db  =  length  /  N; 

K2  =  2  •  epsilon  •  sigma  *  v  *  pow(Tbase,4); 

Tave  =  Tbase  *  powCV,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  •  Tbase  *  sqrt(  2  *  V  )  +  (  (  Tbase  -  (Tbase  *  sqrt(  2  *  V  )))/2); 

else 

Ttip  =  Tave  +  0.0001; 
inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overal l_di f f )  >  (  Tbase  *  error)  J |  i  <=  N  )( 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  edge_width  *  epsilon  *  pow(Ta,  A)  •  depth; 

a  =  b1  -  length; 

white  (i  <=  N  &&  overflow  ==  No)( 
segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  b1  -  length); 
b_width  =  edge_width  *  (  a  +  db  )  /  (  b1  -  length); 

while  (fabs(segment_diff )  >  (  Tb  *  error)  /  N  &S  overflow  ==  No){ 

ua  =  u(a)  ;  ub  =  u(a  +  db); 

A  =  A_();  B  =  B_();  C  =  C_();  0  =  0_(); 
temp  =A’Ta+B*  qa; 
qb  =  C  *  Ta  +  0  *  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  +  segment_di f f ; 

if  (Tb  >  Tbase  ■*  50  ]  |  Tb  <  0)( 

overflow  =  Tes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min”-  Tmax)  /”2; 

Ttip  =  fabs(Tbv_riin  •  inc); 
if  (Ttip  <  TaveT  Ttip  *  Tave;)> 

if  (overflow  ==  No)( 

datalillO)  =  qb;  dataliHD  =  Tb;  datati)  [23  =  A; 
data[i][3]  =  B;  data[i][A}  =  C;  datali)  [5]  =  D; 

Ta  =  Tb;  qa  =  qb;  a  =  a  +  db;  ++i;>> 


if  (number_of_tries  >  29 )( 

clrscrO;  cprintf ("Fatal  Error:  Unable  to  Compute"); 
getchO;  exit(O);) 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  *=  Yes)( 
inc  =  inc  /  2; 

if  (overal l_di f f  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  ♦  inc;> 

else 

Ttip  =  Ttip  -  overalt_diff; 
number_of_tries++; 
gotoxy(27,24); 

cprintfC'N  =  %d  Ta  =  %4.1f  Tb  =  %4.1f  ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 
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for  (i  =  N  -  1;  i  >=  1;  i--){ 

A1  =  A  •  data[i]C2]  ♦  B  *  dataCilCA];  B1  =  A  *  data[iH3]  ♦  B  *  data[iH5]; 

Cl  =  C  •  data[i][2]  +  D  *  data[i][4];  D1  =  C  *  data[i][3]  +  D  *  data[i]t5]; 

A  =  A1;  B  =  B1;  C=C1;  D=D1;) 

if  (S  ==  1){ 

out  =  fopenCf i le_name,  "w"); 
for  (i  =  1;  i  <=  N;  i++) 

fprintffout,  "%e,  %e,  %e\n",  w_L(d,  i  *  length  /  N),  w  TCT.dataCi]  [1] ),  w  H(H, 

datati]  [0])  ); 

fclose(out);> 

Ta  =  Ttip; 

Ap  =  .5  *  length  *  (  edge_width  +  base_width  ); 
mass  =  Ap  •  depth  *  density  *  1000; 

PN  =  2  •  sigma  *  epsilon  *  pow( length,  2)  *  poH(Tbase,3)  /  (k*  base_width); 
eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  •  depth  *  pow(Tbase,4)); 
return;) 

jit  ****************  *«***«1i********»'fr*********fk*A***4t««*4r******itr***ilr**ilr***itr 

OUTPUT  PROCEDURE 

**«**********************»*«***K***«**«*************«*********************  it ! 

void  resultsCvoidX 


cl rscr( ); 
cprintf ("\a\r\n 
cprintf  {"\r\n 
cprintf  (‘‘\r\n 
cprintf  (“ 
cprintf  (" 
cprintf  (•' 
cprintfC 
cprintf  (" 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC 
cprintfC  === 

cprintfC 
cprintfC  A  - 

cprintfC  C  = 

cprintfl"  D 

getchO; 
return; > 


OPTIMUM  PARAMETERS"); 

THE  LONGITUDINAL  FIN  OF  TRAPEZOIDAL  PROFILE"); 


Density 
F in  Depth 
Total  Mass 
Fin  Length 
Emissivity 
Prof i le  Area 
Prof i le  Number 
Base  Heat  Flow 
Fin  Efficiency 
F in  Edge  Width 
Fin  Base  Width 
Base  Temperature 
Edge  Temperature 
Thermal  Conductivity 
Environmental  Parameter 


.2f  kg  /  dm''3\r\n",densi ty); 
;  p_L(L,  depth); 

;  p~M(M,  mass); 

;  p_L(L,  length); 

.2f\r\n",  epsilon); 

;  p  A(X,Ap); 

.3fNr\n'',PN); 

;  p_H(H,  qb); 

.3f\r\n",eta); 

;  p_L(d,  edge_Hidth); 

;  p~L<d,  base  width); 

;  p_T(T,Tb);  ~ 

;  p_T(T,Ta); 

;  p_K(K,  k); 

.2f\r\n",  V); 


====\r\n"); 


Transmission  Parameter  Matrix  \r\n\r\n"); 
%4.4f  B  =  ",  A);  p_Z(Z,B); 

");  p_Y(Y,C); 

=  %4.4f",  D); 


=\r\n"); 


INPUT  PROCEDURE 


****************************************************************************  *  y 


void  inputlvoidX 

clrscrO;  _setcursortype<  NORMALCURSOR); 

cprintfC\r\n  THE'lONGI  TUDI NAL  FIN  OF  TRAPEZOIDAL  PROFILE\r\n"); 

cprintfC\r\n  ===========  =  =  :=============:======:=========================:===i====\r\n"); 

cprintfC  Fin  Depth  ");  depth  =  i_L(L); 

cprintfC  Heat  Dissipation  ");  heat  =  i_H(H); 

cprintfC  Taper  Ratio  =  ");  scanf C%f",&taper_ratio); 

cprintfC  Emissivity  =  ");  scanfCXf",  Sepsilon); 
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cprintfl"  Density  (  kg  /  dm‘J  )  =  “);  scanf(“Xf“,  (density); 

cprintfC  Base  Temperature  *•);  Tbase  *  i_T(T); 

cprintfC  Thermal  Conductivity  '•);  k  =  i_IC{K); 

cprintfC"  Environmental  Parameter  (0  to  .8)  •  '•);  scanfC'Xf",  iv); 

cprintf{"\n===========================================*============*==**=*»==Vr\n*'); 

if  (  epsilon  »  1.0  ||  epsilon  <=  OX 

cprintf ("\r\n\r\nFATAL  ERROR:  Emtssivity  is  not  between  0  and  1"); 
getchO;  exit(0);> 

if  (  V  <  0  ;  I  V  >  0.8X 

cprintf ("\r\n\r\nFATAL  ERROR:  Envi ronn«ntal  Parameter  is  not  between  0  and  0.8") 
getchO;  exit(0);> 

if  (  taper_ratio  >  .99  ||  taper_ratio  <=  OX 

cprintf (“\r\n\r\nFATAL  ERROR;  Taper  Ratio  is  not  between  0  and  1.0"); 
getchO;  exit(O);} 

_setcursortype(_NOCURSOR ) ; 
return;) 


GLOBAL  FUNCTIONS 


*/ 


double  hr(void){ 

float  Tav  =  (To  •*  Ta)  /  2; 

returnCsigma  *  epsilon  *  powtTav,  3)  -  (K2  /  (2  *  Tav)));) 
double  m(voidX 

returnC  hrO  /  (k  *  s  m(  f  i  n_taper ) ) );  ) 

double  uldouble  xX 

returnC  2  *  sqrtCx  *  m( )  )  );) 

double  A_(voidX 

returnC  ua  *  C  MCuaJ'kOCub)  »  10Cub)*KlCua))  );) 
double  B  CvoidX 

return  c“  ub»ub*C 10Cub)*K0Cua ) • lO(ua)*KO<ub) )  /  (  2*k*b_width*depth*mC ) ) ); ) 
double  C  CvoidX 

return  C~2*k*a_width»depth*mO*ub*C  1 1  Cub)*Kl  Cua)  -  1 1Cu3>*Kl  Cub)  )/Cua)  );) 
double  D_Cvoid)C 

returnC  ub*C  I0Cua)*K1 Cub)* 1 1 Cub)*K0Cua) )  );) 
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OPTIMIZE  Analysis  Program 
Longitudinal  Fin  of  Triangular  Profile 

. . . . .  ./ 

tinclude  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  <eonio.h> 

#include  "units. h" 

#include  "bessel.h" 

#define  Yes  1 

lYdefine  No  0 

#def i ne  N  99 

#define  sigma  5.66961e-8 

char  f i le_name[20] ; 

double  length,  edge_width,  base_width,  a_uidth,  b_uidth,  depth; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  ua,  ub,  fin_taper,  a; 

float  V,  K2,  epsilon,  k,  Tbase,  density,  mass,  Ap,  error,  eta.  PN,  heat; 

int  T,  L,  M,  S,  H,  d,  K,  i,  X,  Z,  Y; 

void  inputO,  resuitsi),  compute(); 

FILE  *data_in,  *out; 

/.  . . 

MAIN  PROGRAM 

*«*»*»****«»  *****«**#**»«*********«<»*»******«**«««***********«***«*******•  ig ! 

int  mainCX 

data  in  =  fopenCuni ts.ext","r"); 

fscanf (data_in,"%d  %d  %d  %d  %d  %d  %d  %d  %d  %f",iT,iL,iM, is, &H,&d,&<,&X,&Z,&Y, terror) 

fcloseldata~in); 

error  =  error  /  100; 

i nput ( ) ; 

heat  ~  heat  /  (1  •  V); 

Ap  »  (1,1234  •  pow(heat/depth  ,  3))  /  (  k  •  poH(sigma*epsi lon,2)  *  pow(Tbase,9)); 

base_width  =  (2.4852  *  pow(heat/depth  ,  2))  /  (  k  *  sigma*epsi lon*pow(Tbase,5)); 

computed; 
resultsO; 
ex  1 1 ( 0 ) ; 
return  0;) 

j*  A********************************************************************** 

COMPUTE  ROUTINE 

**************************************************************************  * ^ 

void  compute( )( 

double  A_(),  B_(),  C_(),  D_(),  m(),  hr(),  u(double  x); 

double  length_inc,  bT,  db;~ 

double  inc,  Tave,  overall_diff  =  100,  Tmax=0,  Tbv_min  =  100000,  data[N+1H6]; 

double  A1,  B1,  Cl,  01,  Ttip,  temp,  qa,  segment_dif f ; 

int  overflow,  possible_overf low  =  No,  number_of_tries  =  0; 


length  =  2  •  Ap  /  base_width; 
edge_uidth  =  base_width  /  100; 
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db  s  length  /  (  N  ♦  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*bl  ); 

K2  =  2  *  epsilon  *  sigma  *  V  *  pow(Tbase,4); 

Tave  =  Tbase  •  pow(V,  0.25)  ♦  0.0001; 

if  (V==0) 

Ttip  =  .8625  •  Tbase  •  sqrtl  2  •  V  )  ♦  (  (  Tbase  -  ITbase  •  sqrtC  2  *  V  )))/2); 

else 

Ttip  =  Tave  ♦  0.0001; 
inc  =  .5  •  Ttip;  i=1; 

while  (fabsloveral l_di f f )  >  (  Tbase  *  error)  |j  i  <=  N  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
segment_diff  =  10; 

while  (f abs(segment_di f f )  >  (  Tb  *  error)  /  N  S&  overflow  ==  No)( 

ub  =  u(db); 

temp  =  Ta  /  lO(ub); 

C  ~  2  *  k  •  edge_width  *  depth  *  pow(m(),2)  *  IKub)  /  (ub  •  lO(ub)); 
qb  =  C  *  temp; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  ♦  segment_di f f ; 

if  (Tb  >  Tbase  *50  11  Ib  <  0)( 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_mTn  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  TaveT  Ttip  =  Tave;)) 

if  (overflow  ==  No)( 

dataCiUOl  =  qb;  dataCiUH  =  Tb;  datalilUl  =  1; 
dataCi)C3]  *  0;  dataCiJ  [A)  =  C;  data[i)(5)  *  1; 

Ta  s  Tb;  qa  =  qb;  +♦:;) 

a  =  b1  -  length; 

while  (i  <=  N  S&  overflow  ==  No){ 
segmentdiff  =  10; 

a_width  =  edge_width  •  a  /  (  b1  -  length); 
b_width  =  edge_width  •  (  a  +  db  )  /  (  b1  -  length); 

while  (fabs(segment_diff)  >  (  Tb  •  error)  /  N  SS  overflow  ==  No){ 

ua  =  u(a)  ;  ub  =  u(a  ♦  db); 

A  =  A_();  B  =  B_();  C  =  C_();  D  =  D_(); 
temp  =  A  *  Ta  ♦  B  •  qa; 
qb  =C*Ta+0*  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  ♦  segment_diff ; 

if  (Tb  >  Tbase  +50  ![  Tb  <  0){ 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

if  (overflow  ==  No){ 

data(i][0)  =  qb;  datali}[1]  =  Tb;  data[i][2]  =  A; 
datalilCS)  =  B;  dataCiJCA)  =  C;  dataCiJ  [5]  =  0; 

Ta  =  Tb;  qa  =  qb;  a  =  a  +  db;  ♦♦i;» 
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if  (nutiiber_of_tries  >  29)( 

clrscrO;  cprintf ("Fatal  Error:  Unable  to  Compute"); 
getchO;  exit(0);> 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Tttp; 

overaU_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  Tes)C 
inc  =  inc  /  2; 

if  (overal l_diff  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  ♦  inc;> 

else 

Ttip  =  Ttip  -  overal l_diff; 
numbe^_of_tries++; 
gotoxy<27,24); 

cprintfC'N  =  %d  Ta  =  7.4. If  Tb  =  %4.1f  “,number_of_tries,w_T(T,Ttip),w_T(T,Tb)); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)> 

for  ( i  =  N  -  1 ;  i  >=  1 ;  i -  - )< 

A1  =  A  •  datati][2)  ♦  B  *  data[i]I4);  B1  =  A  •  data[i)t3)  +  B  *  data[i)[5]; 

Cl  =  C  *  data[i][2]  ♦  0  *  data[il[4);  D1  =  C  •  data[i]t3)  +  D  ♦  data[i)[5]; 

A  =  Al;  B  =  B1;  C=C1;  0=01;) 

length  =  length  +  db; 


if  (S  ==  DC 

out  =  fopen(f i le_name,  "w"); 
for  (i  =  1;  i  <=  N;  !•►+) 

fprintf(out,  "7.e,  7.e,  7.e\n",  w  L(d,  i  *  length  /  N),  w  T(T,data[i]  CD ),  w  H(H, 

dataCi]  [0]  )); 

fclose(out);) 

Ta  =  Ttip; 

mass  =  Ap  ♦  depth  •  density  ’  lOQQ; 

PN  =  2  *  sigma  *  epsilon  *  pow( length, 2)  *  pow(Tbase,3)  /  (k*base_width); 
eta  =  qb  /  (  2  *  sigma  *  epsilon  *  length  *  depth  *  pow(Tbase,4)); 
return;) 

*«***********»**********X****«*X*****-*«****H**WHX*****»«******»****»«** 

OUTPUT  PROCEOURE 

********x**********************************#******»rx**********************  * ! 


void  results(void)< 


clrscrO; 
cprintf ("\a 
cprintf ("\r\n 
cprintf ("\r\n 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintfC 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintfC 


OPTIMUM  PARAMETERS"); 

THE  LONGITUOINAL  FIN  OF  TRIANGULAR  PROFILE"); 


Density 
Fin  Depth 
Total  Mass 
F i n  Length 
Emissivi ty 
Prof i le  Area 
Prof i le  Number 
Base  Heat  Flow 
Fin  Efficiency 
Fin  Base  Width 
Base  Temperature 
Edge  Temperature 
Thermal  Conductivity 
Environmental  Parameter 


%4.2f  kg  /  dm''3\r\n",densi ty); 
");  p_L(L,  depth); 

");  p_M(M,  mass); 

");  p_L(L,  length); 

%4.2f\r\n",  epsilon); 

");  p_A()(,Ap); 

%4.3f\r\n",PN); 

»);  p  H(H,  qb); 
%4.3f\r\n",eta); 

");  p  L(d,  base  width); 

");  p~T(T,Tb);  “ 

");  p_T(T,Ta); 

■');  P_K(K,  k); 

%4.2f\r\n",  V); 


\r\n"); 


==\r\n"); 
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cprintfC  Transmission  Parameter  Matrix  \r\n\r\n"); 

cprintfC  A  =  XA.Af  B  =  A);  p  Z<Z,B); 

cprintfC  C  =  p_Y(Y,C); 

cprintfC  D  =  %4.4f‘',  D); 

getchO; 

return;) 

/.  . . . . . . 


INPUT  PROCEDURE 


*/ 


void  input(void}{ 

clrscrC ) ;_set cursor type(_NORMALCURSOR) ; 


cprintfC\r\n  THE  LONGITUDINAL  FIN  OF  TRIANGULAR  PROFILE"); 

cprintfC  Fin  Depth  “);  depth  =  i_L(L); 

cprintfC  Heat  Dissipation  ");  heat  =  i_H(H); 

cprintfC  Emissivity  =  scanfC%f",  &epsilon); 

cprintfC  Density  (  kg  /  dm^3  )  =  ");  scanfCXf",  &density>; 

cprintfC  Base  Temperature  ");  Tbase  =  j_T(T); 

cprintfC  Thermal  Conductivity  ");  k  =  i_K(K); 

cprintfC  Environmental  Parameter  (0  to  .8)  =  ");  scanfCXf",  &V); 

cprintf C\n==================================================================\r\n"); 

if  (  epsilon  >  1.0  ||  epsilon  <=  0)< 


cprintf C\r\n\r\nFATAL  ERROR:  Emissivity  is  not  between  0  and  1"); 
getchO;  exit(O);) 

if  (  V  <  0  1 1  V  >  0.8)( 

cprintf C\r\n\r\nFATAL  ERROR:  Environmental  Parameter  is  not  between  0  and  0.8"); 
getchO;  exit(O);) 

_setcursor type(_N0CURS0R ) ; 
return; > 

f*  *********************************************************************** 

GLOBAL  FUNCTIONS 

****«*************«i>»***«**»**i»*i» ******************************************  * ! 


double  hr(void){ 

float  Tav  =  (Tb  +  Ta)  /  2; 

returnCsigma  *  epsilon  *  pow(Tav,  3)  -  {K2  /  (2  *  Tav)));> 
double  m(void){ 

return!  hr()  /  (k  *  si n( f i ntaper ) ) ); } 

double  uldouble  x){ 

return!  2  *  sqrt!x  *  m!)  )  );> 

double  A_!void){ 

return!  ”ua  *  !11!ua)*K0!ub)  +  I0!ub)*l!l!ua))  );> 

double  B_!void)< 

return  !  ub*ub*! I0!ub)*k0!ua)- I0!ua)*K0!ub))  /  !  2*k*b_width*depth*m! ) ));) 
double  C_!void){ 

return  !  2*k*a_width*depth*m!  )*ub*!  1 1  !ub)*l!1!ua)  -  Il!ua)*K1!ub))/!ua)  );> 
double  0  !void){ 

return!  ub»!  I0!ua)»K1 !ub)+I 1 !ub)*K0!ua))  );> 
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/ 


UNIT  CONVERSION  Program  Header  File 


*/ 


#include  <stdio.h> 
•include  <conio.h> 
•include  <string.h> 
•include  <stdlib.h> 
•include  "units.h" 


float  i_T{int  T)< 
float  v7 
snitch  (T){ 

case  0: 
case  1: 
case  2: 
case  3: 
return(V);) 


cprintf ("(deg  K) 
cprintf("(deg  C) 
cprintf ("(deg  F) 
cprintf ("(deg  R) 


scanf("%f",  SV>;  break; 

scanfC'Xf",  SV);  V  =  V  ♦  273.15;  break; 

");  scanf("%f",  &V);  V  =  5  *  (  V  ♦  «59.4  )  /  9 
");  scanf("Xf".  SV);  V  =  5  *  V  /  9;) 


void  p_T(int  T, 
snitch  (T)< 

case  0: 
case  1: 
case  2: 
case  3: 
return; > 


float  V){ 

cprintf ("%4.2f  deg  K\r\n",  V);  break; 

V  =  V  -  273.15;  cprintf ("%4.2f  deg  C\r\n",  V);  break; 

V  =  9  *  (  V  -  255.222222)  /  5;  cprintf("%4.2f  deg  f\r\n",  V);  break; 

V  =  9  *  V  /  5;  cprintf(“%4.2f  deg  R\r\n",  V);) 


float  n_T( int  T, 
snitch  (T){ 

case  0: 
case  1 : 
case  2: 
case  3: 
return(V);> 


float  V){ 
break; 

V  =  V  -  273.15;  break; 

V  =  9  *  (  V  -  255.222222)  /  5;  break; 

V  =  9  •  V  /  5;) 


void  p_2(int  2,  float  V){ 
sni tch~(Z){ 

case  0:  cprintfC'XA.Af  deg  C  /  natt\r\n",  V);  break; 

case  1:  V  =  0.1760831  *  V  ;  cprintf ("%4. 34  deg  F  hr  /  BTU\r\n",  V);) 

return;) 


float  n_Z(int  Z,  float  V){ 
snitch  TZ){ 

case  0:  break; 

case  1:  V  =  0.1760831 

return(V);) 


V;) 


void  p_Y(int  Y,  float  V)( 
snitch  (Y)( 

case  0:  cprintf ("’44. 4f  natt  /  deg  C",  V);  break; 

case  1:  V  =  V  /  0.1760831  ;  cprintf("%4.4f  btu  /  hr  deg  F",  V);) 

return; > 


float  n_Y(int 

Y,  float  V){ 

snitch  (Y){ 

case  0; 

break; 

case  1 : 

V  =  V  / 

0.1760831;) 

return(V);) 

float  i_L(int 

L){ 

float  V; 

switch  (L){ 

case  0: 

cprintf ( 

"  (m) 

S 

case  1: 

cprintf ( 

"(cm) 

= 

case  2: 

cprintf ( 

"(ft) 

= 

case  3: 

cprintf ( 

"(in) 

= 

case  4; 

cprintf ( 

"(mm) 

= 

scanfC'Xf",  &V);  break; 
scanfC'Xf",  SV);  V  =  V  /  100;  break; 

scanfC'Xf",  SV);  V  =  0.3048  *  V;  break; 

scanfC'Xf",  SV);  V  =  0.0254  *  V;  break; 

scanfC'Xf",  SV);  V  =  V  /  1000;) 


;  break; 


return(V);> 


void  p_L(int  L,  float  V){ 
switch  (DC 

case  0:  cprin 

case  1 :  V  =  1 

case  2:  V  =  V 

case  3:  V  =  V 

case  4:  V  =  1 

return;} 


cprintf ("%A .2f  m\r\n",  V);  break; 

V  =  100  •  V;  cprintf (‘'%4.2f  cm\r\n",  V);  break; 

V  =  V  /  0.3048;  cprintf ("%4.2f  ft\r\n“,  V);  break; 

V  =  V  /  0.0254;  cprintf ("*44. 2f  in\r\n«,  V);  break; 

V  =  1000  •  V;  cprintf("%4.2f  (nm\r\n",  V);> 


float  w_L(int  L,  float  V)C 
switch  (DC 


case 

0 

break; 

case 

1 

V  =  100  *  V;  break 

; 

case 

2 

V  =  V  /  0.3048; 

break, 

case 

3 

V  =  V  /  0.0254; 

break 

case 

4 

V  =  1000  *  V;) 

return(V);) 

float 

i  A(int 

A){ 

float 

V; 

switch 

(A)( 

case 

0 

cprintfC"  (m"2) 

= 

"); 

scanfC'Zf", 

case 

1 

cprintf ("(cm"2) 

"); 

scanfC'Zf", 

case 

2 

cprintf("(ft*2) 

= 

“); 

scanfC'Zf", 

case 

3 

cprintf("(in*2) 

= 

scanfC'Zf", 

case 

4 

cprintf  ("(tm'2) 

= 

scanfC'Zf", 

SV);  V  =  V  /  10.76;  break; 


return(V);) 

void  p_A(int  A,  float  V)C 
switch  (A)C 

case  0;  cprin 

case  1 :  V  =  1 

case  2:  V  =  V 

case  3:  V  =  V 

case  4:  V  =  1 

return;) 


cprintf ("%4. 2  m  '  2\r\n",  V);  break; 

V  =  10000  •  V;  cprintf(''%4.2f  cm  *  2\r\n",  V);  break; 

V  =  V  •  10.76;  cprintf("%4.2f  ft  '  2\r\n",  V);  break; 

V  =  V  »  1550;  cprintf (•'%4.2f  in  *  2\r\n’’,  V);  break; 

V  =  1000000  *  V;  cprintf("%4.2f  mm  *  2\r\n",  V);) 


float  w_A(int  A,  float  V){ 
switch  (A)C 

case  0:  break; 

case  1:  V  =  10 

case  2;  V  =  V 

case  3:  V  =  V 

case  4:  V  =  10 

return(V);> 


break; 

V  =  10000  *  V;  break; 

V  =  V  •  10.76;  break; 

V  =  V  *  1550;  break; 

V  =  1000000  »  V;) 


float  i_M( int  H)C 
float  V; 
switch  (M)C 

case  0: 
case  1 : 
case  2: 
case  3; 
returnCV);) 


cprintfC  (kg)  =  “);  scanfC’Xf",  &V);  break; 

cprintfC  (gm)  =  ");  scanf(“%f“,  8V);  V  =  V  /  1000;  break; 

cprintfC'dbm)  =  ");  scanfC'Zf",  &V);  V  =  V  /  2.204622;  break; 

cprintfC"  (oz)  =  ");  scanfC'Zf",  &V>;  V  =  V  /  35.27396;  break;) 


void  p_M(int  M,  float  V)C 
switch  (M)< 

case  0:  cprin 

case  1 :  V  -  1 

case  2:  V  ^  2 

case  3:  V  =  3 

return;) 


cprintf ("Z4.2f  kg\r\n",  V);  break; 

V  -  1000  *  V;  cprintf("Z4.2f  grams\r\n",  V);  break; 

V  s  2.204622  *  V;  cprintf ("Z4.2f  lbm\r\n",  V);  break; 

V  =  35.27396  *  V;  cprintf ("Z4.2f  oz\r\n",  V);) 


float  w_M(int  M,  float  V)( 
switch  {M){ 

case  0: 
case  1 : 


break; 

V  =  1000  '■  V;  break; 
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case  2:  V  =  2.20A622  •  V;  break; 

case  3:  V  =  35.27396  •  V;) 

return(V);) 


float  H)< 

float  vj 
switch  (H){ 

case  0: 
,:ase  1 : 
case  2: 
case  3: 
return(V);) 


cprintf (“(watts)  =  “);  scanf(”)lf",  iV);  break; 
cprintfC  (kw)  =  “);  scanf("%f".  SV);  V  =  1000  •  V;  break; 
cprintf("(BTU/sec)  =  ");  scanf ("Xf",  iV);  V  =  V  •  1054.3502;  break; 
cprintfC  (BTU/hr)  =  “);  scanfCXf",  SV);  V  =  V  /  3. 414425;) 


void  p_H(int  H, 
switch  (H){ 

case  0: 
case  1 : 
case  2: 
case  3: 
return;) 


float  V)( 

cprintf ("%4.2f  watts\r\n",  V):  break; 

V  =  V  /  1000;  cprintf CX4.2f  kw\r\n",  V);  break; 

V  =  V  /  1054.3502;  cprintf CX4.2f  BIU  /  sec\r\n",  V);break; 

V  =  3.414425  *  V;  cprintf CX4.2f  BIU  /  hr\r\n",  V);) 


float  w_H(int  H, 
switch  (H)< 

case  0: 
case  1 : 
case  2: 
case  3: 
return(V);) 


float  V){ 
break; 

V  =  V  /  1000;  break; 

V  =  V  /  1054.3502;  break; 

V  =  3.414425  *  V;) 


float  i_K(int  K){ 
float  V; 
switch  (!(){ 

case  0;  cprintf ("(watts  /  m  deg  K)  *  “);  scanf("%f",  SV);  break; 

case  1:  cprintf ("(BTU  /  ft  hr  deg  F)  »  '•);  scanfCXf",  SV);  V  =  V  /  1.731;) 

return(V);) 


void  p_K(int  K,  float  V){ 
swi tch~(K)( 

case  0:  cprintf CX4.2f  watts  /  m  deg  l(\r\n",  V);  break; 

case  1;  V  =  1.731  *  V;  cprintf CX4.2f  BTU  /  ft  hr  deg  F\r\n",  V);) 

return; ) 
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y*  ************************************************************************** 


BESSEL  FUKCTION  Program  Header  File’® 

y*  *************************************************************************** 

#inctude  <stdio.h> 

^include  <math.h> 

#include  "bessel.h" 

double  KOCfloat  x}{ 

double  y,  z; 

double  P1=-0. 57721566,  P2=0. 4227842,  P3=0. 23069756,  P4=0. 0348859,  P5 =0.00262698, 

P6=0. 0001075,  P7=0. 0000074; 

double  01=1.25331414,  Q2=-0. 07832358,  03=0.02189568,  04=-0. 01062446,  05=0.00587872, 

06= -0.0025 1546; 
double  07=0.00053208; 

if  (  X  <  2  )< 

y=(x*x)/4; 

z  =  (-log(x/2j*l0(x))+(P1+y*(P2*y»<P3+y*(P4+y*(P5*y*(P6-ry*P7))  ))>>;> 

else{ 

y  =  2  /  x; 

z  =  (exp(  -x)/sqrt  (x)  )"(Q1+y*<Q2+y*(Q3'»y’(Q4'*y*(Q5»y*(06»y*Q7) )  ))));> 
returnlz);) 

double  KKfloat  x)L 

double  y,  z; 

double  P1=1.0,  P2=0. 15443144,  P3=- .67278579,  P4=-0. 18156897,  P5=-0. 01919402,  P6=-0. 001 10404, 

P7=- 0.00004686; 

double  01=1.25331414,  02=0.23498619,  03= -0.03655620,  04=0.01504268,  05= -0.00780353, 

06=0.00325614; 
double  07=-0.00068245; 

if  (  X  <  2  ){ 

y  =  (  X  *  X  )  /  4; 

z  =  <log(x/2)*I1(x))-<-(1/x)»(Pl*y*(P2*y»(P3+y»(P4*y»(P5+y»(P6+y*P7))))));} 

else{ 

y  =  2  /  x; 

Z  =  (exp(-x)/sqrt(x))*(0l+y*(Q2fy*(03»y*(04+y»(05+y*(06+y*07))))));> 
return(z);> 


double  lOlfloat  x){ 
double  y,  z,  ax; 

double  P1=1,  P2=3. 5156229,  P3=3. 0899424,  P4=1 .2067492,  P5=0. 2659732,  P6=0. 0360768, 

P7=0. 0045813; 

double  01=0.39894228,  02=0.01328592,  03=0.00225319,  04=-0. 00157565,  05=0.00916281, 

O6=-0. 0205 7706; 

double  07=0.02635537,  08=-0. 01647633,  09=0.00392377; 

if  (fabsCx)  <  3.75  ){ 

y  =  (  X  /  3.75  )  •  (  X  /  3.75  ); 
z  =  Pl+y*(P2+y*(P3-ry*(P4+y*(P5+y*(P6+y*P7)))));> 

else{ 

ax  =  fabsCx); 
y  =  3.75  /  ax; 

z  =  (exp(ax)/sqrt(ax))*l01+y*<Q2+y*(03+y*(04+y*(05»y*(06*y*(07+y*(08+y*09) ) ))))));> 
returnlz); > 

double  IKfloat  x)< 

double  y,  z,  ax; 

double  P1=0.5,  P2=0. 87890594,  P3=0. 51498869,  P4=0. 15084934,  P5=0. 02658733,  P6=0. 00301532, 

P7=0. 000324 11; 

double  01=0.39894228,  02= -0.03988024,  O3=-0. 00362018,  04=0.00163801,  05=-0. 01031555, 


116 


06=0.02282967; 

double  Q7=-0. 02895312,  08=0.01787654,  Q9=-0. 00420059; 

if  (fabs(x)  <  3.75  ){ 

y  =  (  X  /  3.75  )  *  (  X  /  3.75  ); 

2  =  x»(P1+y»(P2»y*(P3+y»{P4+y*(P5+y»(P6»y*P7))))));) 

else{ 

ax  =  fabs(x); 
y  =  3.75  /  ax; 

2  =  <exp(ax)/sqrtCax))*(01+y»(02+y*(Q3+y*(04+y*(05+y*{06+y*(Q7+y*(08*y*09))))))));> 
return(2);> 
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APPENDIX  B  COMPUTER  CODE  TO  DERIVE  OPTIMUM  GEOMETRY 


*************************************************************************** 

profile  area  Analysis  Program  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Rectangular  Profile 

*********************************  AM*****************************************  * ! 

^include  <time.h> 
tHinclude  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  <conio.h> 

#define  Yes  1 

#def i ne  No  0 

#define  N  100 

#define  sigma  5.66961e-8 

double  length,  width_base,  raa;;_len9th; 

double  A,  B,  C,  D,  Ta,  Tb,  qb=0,  Heat,  qb  max,  db.  Continue; 

float  V,  tC2,  epsilon,  k,  Tbase,  Ap,  taper_ratio,  depth,  erronO.Ol; 

int  i,  j,  number_of _p3sses ,  failures; 

void  ComputeO; 

. . . . . 

MAIN  PROGRAM 

**********************.*.*..**......*...**.. . . . . . *,..*  */ 

int  mainOf 
float  dq; 

FILE  ‘out; 
randomizef); 

-  error  /  100; 
cliv,cr(  ); 

I  or  (j  =  1;  j<102; 

nuniLicr_of_passes  =  0; 

epsilon  =  (float)  <  random(98)  ); 

epsilon  =  .01  +  epsilon  /  100; 

depth  =  (float)  (20  +  random(200)); 

depth  =  depth  /  100; 

k  =  10  ♦  random(200); 

Tbase  =  100  ♦  random(500); 

Heat  =  (float)  (25  ♦  random( 1000)); 

Heat  =  Heat  /  depth; 

V  =  (float)  (  random(80)); 

V  =  .01  +  V  /  100; 

V  =  0; 

taper_ratio  =  (float)  (random(98)); 
taper_ratio  =  0.01  ♦  taperratio  /  100; 
taper_ratio  =  1; 

Continue  =  Yes; 
qb_max  =  - 1 ; 
qb  =  0; 
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dq  =  0; 

Ap  =  1.635  *  pow(Heat,3)  /  (  k  *  pow(sigma,2)  *  powCepsi lon,2)  *  pow(Tbase,9)); 

length  =  .87  •  Heat  /  (  sigma  •  epsilon  »  pOM(Tbase,4) ); 

db  =  length  /  20016; 

length  =  length  -  db; 

width_base  =  Ap  /  length; 

if  (depth  <  10*uidth_base)  Continue  =  No; 
if  (length  <  10*width_base)  Continue  =  No.- 

while  (Continue  ==  YesX 

length  =  length  +  db; 
width_base  =  Ap  /  length; 
nuniber_of_passes++; 
failures  =  0; 

computeO; 

while  (qb  ==  0  &&  Continue  ==  Yes)( 
fai  lures+-*'; 

length  =  length  +  db  /  2; 
width_base  =  Ap  /  length; 
computeO; 
if  (fai lures  ==  3)< 

Continue  =  No; 
nuniber_of_passes  =  1;)) 


dq  =  qb  -  dq; 

cprintf (“\r\nwidth_base  =  /i4.3f  length  =  XA.bf  heat  =  %4.6f  dq  =  %4.6f", 
width_base*1000,  length,  qb,  dq); 
dq  =  qb; 


if  (qb_max  <  qb  )( 

max_length  =  length; 
qb_max  =  qb;> 

else 

Continue  =  No;) 

if  (number_of_passes  >  2){ 
length  =  max_length; 
width_base  =  Ap  /  length; 
computeO; 

cprintf ("\r\nWri t ing  Rec  #  %d  d  =  %4.3f  b  =  %4.6f  qb  =  %4.3f  Tb  =  %4.1f  Ap  *  514. If", 
j,  width_base,  length,  qb,  Tbase,  Ap); 
out  =  fopen("OPT_RECT.MAT","a"); 

fprintfCout,  "7.d,  %e,  7.e,  %e,  %e,  %e,  %e,  7.e,  %e,  %e,  Xe,  Xe\r\n",  j, 

depth,  length,  Ap,  k,  Tbase,  epsilon,  V,  qb,  taper_ratio,  width_base,  width_base  ); 
fclose(out);) 

else( 

cprintf ("\r\nIgnoring  This  Set  Of  Parameters"); 

return  0; > 

***********************«*******«**********'*****«**««*«**««**««***1l«*«X* 

COMPUTE  ROUTINE 

*** »**********************************************ix***********************  * f 


void  co(npute(void){ 

double  Tav,  inc,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  =  100000; 

double  A,  B,  C,  D,  m,  hr.  To,  Ttip,  temp,  qa,  Tave,  segment_di f f ; 

int  oversow,  possible_overf low  =  No,  number_of_tries=0;” 

width_base  =  Ap  /  length; 
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K2  =  2  *  epsilon  •  sigma  *  V  *  pow(Tbase,4); 

Tave  =  Tbase  *  pow(V,  0.25)  +  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  •  V  )  ♦  (  (  Tbase  -  (Tbase  •  sqrt(  2  *  V  ))>/2); 
else  • 

Ttip  =  .0001  +  Tave; 

inc  =  .5  *  Ttip;  i  =  1; 

uhile  (tabs(overal  l_di  f  f )  >  {  Tbase  *  error)  j i  <=  N  ){  ' 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  *  width_base  *  epsilon  *  powCTa,  4)  *  depth; 

while  (i  <=  N  &&  overflow  ==  No){ 

segment_diff  =  10; 

while  (fabs(segment_diff )  >  (  Tb  *  error)  /  N  8.&  overflow  ==  No)< 

Tav  =  (Tb  +  Ta)  /  2; 

hr  =  sigma  *  epsilon  •  pow(Tav,  3)  -  (K2  /  (2  *  Tav)); 
m  =  -qrt(  (  2.0  *  hr)  /  (  k  *  width_base)  ); 

Yo  -  sqrtl  2.0  *  hr  *  k  *  width_base)  *  depth; 

A  =  cosh<  m  *  length  /  N  ); 

B  =  sinh{  m  *  length  /  N  )  /  Yo; 

C  =  sinh(  m  *  length  /  N  )  *  Yo; 

D  =  A; 

temp  ::A*Ta  +  B*  qa; 
qb  =C*Ta+D*qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  »  segment_di f f ; 

if  (Tb  >  Tbase  *50  ||  Tb  <  0)( 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_mTn  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;>> 

if  (overflow  ==  No)( 

Ta  =  Tb;  qa  =  qb;  ++i;)> 

if  (nurober_of_tries  >  29)( 

Ta=0;Tb=0;qb=0; return; ) 

if  (overflow  ==  No)( 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase  )  Tmax  =  Ttip; 

overall_diff  =  (Tb  •  Tbase)  /  2; 

if  (  possible_overf low  ==  Yes){ 
inc  =  inc  /  2; 

if  (overal l_di f f  >0)  Ttip=Ttip-  inc; 
else  Ttip  =  Ttip  +  inc ;> 

else 

Ttip  =  Ttip  -  overal l_di ff ; 
number_of _t  r i es++ ; 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

return; ) 
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I*  ***************  ************************************************************* 

PROFILE  AREA  Analysis  Program  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Trapezoidal  Profile 

*************  **lk4r***)lr***l****iktt  ***********************************  **itritr**ilritr***  *^ 

^include  <time.h> 

#include  <stdio.h> 

^include  <stdlib.h> 

# include  <math.h> 

#include  <conio.h> 

Oinclude  "bessel.h" 

#define  Yes  1 

#define  No  0 

(Kdefine  N  100 

^define  sigma  5.66961e-8 

double  max_length,  length,  depth,  edge_width,  base_width,  Ta,  Tb,  qb=0,  inc,  Tave; 

float  dq,  V,  qb_max,  db,  epsilon,  taper_ratio,  k.,  error=0.01,  Tbase,  K2,  Heat,  Ap; 

int  i.  Continue,  number_of_passes,  j.  Exit,  failures; 

void  computeO; 

j-k  **********«******************«!»**********«*«**«****«******************* 

MAIN  PROGRAM 

***i»**«****»***»**W»**»****»**»**«#»***««K**ik**tt*«*«****it*ik**»*ik*««***lt***  * ! 

int  mainOf 

FILE  *out; 
random! ze( ); 
error  =  error  /  100; 

ClrscrO; 

for  (j=1;  j  <  102;  j»»)( 

number_of_passes  =  0; 
epsilon  =~(float)  <  randomC98)  ); 
epsilon  =  .01  »  epsilon  /  100; 
depth  =  (float)  (20  t  random(200)); 
depth  =  depth  /  100; 
k  =  10  +  random(200); 

Tbase  =  100  f  random(500); 

Heat  =  (float)  (25  +  random( 1000) ); 

Heat  =  Heat  /  depth; 

V  =  (float)  (  random(80)); 

V  =  .01  ♦  V  /  100; 

V  =  0; 

taper_ratio  =  (float)  (random(98)); 
taper_ratio  =  0.01  ♦  taper_ratio  /  100; 

Continue  =  Yes; 

Exit  =  No; 
qb_tnax  =  -1; 
dq  =  0; 
qb  =  0; 

Ap  =  (.63077  *  taper_ratio  +  1.31723)  *  pow(Heat,3)  /  (  k  *  pow(sigma,H)  *  powCepsi lon,2) 
pow( Tbase, 9) ) ; 

length  =  .82  *  (-.0315113  *  taper_ratio  +  1.0A25113)  •  Heat  /  (  sigma  *  epsilon  • 
pow(Tbase,4)); 
db  =  length  /  20016; 

length  =  length  -  db; 

base_width  =  2  *  Ap  /  (  length  *  (  taper_ratio  ♦  1)); 
edge_width  =  taper_ratio  *  base_width; 


if  (depth  <  10*base_Hidth )  Continue  =  No; 
if  (length  <  10*base_nidth)  Continue  =  No; 

while  (Continue  ==  Yes)( 

length  =  length  ♦  db; 

base_width  =  2  *  Ap  /  (  length  *  (  taper_ratio  »  1)); 
edge_width  =  taper_ratio  *  base_width; 

failures  =  0; 
nuniber_of_passes*+; 

computet ); 

while  (qb  ==  0  &&  Continue  ==  Yes)( 
fai lures**; 

length  =  length  *  db  /  2; 

base_width  =  2  *  Ap  /  (  length  •  (  taper_ratio  *  1)) 
edge_width  =  taper_ratio  *  base_width; 
computet ); 
if  (failures  ==  3)( 

Continue  =  No; 
nutnber_of_passes  =  1;)> 


dq  =  qb  -  dq; 

cprintf ("\r\nwidth  =  %^.3(  length  =  XA.6f  heat  =  dq 

length,  qb,  dq); 
dq  =  qb; 


if  (qb_max  <  qb  )( 

max_length  =  length; 
qb_max  =  qb; ) 

else 

Continue  =  No;  ) 

if  (number_of_passes  >  2)( 
length  =  max_length; 

base_width  -~2  *  Ap  /  (  length  *  (  taper_ratio  *  1)); 
edge_Hidth  =  taper_ratio  »  base_width; 
computeO; 

cpr  int f  ( "\r\nwr i  t  i ng  Rec  #  %d  d  =  %4.3f  b  =  %4.6f  qb  = 
j,  base_width,  length,  qb,  Tbase,  Ap); 
out  =  fopen("OPTrRAP.MAT’‘,"a"); 

fprintftout,  "%d,  %e,  %e,  %e,  Se,  %e,  %e,  %e,  Xe, 

depth,  length,  Ap,  k,  Tbase,  epsilon,  V,  qb,  taper_rati( 
f closetout ) ; ) 

else{ 

cprintf ("\r\nlgnoring  This  Set  Of  Parameters"); 

return  0;> 

void  compute(void){ 

double  m,  hr,  n,  ua,  ub.  A,  B,  C,  0,  fin_taper,  I1UA,  KOUB,  lOUB, 

double  qa,  b1 ,  dH,  a,  a_width,  b_width,  Ttip,  temp,  segment_di f f , 

double  b,  inc,  Tav,  Tave,  overall_diff  =  100,  Tmax  =  0,  Tbv_min  = 

int  overflow,  possible_overf low  =  No,  number_of_tries  =  0; 

b1  =  length  /  (  1  -  taper_ratio  ); 
fin_taper  =  atan2(  base_width,  2  *  b1  ); 
dH  =  length  /  N; 

K2  =  2  •  epsilon  *  sigma  *  V  *  pow(Tbase,A); 

Tave  =  Tbase  *  pow(V,  0.25)  *  0.0001; 

if  (V==0) 

Ttip  =  .8625  *  Tbase  *  sqrt(  2  *  V  )  *  (  (  Tbase  -  (Tbase  *  sqrtt 

else 


=  XA.bf",  base_width*1000, 


3f  Tb  =  XA.If  Ap  =  X4.1f", 


Xe,  Xe,  Xe\r\n",  j, 

>,  base_width,  edge_width  ) 


KlUA,  KOUA,  lOUA,  KlUB; 
HUB; 

100000; 


2  •  V  )))/2); 
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Ttip  =  Tave  »  0.0001; 
inc  =  .5  *  Ttip;  i=1; 

while  (fabsCoveral  l_diff )  >  (  Tbase  *  error)  j|  i  <=  *|  ){ 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 

qa  =  sigma  •  edge_width  •  epsilon  •  pow(Ta,  4)  *  depth; 

a  =  bl  -  length; 

while  (i  <=  N  &S  overflow  =-  No}{ 
segment_diff  =  10; 

a_width  =  edge_width  *  a  /  (  bl  -  length); 
b_width  =  edge_width  •  {  a  ♦  dH  )  /  (  bl  -  length); 

while  (fabs(segment_di f f )  >  <  Tb  *  error)  /  N  overflow  ==  NoX 

Tav  =  (Tb  ♦  Ta)  /  2; 

hr  =  sigma  *  epsilon  •  pow(Tav,  3)  -  (K2  /  (2  •  Tav)); 
n  =  hr  /  (k  •  sin(f in_taper)); 
ua  =  2  *  sqrt(a  •  n  ); 
ub  =  2  *  sqrt(  (a  ♦  dH)  *  n  ); 

IlUA  =  IKua); 

KOUB  =  KOlub); 
lOUB  :  I0<ub); 
klUA  =  K1(ua); 

KOUA  =  k0{ua); 
lOUA  =  lO(ua); 

KlUB  =  KKub); 

HUB  =  l1<uD); 

A  =  ua  *  (HUA'KOUB  ♦  I0UB*KlUA); 

B  =  ub*ub*(lOUB»KOUA  ■  I0UA*K0UB)  /  (  2*k*b  width*depth*n); 
C  =  2*k*a  w)dth*depth*n*ub»( I lUB'KIUA  -  nUA*K1UB)/(ua); 

0  =  ub*(  T0UA*K1UB  +  nUB*K0UA); 
temp  =  A  *  Ta  ♦  B  •  qa; 
qb  =  C  *  Ta  +  D  •  qa; 
segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  ♦  segment_diff ; 

if  (Tb  >  Tbase  ♦  50  ; I  Tb  <  0){ 

overflow  =  Yes;  possible_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  •  inc); 
if  (Ttip  <  Tave)  Ttip  =  Tave;)> 

if  (overflow  ==  No)( 

Ta  =  Tb;  qa  =  qb;  a  =  a  ♦  dH;  »+i;)> 

if  (number_of_t  r  I  es  :>  29){ 

Ta=0; Tb=0;qb=0; return; ) 

if  (overflow  ==  No){ 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  Yes){ 
inc  =  inc  /  2; 

if  (overal l_di f f  >  0  )  Ttip  =  Ttip  -  inc; 
else  Ttip  =  Ttip  +  inc;> 

else 

Ttip  =  Ttip  •  overal l_di ff; 
number_of_t  r i es++ ; 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

return;) 
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/* 


PROFILE  AREA  Analysis  Program  Tb  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Triangular  Profile 

****************************************************************************  * y 


#include  <time.h> 
#include  <stdio.h> 
^include  <stdl ib.h> 
#include  <math.h> 
iVinclude  <conio.h> 
#include  “bessel.h" 


#def ine 
#def i ne 
#def i ne 
#def ine 


Yes  1 

No  0 

N  99 

sigma  5.66961e*8 


double 
doub I e 
float 
int 
void 


length,  depth,  edge  width,  base_width,  max  length; 

Ta,  Tb,  qb=0; 

V,  qb_max,  db,  epsilon,  k,  error=0.01,  taper_ratio,  Tbase,  K2,  Heat,  Ap; 

i.  Exit,  Continue,  numfcer_of_passes,  j,  failures; 

computeO; 


<r4t*«**»»*****ir***«*»*tt»*****«K«K*«**«KK'<t«***«***«*««*r'«*<tli***li'k«***«‘ft««‘ft 


MAIN  PROGRAM 

****«***«)t*ik*»»»««r**»*ir»»»KK»»«**«*K»»*it*iir««***K**tt***«ttir*******«**»*«**ik*  *y 

int  main(){ 
float  dq; 

FILE  *out; 
randomizel ); 
error  =  error  /  100; 
clrscrO; 

for  (j=1;  j  <  102;  ]♦♦)< 

number_of_passes  =  0; 

epsilon  =  (float)  (  random(98)  ); 

epsilon  =  .01  +  epsilon  /  100; 

depth  =  (float)  (20  +  random(200)); 

depth  =  depth  /  100; 

k  =  10  +  random(200); 

Tbase  =  100  ♦  random(500); 

Heat  =  (float)  (25  +  random( 1000)); 

Heat  =  Heat  /  depth; 

V  =  (float)  (  random(80)); 

V  =  .01  +  V  /  100; 

V  =  0; 

taper_ratio  =  (float)  (random(98)); 
taper^ratio  =  0.01  +  taper_ratio  /  100; 
taper_ratio  =  0; 

Exit  =  No; 

Continue  =  Yes; 
qb_max  =  - 1 ; 
qb”=  0; 
dq  =  0; 

Ap  =  1.304  •  pow(Heat,3)  /  (  k  •  pow(sigma,2)  *  pow(epsi lon,2}  •  pow(Tbase,9}); 

length  =  .95  *  Heat  /  (  sigma  •  epsilon  *  pow(Tbase,4)}; 

db  =  length  /  20016; 

length  =  length  -  db; 

base_width  =  2  *  Ap  /  length; 
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if  (depth  <  10*base_width)  Continue  =  No; 
if  (length  <  10*base_width)  Continue  =  No.- 

while  (Continue  ==  YesX 

length  =  length  ♦  db; 
base_width  =  2  *  Ap  /  length; 

nLmber_of_passes+*; 
fai lures  =  0; 

coinputeO; 

while  (qb  ==  0  &&  Continue  ==  Yes){ 
fai lures++; 

length  =  length  +  db  /  2; 
base_width  =  2  *  Ap  /  length; 
conputel  ); 
if  (fai lures  ==  3){ 

Continue  =  No; 
nufnber_of_passes  =  1;>> 


dq  =  qb  -  dq; 

cprintf("\r\ndelta  =  %A.3f  length  =  %A.6f  heat  =  W.6f  dq  =  X4.6f‘',  base_width*1000, 
length,  qb,  dq); 
dq  =  qb; 


if  (qb_niax  <  qb  )( 

max_length  =  length; 
qb_fnax  =  qb;> 

else 

Continue  =  No;  ) 

if  (nL»iiber_of_passes  >  2){ 
length  =  max_length; 
base_width  =~2  *  Ap  /  length; 
computeC ); 

eprintf(''\r\nUriting  Rec  #  %d  d  =  544. 3f  b  =  544. 6f  qb  =  544. 3f  Tb  *  544. If  Ap  =  544. If", 
j,  b-’se  width,  length,  qb,  Tbase,  Ap); 
out  =  fopen(*'OPTrRi.«Ar","a’'); 
edge  width  =  0; 

fprintf(out,  "54d,  54e,  54e,  54e,  %e,  %e,  %e,  54e,  %e,  54e,  Sle,  S4e\rSn",  j, 

depth,  length,  Ap,  k,  Tbase,  epsilon,  V,  qb,  taper_ratio,  base_width,  edge_width  ); 
fcloseCout ); ) 

else{ 

cprintf ("\r\nlgnoring  This  Set  Of  Parameters"); 

return  0;> 

void  computelvoidX 

double  m,  hr.  A,  B,  C,  D,  I1UA,  KOUB,  lOUB,  KIUA,  KOUA,  lOUA,  K1UB,  HUB; 

double  fin_taper,  b1,  db,  a,  ua,  ub,  a_width,  b_width; 

double  inc,  Tave,  overall_diff  =  100,  Tmax  =  0,~Tbv_nin  =  100000; 

double  Tav,  Ttip,  temp,  qa,  segment_diff;  ” 

int  overflow,  possible_overf low  =  No,  number_of_tries  =  0; 

edge_width  =  base_width  /  100; 

db  s~length  /  (  N~+  1  ); 

b1  =  length; 

length  =  length  -  db; 

fin_taper  =  atan2(  base_width,  2*bl  ); 

K2  =  2  •  epsilon  *  sigma  *  v  •  pow(Tbase,4); 

Tave  =  Tbase  •  pow(V,  0.25)  0.0001; 

if  (V==0) 

Ttip  =  .8625  •  Tbase  •  sqrt(  2  *  V  )  ♦  (  (  Tbase  -  (Tbase  •  sqrt(  2  •  V  )))/2); 
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else 

Ttip  =  Tave  ♦  0.0001; 
inc  =  .5  *  Ttip;  i=1; 

while  (fabs(overall_diff >  >  (  Tbase  *  error)  | |  i  <s  n  )( 

overflow  =  No;  i  =  1;  Ta  =  Ttip;  Tb  =  Ta; 
seginent_<liff  =  10; 

while  (fabs(se9iTient_di ff )  >  (  Tb  •  error)  /  N  £&  overflow  ==  No){ 

Tav  =  <Tb  ♦  Ta)  /  2; 

hr  =  sigma  •  epsilon  •  powlTav,  3)  -  (K2  /  (2  •  Tav)); 
m  =  hr  /  (k  •  sinCf in_taper)); 
ub  »  2  •  sqrt{  db  •  m”); 

lOUB  =  lO(ub); 

K1UB  =  Kl(ub); 

HUB  =  I1(ub); 

C  =  2*k*edge_width*  depth*pow(m,2)*I1UB/(ijb  •  lOUB); 


temp  =  Ta  /  lOUB; 
qb  =  C  *  temp; 

segment_diff  =  (  temp  -  Tb)  /  2; 

Tb  =  temp  ♦  segment_dif f ; 

if  (Tb  >  Tbase  *50  II  Tb  <  0){ 

overflow  =  Yes;  possibte_overf low  =  Yes; 
if  (Ttip  <  Tbv_min  )  Tbv_min  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /  2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave)  Ttip  *  Tave;)) 

if  (overflow  *=  No){ 

Ta  =  Tb;  qa  =  qb;  ♦♦i;) 

a  =  b1  -  length; 

while  (i  <*  N  &&  overflow  --  No){ 
segment_diff  =  10; 

a_width  s  edge_width  •  a  /  (  b1  -  length); 
b_width  =  edge^width  *  (  a  ♦  db  )  /  (  b1  •  length); 

while  (fabs(segment_diff )  >  (  Tb  •  error)  /  N  £&  overflow  No){ 

Tav  =  (Tb  ♦  Ta)  /  2; 

hr  =  sigma  *  epsilon  *  pow(Tav,  3)  •  (K2  /  (2  •  Tav)); 
m  =  hr  /  (k  *  sin(f in_taper)); 
ua  =  2  *  sqrtCa  *  ro  )7 
ub  =  2  *  sqrt(  (a  ♦  db)  •  m  ); 

JlUA  =  I1(ua); 

K0U8  =  KOlub); 
lOUB  =  lO(ub); 

K1UA  I  Kl(ua); 

KOUA  =  KO(ua); 
lOUA  =  lO(ua); 

K1UB  =  Kl(ub); 

HUB  =  Il(ub); 

A  =  ua  *  (nuA^kOuB  ♦  tOUB*K1UA); 

B  =  ub»ub*(IOUB»KOUA  -  I0UA*K0UB)  /  (  2«k*b  width*depth*tn); 
C  =  2*k*a  width*depth*m»ub*(llUB*KlUA  -  I 1UA»K1UB)/(ua); 

0  =  ub*(  T0UA»K1UB  ♦  nUB*K0UA); 

temp  =  A  •  Ta  ♦  B  •  qa; 
qb  =C*T8*0»  qa; 

8egment_diff  »  (  temp  -  Tb)  /  2; 

Tb  =  teiSp  ♦  segment_diff; 
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if  (Tb  >  Tbase  *50  [|  Tb  <  OX 

overflow  =  Yes;  possibte_overf low  =  Yes 
if  (Ttip  <  Tbv_mir\  )  Tbv_mTn  =  Ttip; 
inc  =  (Tbv_min  -  Tmax)  /~2; 

Ttip  =  fabs(Tbv_min  -  inc); 
if  (Ttip  <  Tave7  Ttip  =  Tave;» 

if  (overflow  ==  NoX 

Ta  =  Tb;  qa  =  qb;  a  =  a  ♦  db;  ++1;)) 


if  (nLfnber_of_tries  >  29X 

Ta=0;  Tb=0;  qb=0;  return;) 

if  (overflow  ==  NoX 

if  (Tmax  <  Ttip) 

if  (Tb  <  Tbase)  Tmax  =  Ttip; 

overall_diff  =  (Tb  -  Tbase)  /  2; 

if  (  possible_overf low  ==  YesX 
inc  =  inc  /  2; 

if  (overall_diff  >  0  )  Ttip  =  Ttip  -  inc 
else  Ttip  =  Ttip  ♦  inc;) 

else 

Ttip  =  Ttip  ■  overal t_diff; 
number_of  _t  r  i  es+*  ; 
if  (Ttip  <  Tave)  Ttip  =  Tave;)) 

length  =  length  ♦  db; 
return;) 
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APPENDIX  C  MATLAB  EVALUATION  M  FILES 


MATLAB  EVALUATION  H  FILE  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Rectangular  Profile 


load  rect; 

0  =  rect; 

N  =  100; 

sigma  =  S.66961e-8; 


L 

S 

D(:.2); 

b 

= 

DC 

.3); 

A 

S 

D(:,4); 

k 

= 

0( 

.5); 

To 

s 

D( 

.6); 

e 

= 

D( 

.7); 

V 

s 

0( 

.8); 

qo 

= 

0(:,9); 

TR 

= 

D( 

,10) 

db 

s 

D( 

.11) 

de 

s 

D( 

.12) 

I 

= 

[1:101] 

C3  =  A.*  k.*  ({(sigma*e).*2).»  To. ''9).*  (L./qo).* 
Cl  =  b.*  sigma.*  e.*  (To. *4).*  (L./qo); 

C2  =  db.*  k.*  sigma.*  e.*  (To. ‘5).*  (L./qo). '2; 


3; 


dispCRECT  V  =  O') 
dispC  ') 

disp( 'Constant  One:  Polynomial  p  is:') 

p  =  polyfit(TR,  Cl,  1) 

fiti  ‘  polyval{p,TR); 

dispCMean  Square  Error:') 

std_err  =  sqrt(sum((C1-fit1).'2)/(N-1)) 

pause 

dispC  ') 

disp( 'Constant  Two:  Polynomial  p  is;') 

p  =  polyfit(TR,  C2,  1) 

fit2  =  polyval(p,TR); 

dispCHean  Square  Error:') 

std_err  =  sqrt(sum((C2-fit2).'2)/(N-1)) 

pause 

dispC  ') 

disp( 'Constant  Three:  Polynomial  p  is:') 

p  =  polyf it(TR,  C3,  1) 

fit3  =  polyvalip.TR); 

dispCMean  Square  Error;') 

std_err  =  sqrt(sum((C3-fit3).*2)/(N-1)) 


¥ 
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/* 


MATLAB  EVALUATION  M  FILE  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Trapezoidal  Profile 
. . ! 


load  trap; 

0  =  trap; 

N  =  100; 

sigma  =  5.66961e-8; 

L  =  0(:,2); 
b  =  D(:.3); 

A  =  0(:,A); 
k  =  D(:.5); 

To  =  D(:,6); 
e  =  D(:,7); 

V  =  D(:.8); 
qo  =  0(:,9); 

TR  =  D(:,10); 
db  =  0(:,11); 
de  =  D(:,12); 

I  =  [1:1011; 

C3  =  A.*  k.»  <((sigma*e).*2).»  To. '9).*  (L./qo).'3; 
C1  =  b.»  sigma.*  e.*  (To. '4).*  (L./qo); 

C2  =  db.*  k.*  sigma.*  e.*  (To. *5).*  (L./qo).*2; 

dispCTRAP  V  =  O') 
dispC  '} 

disp( 'Constant  One:  Polynominal  p  is:') 

p  =  polyfit(TR,  Cl,  1) 

fiti  »  polyval(p,TR); 

dispCHean  Square  Error:') 

std_err  =  sqrt(sum((C1-fit1).^2)/(N-1)) 

pause 

dispC  ') 

disp( 'Constant  Two:  Polynominal  p  is:') 

p  =  polyfit(TR,  C2,  2) 

fit2  =  polyval(p,TR); 

dispCMean  Square  Error:') 

std_err  =  sqrt(sum((C2-f it2).'2)/(N-1)) 

pause 

dispC  ') 

disp( 'Constant  Three:  Polynominal  p  is:') 

p  =  polyfit(TR,  C3,  1) 

fit3  =  polyvaUp.TR); 

dispCMean  Square  Error:') 

std_err  =  sqrt(sum((C3-fit3).*2)/(N-1)) 

pause 


V 
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/* 


MATLAB  EVALUATION  M  FILE  To  Determine  The  Optimum  Geometry 
Longitudinal  Fin  of  Triangular  Profile 

**************************************************************************** 


load  tri; 

D  =  tri;  J 

N  =  100; 

sigma  =  5.66961e-8; 

L  =D(:.2); 
b  =  D(:.3); 

A  =  D(:,4); 
k  =  D(:.5); 

To  =  D(:.6); 
e  =  D(:,7); 

V  =  0(:,8); 
qo  =  D(:,9); 

TR  =  D(:,10); 
db  =  DC.ll); 
de  =  D(:,12); 

I  =  11:101]; 

C3  =  A.*  k.*  ((<sigma»e).*2).*  To.'9).»  {L./qo).'3; 

Cl  =  b.*  sigma.*  e.*  (To. “4).*  (L./qo); 

C2  =  db.*  k.*  sigma.*  e.*  (To. '5).*  (L./qo).'2; 

dispCTRI  V  =  O') 
dispC  ') 

disp( 'Constant  One;  Polynominal  p  is:') 

p  =  polyfit(TR,  Cl,  1) 

fiti  =  polyvaUp.TR); 

dispCHean  Square  Error:') 

std_err  =  sqrt(sum((C1-f itl).*2)/(N-1)) 

pause 

dispC  ') 

disp( 'Constant  Two:  Polynominal  p  is:') 

p  =  polyfit(TR,  C2,  1) 

fit2  =  polyvaKp.TR); 

dispCHean  Square  Error:') 

std_err  =  sqrt(sum( (C2-f i t2) . *2)/(N- 1 )) 

pause 

dispC  ') 

dispCConstant  Three:  Polynominal  p  is:') 

p  =  polyf i t(TR,  C3,  1 ) 

fit3  =  polyvaUp.TR); 

dispCHean  Square  Error:') 

std_err  =  sqrt(sum({C3-fit3).‘2)/(N-1)) 


» 
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